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ABSTRACT 



Previous investigations of longshore currents have included various simplifying as- 
sumptions and restrictions including a planar beach, a steady and depth uniform flow, 
spatiallv-invariant bed shear stress and turbulent momentum exchange, and the exclu- 
sion of surface wind stress. These assumptions are quantitatively investigated by calcu- 
lating the relative importance of each term in the longshore momentum balance with an 
emphasis on the relative importance of wind forcing across the barred nearshore. Wind 
and wave forcing of longshore currents across a barred beach are examined using both 
a numerical model and field measurements. A local momentum balance was measured 
at various locations across the surf zone during the SUPERDUCK experiment held at 
the U'.S. Army Corps of Engineers, CERC Field Research Facility. Duck. N.C.. U’.S.A. 
in October 19S6. A moveable sled was instrumented with pressure, current, and wind 
sensors to measure the various terms in the longshore momentum equation. Stability- 
dependent atmospheric drag coefficients for the surf zone are determined from wind 
stress measurements acquired just beyond the surf zone and wind speed measurements 
acquired from an anemometer atop the 9 m sled mast. Breaking waves were visually- 
identified and electronically marked on the data tapes. Mean current velocities at three 
levels were measured across the barred surf zone to examine the effects of a bar on the 
vertical structure of the How. A modified longshore current model is used to study the 
effects of barred topography. The predicted longshore current profiles are compared with 
field data and bed shear stress coefficients are determined from model-fitting. Significant 
spatial variability of longshore current [(7(1.0 ms -l )3 is observed across the bar with 
depth uniform longshore flow. Wind force is found to be a first-order term along with 
wave force under certain wind and wave conditions. Spatially-dependent bed shear 
stress coefficients relative to the bar are determined from a local momentum balance. 
Mean bed shear stress coefficients were 0.004 ± 0.0010, 0.002 + 0.0006, and 0.001 + 
0.0003 for offshore the bar, on top and immediately before the bar, and in the trough. 
Mean bed shear stress coefficients determined from model-fitting and from the momen- 
tum balance are in excellent agreement. Longshore current velocities are modeled within 
± 20% of observed velocities across the barred beach at SUPERDUCK with a mean 
least square error of 7-15%. 
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1. INTRODUCTION 



Previous investigations of longshore currents have utilized many simplifying re- 
strictions and assumptions. Major assumptions and restrictions include a planar beach, 
a steady and depth uniform longshore current, spatially-invariant bed shear stress and 
turbulent mixing, and the exclusion of surface wind stress. These assumptions are often 
too restrictive to permit realistic results. This research quantitatively investigates these 
assumptions by measuring the relative importance of each term in the longshore mo- 
mentum balance equation (see equation (1)). 

Modern theories for mean longshore currents are based on a longshore momentum 
balance (Bowen. 1969; Thornton, 1970; and Longuet-Higgins, 1970a, 1970b) with a ra- 
diation stress approach (Longuet-Higgins and Stewart, 1964). For a plane-parallel 
beach, the complete longshore momentum equation is: 
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where a right-handed Cartesian system is adopted with x increasing offshore (see 
Figure 1), and 

term 1 = temporal change of the time-averaged depth-integrated mean mo- 
mentum per unit area in the longshore direction due to both steady (M y ) 
and unsteady flow (My), 

term 2 = cross-shore gradient of the wave-induced longshore momentum flux 
(also called the wave-induced radiation stress gradient), 
term 3 = bottom shear stress in the longshore direction, 
term 4 = surface wind stress in the longshore direction, and 
term 5 = cross-shore gradient of the longshore momentum flux due to turbu- 
lence. 



Most investigators (Longuet-Higgins, 1970a; Jonsson et al., 1974; Kraus and Sasaki, 
1979; Thornton and Guza, 19S6; Losada et al., 1987; Baum and Basco, 1987; among 
others) apply the previously mentioned assumptions and restrictions and then reduce the 
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Wind Direction (0) and Wave Incident Angle (a) 



Figure 1. Section and plan definition for the SUPERDUCK e.xperiment (modi- 
fied from Birkemcier and Dalrymplc, 1975). 



equation for longshore momentum to a simplified balance between the wave-induced 
radiation stress gradient, the bottom shear stress, and the turbulence-induced radiation 
stress gradient. The most restrictive of these assumptions is that of a planar beach, since 
barred beaches are more common than planar beaches (Wright et al., 19S6). Temporal 
variability of longshore currents, often exceeding ± 100% of the mean flow, have been 
measured in the field (Dette and Fuhrboter, 1974; Wood and Meadows, 1975; Meadows, 
1976). The depiction of spatially- invariant bed shear stress and turbulent mixing across 
the dynamic surf zone, especially over a barred beach, may be unrealistic. Additionally, 
recent investigations addressing combined currents and waves indicate that the bed shear 
stress coefficient, which is directly proportional to bed shear stress itself, may be 1.5 to 
3.5 times larger in magnitude (Grant et ah, 1984) than the invariant values universally 
used by longshore current investigators. 

Longshore current forcing is studied utilizing the data sets acquired during the 
SL'PERDL'CK experiment. This experiment was held at the U.S. Army Corps of Engi- 
neers. Waterways Experiment Station's Coastal Engineering Research Center (CERC), 
Field Research Facility (FRF) in Duck. NC in September and October 19S6. A unique 
instrumentation system designed specifically for this experiment was utilized to measure 
a local momentum balance at various locations across the surf zone of this single-barred 
beach under various wave, current, and wind conditions. Instruments were mounted on 
a mobile sled to measure the various terms in the momentum balance of equation (1). 
Two triangular, differential pressure sensor slope arrays were spaced approximately 2 m 
apart. Coupled with an absolute pressure reference, the slope arrays can be used to 
measure the gradient of the wave-induced momentum flux. Marsh-McBirney 
electromagnetic (EM) current meters located at the centroid of the slope arrays can be 
used to determine both mean momentum as well as the total, or combined, momentum 
flux due to the wave-induced and turbulent components. The difference between the 
current meter determination of the total momentum flux ( S vx ) and the slope array de- 
termination of the wave-induced momentum flux ( S yx ) yields a measure of the depth- 
integrated Reynolds stress, or turbulent component of the momentum flux. The surface 
stress can be inferred from the wind velocity measured using a wind anemometer 
mounted atop the 9 m sled mast. 

Three methods for calculating the bed shear stress coefficient are attempted -- 
measuring the momentum balance, model-fitting to field data, and the mean profile 
method. By operating in areas where a plane-parallel beach assumption is valid and by 
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precisely orienting dynamic measurements perpendicular to the local bottom contours, 
all the terms contributing to the longshore momentum can be measured. Having meas- 
ured all the terms in the balance, the unspecified bed shear stress coefficient is then de- 
termined. Theoretical mean longshore current distributions generated by a numerical 
model are compared with field-measured longshore current speeds. The bed shear stress 
coefficients in the numerical model are then adjusted until the modeled and observed 
current distributions coincide in a least square sense. Preliminary indications from an 
on-going surf zone investigation suggest that the mean longshore velocity profile in the 
vertical appears to be logarithmic (personal communication. R. T. Guza, 1987). Hence, 
the bed shear stress coefficient might also be determined by the profile method of 
boundary layer physics. However, logarithmic profiles were not found in the SUPER- 
DUCK data. 

The sled was designed to make transects anywhere along the shoreline and as far 
offshore as the mast height would allow. Sled mobility was provided by U.S. Army 
Corps of Engineers equipment. The sled was tethered by a 300 m chain to a heavy-duty 
forklift on the beach. The sled was towed offshore by the Coastal Research Amphibious 
Buggy (CRAB) to the desired location beyond the surf zone and then detached from the 
sled. The forklift dragged the sled onshore along its pre-planned transect of data col- 
lection locations. 

A previous difficulty with directly measuring momentum flux in the field was the 
precise orientation required of tire instrumentation (Thornton and Guza, 1981). In the 
SUPERDUCK experiment, precise orientation of the sled (and attached slope arrays) 
was achieved to within 0.5° using laser ranging to triangulate on two prisms mounted 
2.4 m apart on a mast spreader above the sled. A daily calibration of current meter 
orientation was accomplished by initially positioning the sled beyond the surf zone for 
the first data run of a surf zone transect. Based on conservation of radiation stress 
outside the surf zone, the current meters were numerically oriented so that their respec- 
tive radiation stress values were identical. This orientation was then held constant for 
the remaining sled locations along that transect. This calibration procedure is essential 
to obtaining realistic physical results from the data and is a distinguishing feature of the 
research. 

A primary improvement between the SUPERDU'CK experiment and previous 
large-scale nearshore experiments was the extensive bathymetric data. Daily 
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bathymetric profiles of the entire experiment area were accurately measured with the 
CRAB and the laser ranging system. 

Pertinent literature is reviewed in the following chapter. A numerical model is de- 
veloped in Chapter 111 to compare wind and wave forces for various wind speeds and 
directions. The SUPERDUCK field experiment is described in Chapter IV. Data anal- 
ysis of the field experiment is extensive and therefore is addressed in separate chapters. 
Preliminary data analysis is described in Chapter V and includes an overview of the ob- 
served mean current structure over the bar. Since surface wind stess is of major impor- 
tance to this research and few wind stress studies have been conducted in the surf zone. 
Chapter VI is entirely devoted to the parameterization of the surf zone surface wind 
stress for this field experiment. Calculation of the longshore momentum balance, de- 
termination of the relative contributions of each term in the balance, and calculation of 
spatially- variable bed shear stress coefficients are addressed in Chapter YU. A longshore 
current model for planar beaches is modified to account for barred topography and is 
applied to the data. The model's longshore current results are compared to observations 
and bed shear stress coefficients are determined by model-fitting in Chapter Y1I1. 
Summary and conclusions are presented in Chapter IX. Because of the large quantity 
of data presented graphically for each day of data analysis, only a single typical day 
(October 16) is illustrated within the text. Data analysis figures for other days are ref- 
erenced in the text by figure number but are collectively grouped in Appendix A. 

The data obtained from this field experiment represents the first quantitative data 
set, including directional wave, surface wind stress, bathymetry, and breaking wave 
measurements, for the study of longshore current over barred beaches. In most previous 
studies, wind forcing of longshore currents has not been considered due to the absence 
of quantitative data. This research quantitatively determines the role of wind forcing in 
driving longshore currents from both field measurements and a simple numerical model. 
Wind forcing, as well as a variable bed shear stress coefficient, can represent measurable 
and reasonably significant terms in the momentum balance equation, and therefore must 
be taken into account in any quantification of longshore currents. 
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II. LITERATURE REVIEW 



A. INTRODUCTION 

One of the earliest works on coastal hydrodynamics and processes was that of 
Johnson (1919). He qualitatively discussed wave set-up, longshore currents, and rip 
currents (undertow) as hydraulic currents due to waves. He utilized incident wave angle 
as his criterion for rip currents (waves approaching perpendicular to shore) and for 
longshore currents (waves approaching oblique to shore). 

In World War II, amphibious landing craft had difficulty landing their troops in the 
presence of a longshore current. The current often represented a greater hazard than the 
breakers themselves by causing the craft to broach inside the breaker zone. The U.S. 
Navy's concern for this effect prompted Putnam et al. (1949) to conduct one of the first 
applied research studies on longshore currents. 

Since that time, various theories have been postulated to account for the basic 
physics of longshore currents and their predicted velocity distribution across the surf 
zone. These theories include conservation of energy, conservation of mass, empirical 
formulation, and conservation of momentum . 

B. RADIATION STRESS 

Modern longshore current theory is based on conservation of momentum. One of 
Putnam et al.'s (1949) approaches was the application of conservation of momentum. 
The momentum flux of an obliquely incident wave was averaged over a wave period to 
give a mean momentum flux entering the surf zone. The longshore component of the 
change (due to wave breaking) in this momentum flux provided a driving force for 
longshore currents. They assumed a steady-state longshore current, and therefore the 
driving force was balanced by a bottom frictional force. 

The momentum flux approach was formalized and the physics clarified when the 
concept of radiation stress was introduced. Longuet-Higgins and Stewart (1964) defined 
radiation stress as the "excess flow of momentum due to the presence of the waves." 
Elaborations on the concept of radiation stress are contained in Longuet-Higgins and 
Stewart (1960, 1961, 1962, 1963, 1964), Bowen (1969), Longuet-Higgins ( 1970a, 1970b, 
1972a, 1972b), and Thornton (1970). 

Theory was further formalized with the derivation of time-averaged and depth- 
integrated conservation of mass and conservation of momentum which separates the 
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mean and fluctuating momentum contributions {Phillips, 1977; Mci, 19S3). This ap- 
proach has been used extensively in the study of wave set-up and set-down, horizontal 
longshore current profiles, rip currents, and surf beat generation. This approach forms 
the foundation of modern longshore current theory, as well as the foundation for the 
momentum balance of equation (1) on which this dissertation is based. 

When there is an oblique wave approach to the shoreline, the longshore component 
of radiation stress (or longshore momentum flux due to unsteady motion) is given by 
(Longuet-Higgins, 1970a): 



u'.v' = fluctuating (x. y) components of the flow velocity, and 

the overbar denotes a time-averaging operator. 

The basic concept is that the longshore current velocity is related to the change in in- 
coming wave momentum, which in turn relates to the radiation stress. This change of 
momentum flux in the cross-shore direction is proportional to the rate of wave energy 
dissipation, where dissipation is due primarily to surface breaking but can include bot- 
tom, and internal (turbulent) dissipation. By using linear wave theory, radiation stress 
for monochromatic waves can be related to the incident wave angle (a) at the measure- 
ment location by: 

S yx =En sin a cos a (3) 

where E is total wave energy density, a is the wave incident angle, or angle between the 
wave crest and a line parallel to the shoreline, and n is the ratio of wave group velocity 
to wave phase velocity expressed in terms of wave number (k) as: 




( 2 ) 



where ?; = the sea surface. 

h = water depth as measured from the still water line (SWL). 
p = water density, 




( 4 ) 
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Radiation stress has also been applied to random waves (Collins, 1972; Battjes, 1972a; 
Thornton and Guza, 19S6) and non-linear wave theory (James, 1974a, 1974b) to solve 
for the longshore current. 

The conservation laws of mass, horizontal momentum flux, and radiation stress by 
Phillips (1977) and Mei (19S3) are both widely applied and similar. Both develop the 
equations for mean water motion by integrating over depth and averaging over a wave 
period. They differ in their approach to the mean flow due to unsteady motion (i.e. 
waves). Separating the mean and fluctuating horizontal velocity components (w z ) as: 

u t = I) + u t i= 1 , 2 , (5) 

Mei(19S3) defines the mean flow as: 

b'/jrj-.O = - - ■■■■-■ f u t dz i— 1,2 (6) 

}] + h i~ h 



and 



n U: dz = 0 / = 1.2 (7) 

J— ft ' 

Thus he incorporates the mean flow due to unsteady motion in the mean velocity term, 
U h resulting in a mean momentum per unit area defined as: 

m. = JJ. = pu-ti + h) (S) 

Phillips (1977) utilizes equation (5) and (6), however he does not include the mean 
flow due to unsteady motion in the term U t . He evaluates the time-averaged depth- 
integrated current flow by partitioning the integration of u (equation 7) into two sepa- 
rate regions: 



Uj dz = 
-h ‘ 



’° ^ 

U: dz + Uj dz 

-A J° 



( 9 ) 



where the first term on the right-hand side is zero by definition. This results in a mean 
momentum per unit area of 
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i = 1,2 



(10) 



Mi = p L) (>i + h) + p Uj dz 

J° 

Phillips approach is adopted for this research. 

Since the momentum balance of equation (1) represents the basis for this research, 
each remaining term -- surface wind stress, bottom shear stress, turbulent momentum 
flux, and temporal variability — will now be examined. 

C. SURFACE WIND STRESS 

Surface wind stress is a shear stress imparted to the water column when the wind 
flows across the water surface. For simplicity, surface wind stress has been historically 
neglected as a forcing function for longshore currents. Although no physically-based 
model has specifically evaluated the wind stress contribution to longshore currents, there 
is empirical evidence suggesting wind stress can be important, and therefore should be 
included as a driving force for longshore currents. Nummedal and Finley (197S) used a 
stepwise multiple regression procedure utilizing surf zone environmental parameters to 
explain the variability of longshore current velocities. Although this procedure was 
empirical and not theoretical, the longshore current component of the wind velocity was 
found to account for more of the observed current variance than any single parameter, 
or combination of parameters, descriptive of the breakers. 

Fox and Davis (1973) were able to relate longshore currents with barometric pres- 
sure at Lake Michigan. Although wind is certainly related to barometric pressure, they 
did not expand their work to directly deal with wind. Komar (1976) mentioned that 
winds blowing in the longshore direction may contribute to the observed flow. 

Birkemeier and Dalrymple (1975) are one of the few nearshore modelers who did 
incorporate a wind forcing term (albeit, an underweighted term) in their longshore cur- 
rent model. Notwithstanding the diflerent spatial scales involved, it is noted that lake 
(e.g., Sheng et al. , 1978) and offshore modelers do consider wind stress whereas 
nearshore modelers do not. Nearshore modelers either assume that wave stress domi- 
nates over surface wind stress or that no wind is present. 

Hubertz (19S6) calculated the spectra of the mean longshore wind component and 
mean longshore current for a four-month period at the FRF in Duck, N.C. Longshore 
current spectral peaks were present at periods of 4.6 days, 23.8 hours, and 12.5 hours. 
Longshore wind component spectral peaks were present at 5.9 days and 24.4 hours. He 
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related the two longest period current peaks to wind forcing by synoptic weather systems 
and by sea breeze conditions, respectively. The 12.5 hours peak was attributed to tides. 

There arc several methods for calculating surface wind stress on the water surface. 
These methods are: 

• the direct method (also called the covariance or eddy correlation method), 

• the profile method, 

• the dissipation method, and 

• the drag coefficient method (also called the bulk aerodynamic method). 

Each method has its own inherent advantages, disadvantages, and accuracy. The drag 
coefficient method was selected for this research based on equipment availability, the 
rough environment and platform involved (surf zone measurements by a mobile sled), 
and the availability of wind stress measurements acquired by other field investigators at 
the SLPERDUCK experiment. 

The drag coefficient method formulation for surface wind stress (Taylor, 1916) in the 
longshore direction is: 

r; = p a C a \llVV y ( 11 ) 

where tJ? = time-averaged surface wind stress vector in the longshore (y) direction 
(hereafter referred to as surface wind stress), 
p a = atmospheric density, 

Q = an atmospheric drag coefficient (generally of 0(1Q~ 3 ) for a height of 10 
m above the water surface) for the specific height above the water surface 
at which W is measured, 

W = the total wind speed, and 

JV y = the longshore component of the wind velocity. 

The absolute value is used in equation (1 1) to retain directionality. 

A large quantity of research has been directed at determining values for C d for land 
and the open ocean. Garratt (1977) provides an excellent review for investigations of 
Q for land and the open ocean prior to 1977. More recent investigations include Wu 
(19S0, 1982), Safaie (1984), Geernaert and Katsaros (19S6), Geernaert et al. (1986, 
19S7a, 19S7b), and Hsu (1986a, 1986b, 19S7). Unfortunately, few investigations have 
dealt with determining Q near the land-sea interface and in the surf zone. C d magnitudes 
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of 0(1.7 x 10 -3 to 2.4 x 10“ 3 ) were determined for nearly calm conditions outside the 
surf zone of a Chukchi Sea barrier island (Walters, 1975) and C d magnitudes of 0( 
2.0 x 10 -3 to 7.0 x 10 -3 ) for the upper swash zone of a Florida beach (Hsu. 1970). 
Walter's (1975) investigation was conducted under adiabatic and onshore wind condi- 
tions. whereas Hsu's (1970) investigation included adiabatic as well as non-adiabatic 
conditions with onshore wind conditions. 

Hsu (1973) determined surf zone wind stress values for onshore winds of 0.5 to 8.5 
ms -1 under adiabatic conditions for a two-week period at Santa Rosa Island, FL. C d 
magnitudes of 0(1.0 x 10 -2 to 1.5 x 10~ 3 ) are obtained upon conversion of his surface 
wind stress values. No breaker wave height values were indicated; however, it is assumed 
that they were small due to the experiment's location (northwest Florida coast), month 
(June), wind speeds, and a photograph accompanying the published results. 

Yugts and Cannemeijcr (1981) calculated C d for non-adiabatic and onshore wind 
conditions from Field measurements at a Danish beach. They instrumented two 20 m 
masts with multi-level cup anemometers and temperature sensors and placed one at the 
high-tide line and one directly inland at a distance of 127 m. They concluded from ap- 
proximately 100 half-hour profiles that the drag coefficient during unstable conditions 
over the surf zone is the same as over the ocean. They calculated Q magnitudes of 0( 
1.3 x 10 -3 to 1.5 x 1 0 — 3 ) for wind speeds at 10 m height from 3-10 ms -1 . 

Geernaert ct al. ( 1 9 S 7 ) measured wind stress using a sonic anemometer from a re- 
search platform in the North Sea. They found an empirical relation between the neutral 
drag coefficient (Q„) and the ratio c 0 /u>, where c 0 is the dominant long wave phase speed 
and u. is the surface wind friction velocity. They applied this relationship to other data 
sets taken in varying water depths. They found increased C dn values for shallower depths 
(Figure 2). As wind speed increased to 25 ms -1 , the C dn value for 10 m depth was in- 
creased by a factor of 2 (to 0(3.0 x 10 -3 )) over that for deep water. 

D. BOTTOM SHEAR STRESS 

Putnam and Johson (1949) conducted one of the earliest investigations of nearshore 
bottom stress or bottom friction. They assumed sinusoidal oscillatory bottom motion, 
an invariant friction coefficient, negligible perpendicular flow at the bottom due to 
percolation, a planar bottom slope, and a constant water density. Based on their inves- 
tigations, bottom shear stress (r*) may be represented by: 

z b = pc f \ u t | u, (12) 



11 




Figure 2. Calculated dependence of C ^ on the 10 in height nind speed for five 
depths (from Geerneart et al., 19S7b), 

where Cj = the dimensionless bed shear stress coefficient due to waves and currents 
and generally of 0(10"-), and 

u t = resultant vector sum of the mean longshore current velocity (V) and the 
instantaneous wave orbital speed near the bottom (see Figure 3). 

As with wind stress, the absolute value is used in equation (12) to retain directionality. 
This equation has essentially formed the basis for modern bottom friction studies. In 
general, later investigators have worked toward its modification as well as solving for 
empirical values for cj. 

One can decompose u and v into a steady or mean velocity (U, V) and an unsteady 
or fluctuating velocity (u\ v'). The unsteady velocity may be further separated into 
wave-induced (u, r) and turbulence-induced (u", v") components (Mei, 19S3): 

U=lj J ru'=lj J rU J rU (13) 

V = V + v' = V + 7 + v" (14) 



12 



Incorporating the above decomposition, the bottom shear stress is represented as a 
vector oscillating in both direction and magnitude due to the interaction of a quasi- 
steady current with superimposed wave motion. From conservation of mass flux and the 
assumption of straight-and-parallel bathymetric contours, the mean cross-shore current 
velocity (U) is assumed zero. u t is then expressed as: 

u t — u x i + (V+ u y ]j = ui + vj = u + v (1 5) 

where u x = u w cos a, 
u y = u w sin a, 

u w = wave orbital speed just above the bottom and immediately outside the 
bottom boundary layer, and 

A 

i , &ij? unit vectors in the x and y directions, respectively. 



Bottom shear stress (t 6 ) can be written as: 



0 V • , o • 

T =T X I + T J 



(16) 



and from equations (12) and (15): 

r b = pry | u t \ u + pcy | u. | v (17) 

The longshore component of time-averaged bottom shear stress can then be written as: 

(18) 



T y = pCf\ U, | V = pCj{u 2 + V 2 ) 2 V 



The time-averaging procedure must be performed after evaluating the velocity compo- 
nents. 

The bed shear stress coefficient is a function of bottom roughness and flow intensity- 
near the bed. Bottom roughness is a function of both sand characteristics and bed form, 
where density and grain size define sand characteristics and bed form is described by the 
existence of ripples, which are in turn, a function of the flow field. Hence, the bed shear 
stress coefficient, bottom roughness, and flow intensity are all interrelated. 

Often in geophysical research, laboratory results complement field measurements 
and assist in building a particular data base as well as increasing our knowledge of the 
phenomena being studied. Unfortunately, very- few of the large number of laboratory- 
studies of ty are appropriate for extrapolation to the field. The flow field in nature will 
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Figure 3. Vector Relationship of 7i, due to presence of naves and currents. 
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almost always be in the rough turbulent flow regime (Jonsson, 1966). Since the 
Reynolds number is the appropriate dynamical similitude scale parameter, large 
Reynolds numbers are required to be in this flow regime. Hence, only laboratory studies 
performed in large wave flumes or oscillating water tunnels, where longer wave periods 
can be obtained, appear appropriate. 

Numerous investigators have solved for Cy under waves (Table 1) with results rang- 
ing from O(10 -1 to 10“ 3 ) for different beaches. Most have sought an invariant value for 
Cf based on a given set of conditions (e.g., bottom roughness, beach slope, or for only 
one particular beach), without regard for any spatial or temporal variability. Vitale 
(1979) states that a common engineering practice is to assume Cy = 0.01. Johns and 
Jefferson (1980) showed the variation of bottom stress within a wave cycle in a numer- 
ical analysis. They found increasing cy with increasing bottom roughness and decreasing 
wave period. 

There are four methods of bottom shear stress, and subsequent cy , determination: 
- momentum balance, model-fitting, the profile method, and direct shear measurements. 
The first three are indirect methods and are addressed below. The direct shear meas- 
urement method involves directly measuring the shear force, such as with pressure plates 
fixed on the bottom or measuring the Reynolds stress at the bed. This has not been 
done in a surf zone and is therefore omitted from further discussion. 

1. Momentum Balance 

Forcing terms and assumptions are specified (often based on data availability) 
to form an appropriate longshore momentum balance equation. The bottom stress term 
and its accompanying bed shear stress coefficient are then determined as residual values. 
The strength of this method depends on how complete the momentum balance is and 
how realistic are the restrictions and assumptions. Thornton and Guza (19S1) utilized 
this method at Torrey Pines Beach, CA (a near-planar, unbarred beach) in connection 
with the Nearshore Sediment Transport Study (NSTS) experiments. Assuming steady- 
state conditions, neglecting wind forcing, and operating in an area of straight-and- 
parallel contours, they assumed a simple balance where the cross-shore gradient of the 
total longshore momentum flux (both wave-induced and turbulent components) was 
balanced by the mean longshore bottom shear stress. For simplicity, they neglected the 
temporal variation term. They found considerable temporal and spatial variation of Cy 
with a mean value of 0.01 + 0.01 . However they stated that the results were to be used 
with caution due to the difficulties in directly measuring the radiation stress (S yx ) and its 



15 



Table 1. BED SHEAR STRESS COEFFICIENTS PROPOSED AND UTILIZED 



BY VARIOUS INVESTIGATORS 



Investigator 


c f 


Comments 


Bretschneider and Reid (1954) 


0.01 




Kishi 1 (1954) 


0.03 - 0.04 




Iwagaki and Tsuchiya 1 (1966) 


0.01 - 0.04 




Iwagaki and Kakinuma (1967) 


0.01 - 1.16 


various data sets 


Jonsson (1966) 


0.01 




Hasselman and Collins (1968) 


0.015 




Longuet-Higgins (1970a) 


0.01 




James 2 (1972) 


0.0025 - 0.001 




Jonsson et al. ( 1974) 


0.05 - 0.10 




Kamphius (1975) 


0.004 - 0.250 


f(bottom roughness) 


Komar 2 ( 1975, 1976) 


0.017 + 0.0225 




Huntley ( 1 976) 


0.0026 + 0.0006 


non-linear eqn. 


Liu and Dalrymple ( 197S) 


0.04 


Hsiao and Shemdin (1978) 


0.002 - 0.100 


various data sets 


Madsen et al. 2 ( 1 9 7 S ) 


0.01 




Kraus and Sasaki (1979) 


0.011 - 0.024 




Vitale (1979) 


0.003 - 0.250 




Johns and Jefferson (1980) 


0.0096 - 0.00319 




Grant et al. (1984) 


0.005 - 0.011 




Wind and Vreugdenhil (19S6) 


0.0015 




Guza et al. (1986b) 


o.oos 




Thornton and Guza (1986) 


0.006 + 0.0007 




• if required, the investigator's friction factor was converted to a bed shear stress 


coefficient from c,-— (friction 

1 from Sonu (1975) 

2 from Basco (1982) 


factor 2). 





gradients. These difficulties were due to angularity errors caused by lack of instrument 
resolution, installation misalignment, and errors associated with choosing the correct 
alongshore direction. All of these errors prevented accurate directional resolution of 
waves with small angles of incidence (< 3°). Resolution of these difficulties was ad- 
dressed in the SUPERDUCK experiment by utilizing a precise surveying system and 
selecting days for analysis when the mean wave incident angle was large (> 8°). In the 
following determination of Cj by the momentum balance method, both the temporal 
variability and wind stress terms in equation (1) are included. 
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2. Model-fitting 

Theoretical horizontal mean longshore current distributions generated by a nu- 
merical model are compared with field-measured longshore current speeds. The cj in the 
numerical model is then adjusted so that the two profiles agree. The resulting numerical 
model values for Cj can then be assumed to apply for the field measurements. Selecting 
the appropriate physics for the numerical model is obviously important for this method. 
Thornton and Guza (19S6) used this method with a narrow-band, random-wave trans- 
formation model for data collected in another NSTS experiment at Leadbetter Beach. 
CA. A hierarchy of Cy solutions were derived for linear and non-linear formulations for 
bottom shear stress, and with and without turbulent momentum exchange. Thornton 
and Guza (19S6) found a least-squares "best fit" value of cy= 0.006 + 0.0007 using the 
full non-linear, and presumably more correct, formulation for the bed shear stress. 

3. Profile Method 

The profile method (e.g., see Roll, 1965; Schlichting, 1968) has been utilized in 
boundary layer physics for years to determine friction or shear velocities and roughness 
heights. It assumes a constant stress layer with a mean logarithmic vertical flow struc- 
ture and requires at least two simultaneous speed values at different heights above the 
same location. 

Grant et al. (19S4) conducted high-quality near-bottom boundary layer meas- 
urements at a continental shelf location in 90 m water depth. Even at this substantial 
depth, they found the near-bottom flow field to be composed of oscillatory currents due 
to swell superimposed on mean currents, a situation similar to Figure 3. The mean, 
near-bottom (<2 m) velocity profiles were found to be logarithmic. Using the profile 
method, mean stress values were 3 to 7 times larger than expected from predictions of 
mean current alone using a typical smooth-bottom drag coefficient, and 1.5 to 3.5 times 
larger than expected for predictions using a drag coefficient based on the observed rough 
bottom. This significant increase in mean stress magnitudes was attributed to the pres- 
ence of the swell wave oscillatory currents. These results may apply to the surf zone, if 
the surf zone's entire vertical extent of chaotic flow may be likened to a fully developed, 
neutral, turbulent boundary' layer in the longshore direction. 

Analogous to surface wind stress (e.g., see Roll, 1965; Schlichting, 196S), the 
longshore bottom shear stress may be written as: 

T b y = p\\ (19) 
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where v» is the longshore shear or friction velocity associated with the mean How. 
Combining equations (IS) and (19): 



i 



2 

r. 



( 20 ) 



(UJb) ( W 2 + V 2 )ir 



The time-averaged velocity profile in a neutral turbulent boundary layer over a 
topographically simple bottom is logarithmic and may be expressed as: 



where k is von Karman's constant and and is the roughness height. By simultaneously 
measuring current speed at two different heights (V(z t ) and V(z 2 )), two equations in the 
form of (21) may be written. The two unknowns, v, and z 0 , are first solved, and then Cy 
is obtained from equation (20). Grant et al. (1984) employed a similar profile method 
in their deep water (90 m depth) investigation where friction velocity was measured in a 
downstream direction and u and v were measured in the streamwise and cross-stream 
directions, respectively. 

E. TURBULENT MOMENTUM FLUX 

The horizontal transfer of momentum due to turbulent mixing processes is not well 
understood. When the classical analytical solutions for the horizontal profile of the 
longshore current were initially solved using monochromatic breaking waves (Bowen, 
1969; Longuet-1 liggins, 1970a, 1970b; Thornton, 1970), a velocity maximum and dis- 
continuity was observed at the breaker line, with zero longshore velocity at the shoreline 
as well as beyond the breaker line. This occurred since there were no driving forces 
outside the sharply defined breaker line to generate a longshore current and there was 
maximum forcing occurring at the breaker line. However, field and laboratory data 
showed a smooth transition of the longshore current from outside the surf zone to a 
maximum inside the breaker line. This observation was attributed to turbulent mixing 
(lateral shear stress). Analytical solutions for longshore currents on a planar beach were 
proposed with the assumptions of isotropic turbulence, a planar beach, monochromatic 
waves, and the parameterization of a lateral shear stress in terms of an eddy viscosity 
coefficient multiplied by a horizontal velocity gradient. Values for the eddy viscosity 
coefficient were obtained by model-fitting. 




IS 



Whereas the bed shear stress coefficient is a physical quantity related to bottom 
roughness and flow intensity, the eddy viscosity coefficient is a function of the flow field 
and is a parameterization used to remove the longshore current discontinuity at the 
breaker line. An alternative approach to a monochromatic wave model is to describe 
breaking wave heights in terms of a probability distribution. This results in a more re- 
alistic region of wave breaking rather than a distinct breaker line. This spreading out 
of the breaker line results in a smooth horizontal profile of the longshore current, thus 
eliminating the need to include an eddy viscosity (Collins, 1970; Battjes, 1972a). 

Thornton and Guza (19S6) utilized a random w'ave model and model-fitting to de- 
termine values of the bed shear stress coefficient and eddy viscosity coefficient for both 
a linear and non-linear expression for the bed shear stress. Eddy viscosity coefficients ( 
C) were represented by an adjustable coefficient, X, in the range of 0 < X < 0.016 , 
expressed by (Longuet-Higgins, 1970a, 1970b): 



where g is the acceleration due to gravity. Optimal cy and X coefficients were determined 
by least-squares fitting between calculated and measured longshore current speeds 
shoreward of the mean breaker line. Using a non-linear bed shear stress formulation, a 
Cj value of 0.006 was determined. Eddy viscosity was found not to be important for the 
near planar topography of their field experiment. 

The last term in equation (1) represents an expression for horizontal diffusion or 
lateral mixing of turbulent momentum in terms of a depth-integrated Reynolds stress. 
The longshore component of momentum flux due to unsteady motion (equation 1) is the 
depth-integrated covariance between the fluctuating velocity components of u and v, 
measured by the current meters. By partitioning the unsteady motion as that due to 
waves and turbulence (equations 13 and 14), the excess momentum flux (equation 2) can 
be separated into a wave-induced (5 >JC ) and turbulent component (S" yx ): 



It has been assumed for this separation that the characteristic time scales of the two 
fluctuations are vastly different (i.e., u i and u"i are uncorrelated and u"j = 0). inte- 
grated viscous stresses are negligible, and that the horizontal scales of depth and mean 
velocity are comparable. The turbulent component of the excess momentum flux due 



t = A r I .v I (gh) 2 .v > 0, 



( 22 ) 




(23) 
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to unsteady motion may be obtained by comparing the total radiation stress (5 >IX ) 
measured by a current meter to the wave-induced radiation stress (S v;c ) measured by a 
slope array. 

F. TEMPORAL VARIABILITY 

Most previous investigators have assumed steady-state conditions for longshore 
momentum. However, significant temporal variations in longshore current velocity 
measurements have been observed (Putnam et al., 1949; Inman and Quinn. 1951; Inman 
and Bagnold, 1963; Dette and Fuhrboter, 1974; Dette, 1974; Huntley and Bowen, 1975; 
Wood and Meadows. 1975; Meadows, 1976, 1978; Guza and Thornton, 1978; Guza et 
ah, 19S6b; and others). The temporal variability of mean currents (dM y /dt) will be in- 
cluded in the analysis where M y is the time-averaged and depth-integrated mass flux per 
unit area, or momentum per unit area, in the longshore direction due to steady (M y ) and 
unsteady (M y ’) flow; 



My = My + My' = P V( V + II) + 



"i 

pv'dz 

-h r 



(24) 



Inman and Bagnold (1963) ascribed temporal variability of longshore currents to the 
wave groupness. Meadows (1976) concluded that the total observed longshore current 
velocity was composed of three components; 

V — 1 steady "f ^ wave "h ^ long (-5) 



where F steady is the steady longshore component, ave is the coincident wave period 
fluctuating component, and I 7 lon£ is the long period (as related to surf beat or low mode 
edge waves) fluctuating component. He defined the steady state component as the ve- 
locity component whose period of variation was greater than 200 seconds. The selection 
of 200 seconds was not explained. Guza and Thornton (197S) concluded that a temporal 
averaging period for mean longshore currents was not well-understood. Various re- 
searchers utilize different averaging times to determine the steady state component, 
based more on record length or data analysis, than on rationale. 

G. VERTICAL DISTRIBUTION OF MEAN LONGSHORE FLOW 

The equations by Phillips (1977) and Mei (1983) integrate the water motion over 
depth. For simplicity, the longshore current is assumed uniform from the surface to the 
bottom. No longshore current theories account for a depth-dependent longshore cur- 
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rent. Theoretical work on the depth-dependency of longshore currents (Stive and de 
Vriend, 1987) and cross-shore currents (Svendsen et ah, 1987) has only recently been 
undertaken. Most laboratory and field investigations have not measured the depth- 
dependence of longshore currents, and as a result, there is little field data on vertical 
current structure. 

Greenwood and Sherman (1983) measured vertical structure of the longshore cur- 
rent at the Wendake Beach, Ontario field experiment. They calculated mean flows for 
three elevations as well as a depth-averaged mean longshore current velocity. Their 
maximum mean variability was only 10% of the depth-averaged flow. Thus their data, 
as well as Meadows' (1978) and Hubertz's (1986) data, suggests that modeling longshore 
currents as uniform over depth is appropriate. Although they assumed a logarithmic 
boundary layer velocity profile in order to calculate a boundary roughness length (z 0 ), 
their logarithmic plots of only two points arc inconclusive regarding the existence of a 
logarithmic vertical velocity profile in the surf zone. However, their linear plots of ver- 
tical profiles show a "consistent deformation of the velocity field strongly suggestive of 
the influence of bed friction," and hence a logarithmic profile within at least 0.5 - 1.0 m 
of the bed remains a possibility for their data. 

Wright et al. (1986) obtained field data of vertical profiles from 0.1 to 1.0 m above 
the bottom of longshore currents in a bar trough. Although the longshore velocities 
were on the same order of magnitude as the wave orbital velocities, the bottom bound- 
ary layer was dominated by the waves acting in the shore-normal direction, orthogonal 
to the longshore currents. The velocity profiles were not logarithmic and it appeared 
that they were strongly influenced by wave-generated vortices from the bed ripples 
present. 

Theoretical investigations by Lundgren (1972), Smith (1977). and Grant and 
Madsen (1979), laboratory studies by Bakker and Van Dorn (1978) and Grant and 
Madsen (19S2), and field measurements by Grant et al. (1984) assume that the mean 
velocity profile is comprised of an inner region (close to the bed which is directly affected 
by wave-induced turbulence) and an outer region above it, and that the outer region has 
a logarithmic profile which is also affected by wave-induced turbulence. Preliminary 
indications from an on-going surf zone field investigation suggest that the mean 
longshore velocity profile in the vertical might be logarithmic (personal communication, 
R.T. Guza, 1987). The existence of a logarithmic profile would allow bed shear stress 
coefficients to be determined by the profile method. 
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H. BARRED BEACHES 

Although barred beaches are more common (Wright ct al., 19S6), far less is known 
about surf zone currents on barred beaches than on planar beaches. Allender and 
Ditmars (19S1) found little spatial variability in the horizontal longshore current profile 
across the entire surf zone of a steep, barred beach during storm conditions on Lake 
Michigan. No significant cross-shore flow or vertical variability of current flow was 
observed. Greenwood and Sherman (1983, 1986) studied surf zone currents in Lake 
Huron off a steep, multi-barred beach at Georgian Bay. They found significant cross- 
shore variability in the longshore current, and up to a four-fold variation during storm 
conditions. These two contrasting conclusions are indicative of the embryonic state of 
knowledge concerning the surf zone dynamics of barred beaches. The SL'PERDUCK 
experiment was conducted on a barred beach. Data analysis of the surf zone currents 
will be used to improve the understanding of surf zone dynamics at a barred beach. 



III. WIND FORCING OF LONGSHORE CURRENTS FROM A 
FINITE-DEPTH SPECTRAL WAVE MODEL APPROACH 



To assess the importance of the wind contribution to longshore currents, wave and 
wind forcing terms in the longshore momentum equation are integrated across the surf 
zone in a simple model and compared. Wave and wind forcing are coupled by assuming 
the waves driving the longshore currents are generated by the local wind. The TMA 
wind-wave shallow water wave spectrum is used as input to the radiation stress. 

A. THE TMA SPECTRUM 

The TMA spectrum (Bouws et al., 1985) was developed from data obtained from 
two extensive shallow-water wave measurement programs. MARSEN (Marine Remote 
Sensing Experiment at the North Sea) and ARSLOE (Atlantic Ocean Remote Sensing 
Land-Ocean Experiment), and data from a prolonged storm near lightship TEXEL in the 
southern North Sea. The MARSEN and ARSLOE sites were both on the continental 
shelf with depths less than 40 m. The combined data set was named TMA. after the 3 
sub-data sets, and includes water depths from 6-42 m, wind speeds up to 30 ms -1 , 
bottom composition from fine to coarse sands, and beach slopes ranging from 1:150 to 
nearly flat. The uniqueness of the TMA data set to this research is that the ARSLOE 
data, which comprises two-thirds of the nearly 3000 wind-wave spectra, were measured 
at the CERC Field Research Facility in Duck, N.C. 

Bouws et al. (1985) proposed, and then verified, that a finite-depth wind-wave 
spectrum could be approached through application of the similarity principles by 
Kitaigorodskii et al. (1975). The TMA spectrum is an extension of the self-similar 
JONSWAP (Joint North Sea Wave Project) spectral shape (Hasselmann et al.. 1976) to 
include the effects of finite water depth embodied in the Kitaigorodskii et al. (1975) 
finite-depth transformation factor, O a -(o>,/i) : 



StMaWO = £pm(/X,*) *Af,f m , Vj, h) 



(26) 



where 




( 27 ) 



23 



(28) 
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= Phillips constant, 

= 2 ref, 

= frequency, 

= peak spectral frequency, 

= peak enhancement factor, and 

= o a = peak spectral width on low frequency side of f m . 
= o b = peak spectral width on high frequency side of f m . 



(29) 



£ P \j iff m v) is the Pierson-Moskowitz (1964) spectral shape for fully-developed deepwater 
waves, and o a a b ) is the JOXSWAP spectral shape function. 

Thompson and Vincent (1983) developed an approximation for <l> K (<x>,h) using linear 
theory and redefining the two independent variables (co . h) as a single variable (co /; ): 



1 2 

c~2 ( ° h 

O a -(co, h) ^ <l>(co h ) = . 

U-4r (2 -co,) 2 



CO/, < 1 
0) h > 1 



(30) 



where 



co, = 2n (31) 

However, the approximation (equation 30) is incorrect for a) h > 2 and therefore is mod- 
ified here as: 
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<b K (a)Ji) = 



-> a, /i 



- — r> - 



0J h ) 1 < (o h < 2 



(o h > 2 



(32) 



0 ( 60 /,) approaches a value of one in deep water and a value of zero as the depth de- 
creases (Figure 4). 




Figure 4. Kitaigorodskji et al. (1975) finite-depth spectral transformation function, 
<D a (w, h) (from Bouws et al., 19S5), 

Bouws et al. (19S5) developed the following empirical parameter expressions for in- 
put into the TMA data set: 



a = 0.0078 K 049 


(33) 


y j = 2.47 K 0 ' 39 


(34) 
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o a = 0.07 f<f m 

°b = °- 09 f>fm 



(36) 



where 1V l0 = wind speed at z = 10 m and 

k m = wave number for waves at peak frequency. 

These expressions are valid for all water depths in the TMA data set. Thus, the TMA 
spectrum requires specification of iV l0 , h. andj^,. 

B. NUMERICAL APPROACH 

To compare wind and wave forcing for various conditions using finite-depth spectra, 
it is assumed that: 

1) Wave conditions are solely the result of a local wind (i.e. there is no swell 
present), and therefore the wind direction (0) and mean incident wave direction 
are identical, and 

2) Sea conditions are fully arisen (i.e., the sea is not duration-limited). 
Comparisons are made for inputs of gentle and steep planar beach slopes (tan /? = .01 
and .10), wind speeds (10-30 ms -1 ), all wind directions, and peak spectral frequencies 
commensurate with wind speed. 

Surface wind stress is determined from equation (11). Q is determined from meas- 
urements taken by Geernaert et al. (19S7b) in open waters of the North Sea for a depth 
of 10 m (Figure 2). These Q values are increased by 33% to account for surf zone 
roughness. The 33% increase will be addressed later in the Chapter VI. Peak frequency 
is determined for the Pierson-Moskowitz spectrum (Carter, 19S2): 



Peak frequency is assumed not to shift during shoaling (Hughes, 1984). 

After input of tan /? and 1F 10 . peak wavenumber (k m ) is calculated using the linear 
theory relationship starting at an initial water depth of 20 m. The T.MA spectrum is 
determined from equation (26). The RMS wave height is obtained by integrating the 
TMA spectrum: 




(37) 




0.5 



(38) 
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where only sea waves (maximium frequency of 0.5 Hz) are considered. 

Thornton and Guza (1982, 1983) found an envelope curve relating 77 rms to depth in 
the inner surf zone: 



#rms ~ 0.42 h (39) 

The 77 rms determined by equation (38), which includes the effects of variable depth, is 
checked against equation (39) to ascertain if wave breaking due to shoaling had com- 
menced. If the ratio H rm Jh is less than required for surf zone wave breaking, water 
depth is reduced and a new 77 rms determined. This process for calculating breaker depth 
iterates until equation (39) is satisfied. Surf zone width (w ft ) is then determined from: 



h 

tan P 



(40) 



The radiation stress gradient (6S yx l ox) and surface wind stress (t(f) are both inte- 
grated with respect to x b to allow comparison of a wave force ( S yx ) and a wind force 
(t^ x b ), units of force per unit longshore length. The integration is confined to x b 
since the radiation stress gradient is conserved outside the surf zone where there is no 
wave breaking. Momentum exchange between the wind-driven longshore current out- 
side x b and the wind- and wave-driven longshore current inside x b is assumed negligible. 
Using equation (2). wave force is expressed as: 



C ^6 C Sy ^ | ^ 

5 y* = Jo 17" dX = T ff ™Axn(f) PS # COS 0 ( 41 ) 

where the energy is expressed in terms of variance, which is: 

0.5 

^TMAxnff) = 7 Xma(/) n if) A f ( 4 ^) 

0.0 

and n(f) is expressed in a frequency-dependent form. The factor of 1/2 is an adjustment 
factor for broad-banded wind-driven waves which incorporates the effect of short- 
crestedness on the ratio of radiation stress to the average energy content of the waves 
(Battjes, 1972b). 

Wind force is constant over the surf zone and is determined from equation (11): 
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(43) 



T y X b = J o Ty dx = Pa C d \ IV \ W sin 6 x b 

C. RESULTS AND DISCUSSION 

Resultant spectra, wind force, and wave force for both a gentle beach slope and a 
steep beach slope are illustrated as cases A through F in Figure 5 and Figure 6. H rms , 
h , f m , wind speed, Q. and x b for each of the six cases are listed in Table 2. Increasing 
wind speed effect for a gentle beach slope (tan /? = 0.01) is depicted in cases A, B, and 
C, and for a steep beach slope (tan /? = 0.10) in cases D, E, and F. Wind and wave 
direction are relative to normal beach incidence. The three spectra previously described 
(Pierson-Moskowitz, JOXSWAP, and TMA) are included to illustrate the effect of 
shoaling on the wind-wave spectra, as well as the expected shift to lower peak spectral 
frequency as wind speed increases. The H rms and depth (h) at commencement of 
breaking also increase with increasing wind speed. 

Table 2. CASE PARAMETERS FOR WIND VERSUS WAVE FORCE COM- 



PARISONS (see Figure 5 and Figure 6). 



Case 


tan /J 


»ind 

speed 

(ms -1 ) 


Hrms 

(ill) 


h (m) 


x b (m) 


fm (Hz) 


Q 


A 


0.01 


10 


0.7 


1.6 


160 


0.13 


0.0021 


B 


0.01 


20 


3.0 


7.0 


700 


0.06 


0.0033 


C 


0.01 


30 


7.6 


18.0 


1800 


0.04 


0.0053 


D 


0.10 


10 


0.7 


1.6 


16 


0.13 


0.0021 


E 


0.10 


20 


3.0 


7.0 


70 


0.06 


0.0033 


F 


0.10 


30 


7.6 


1S.0 


ISO 


0.04 


0.0053 



The wind force is directly proportional to the width of the surf zone (x b ), which is a 
function of tan /?. As a beach steepens, surf zone width narrows and thus there is a re- 
duced area for surface wind stress to act. Wave force is independent of tan /?. 

Wind direction dependency is inherent in equations (41) and (43) where: 

S : J oc sin d cos 8 = sin(20) (44) 
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Figure 5. Comparison of Mind and wave force for a gentle beach slope (tan /i = 
0.01) Mith associated wave spectra. 
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Figure 6. Comparison of wind and wave force for a steep beach slope (tan /J = 0.10) 
with associated wave spectra. 
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(45) 



t y oc sin 0 

Thus the maximum force occurs at ± 90° (alongshore) for wind and ± 45° for waves. 
Wind force increases more rapidly than wave force for increasing 6. When winds are 
blowing offshore (i.e. absolute wind directions > 90° ), waves do not form because of 
fetch limitations and thus only wind force exists, creating a weak, but totally wind-driven 
longshore current. 

The dependency of wind and wave force on wind speed is not as readily apparent 
because of the interrelationship of wind speed with the various TMA parameterizations 
(equations 33 to 37). The absolute force of both wind and waves increases as wind speed 
increases. The increase in wave force is due to the increase in spectral energy. The in- 
crease in wind force is through the wind speed and the increase in surf zone width caused 
by larger waves (generated from the TMA spectrum) breaking farther offshore. Thus 
the wind force's dependency on wind speed is indirectly dependent on the spectral wave 
energy. 

Wave and wind force are compared for increasing wind speed and constant wave 
incident angle (Figure 7). Since wave force is independent of /? . only wind force is il- 
lustrated for both gentle (solid line) and steep (dashed) beach slopes. The amount of 
increase for both wind and wave force over the wind speed range of 5 ms -1 to 30 ms -1 
is nearly four orders of magnitude. Wind force is seen to increase more rapidly than 
wave force for increasing wind speed and constant incident wave angle. Varying the in- 
cident wave angle results in the same conclusion; therefore only 0 equal to 20° is illus- 
trated. 

The more rapid rate of increase for wind force is quantified in terms of a wind force 
to wave force ratio: 



r *l Y 

~y x b 



J yx 



(46) 



The ratio of forces is depicted in Table 3 (tan /? = 0.01) and Table 4 (tan /? = 0.10) for 
the quadrant of wind direction from 0 — 90°. The ratio is symmetric for wind directions 
of 0 to —90° and is a function of slope, wind speed, and wind direction. For steep beach 
slopes, wind force is much less than wave force and thus may be neglected. However, 
for gentle slopes, wind force increases as the wind becomes more alongshore in direction. 
For a given wind direction, an increase of 10 ms -1 in wind speed results in a 50% in- 



31 




Figure 7. Comparison of ind and wave force for increasing wind speed at constant 
wind direction ( 0 = 2U°). 



crease in the ratio. Although this ratio for gentle beach slope is generally less than unity 
for all wind speed and direction combinations, wind force is never less than an order of 
magnitude compared to wave force. Wind to wave force ratio versus beach slope for a 
typical mid-latitude cold frontal passage is depicted in Figure S. Wind and wave frontal 
conditions are defined as a wind speed of 15 ms -1 and 6 = 20°. The wind to wave force 
ratio increases with decreasing beach slope. The ratio is approximately 10% for a beach 
slope of 0.04 and approaches 40-50% for beach slopes less than 0.01. Therefore wind 
force should be considered in studying longshore currents during strong winds on gentle 
beaches and during offshore winds. In Chapter VII, wind to wave force ratios from the 
field experiment will be compared to the results in Table 3 and Table 4. 



Table 3. RATIO OF WIND FORCE TO WAVE FORCE FOR A GENTLE 
BEACH SLOPE (tan [t = 0.01) 



Wind 

Speed 

(ms -1 ) 


Wind Direction (relative to normal incidence) 


0° 


10° 


20° 


30° 


40° 


50° 


60° 


70° 


80° 


90° 


10 


0.00 


0.15 


0.15 


0.17 


0.19 


0.22 


0.29 


0.42 


0.83 


INF 


20 


0.00 


0.23 


0.24 


0.26 


0.29 


0.35 


0.44 


0.65 


1.2S 


INF 


30 


0.00 


0.34 


0.35 


0.3S 


0.43 


0.52 


0.66 


0.97 


1.91 


INF 



• INF - ratio is infinite since there is no wave force 

• waves are assumed to be locally generated and to approach the shoreline from 
the wind direction 
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Table 4. RATIO OF WIND FORCE TO WAVE FORCE FOR A STEEP BEACH 
SLOPE (tan (i = 0.10) 



Wind 

Speed 

(ms- 1 ) 


Wind Direction (relative to normal incidence) 


0° 


10° 


20° 


30° 


40° 


50° 


60° 


70° 


80° 


90° 


10 


0.00 


0.0 i 


0.02 


0.02 


0.02 


0.02 


0.03 


0.04 


o.os 


INF 


20 


0.00 


0.02 


0.02 


0.03 


0.03 


0.03 


0.04 


0.07 


0.13 


INF 


30 


0.00 


0.03 


0.04 


0.04 


0.04 


0.05 


0.07 


0.10 


0.19 


INF 



• INF - ratio is infinite since there is no wave force 

• waves are assumed to be locally generated and to approach the shoreline from 
the wind direction 
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0.50 



Wind Speed 15 ms 







Beach Slope 



Figure 8. Wind to nave force ratio versus beach slope for a typical mid-latitude cold 
frontal passage. The wind and wave incident angle ( 9 ) is relative to 
beach normal. 
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IV. FIELD EXPERIMENT 



A. INTRODUCTION 

To develop an improved understanding of coastal processes (currents, waves, 
sediment transport, and nearshore geomorphology) under a wide variety of conditions, 
a one month-long field experiment was held in the autumn of 1986 at the U.S. Army 
Corps of Engineers, Waterways Experiment Station's Coastal Engineering Research 
Center (CERC), Field Research Facility (FRF) at Duck, N.C. Duck is on the Atlantic 
Ocean side of the North Carolina Outer Banks and is located approximately 10 miles 
north of Kitty Hawk, N.C. This major field experiment was entitled SUPERDUCK and 
included nearly 20 scientific investigations. 

B. SITE DESCRIPTION 

The FRF site (Figure 9) is located along a 100-km unbroken stretch of shoreline 
of an extensive barrier island formation known as North Carolina's "Outer Banks." The 
barrier island upon which the FRF site is located is bordered by the Atlantic Ocean to 
the east and Currituck Sound to the west. Except for a few fishing piers, there are no 
major coastal structures or littoral barriers along the entire reach of shoreline. The site 
has a tidal range of 0.5 to 2.0 m and regular offshore bathymetry free of features which 
may affect wave climate. The beach is characterized by a mean foreshore slope of 1:10, 
a single longshore linear bar system which can become three-dimensional, and a mean 
slope of 1:100 offshore of the bar. Shoreline excursion due to tidal fluctuations is typi- 
cally 10 m. Mean sediment grain size is approximately 0.5 mm. 

FRF personnel support includes a staff of nine (scientists, technicians, divers, and 
heavy equipment and computer operators). Extensive equipment support includes a 
LARC-V amphibious vehicle, a Coastal Research Amphibious Buggy (CRAB), auto- 
mated data processing (ADP) equipment, a heavy-duty forklift, four-wheel drive and 
all-terrain vehicles (ATV's), laboratories, electronic workshops, a maintenance garage, 
and high-resolution Zeiss Elta-2 electronic surveying systems. 

Several planning workshops were held months prior to the experiment to synergize 
the efforts of all the principal investigators. The FRF staff and SUPERDUCK partic- 
ipants jointly installed more than 70 electronic sensors in the nearshore environment for 
the experiment. The overall set-up of the various instruments at the FRF site is illus- 
trated in Figure 10. The longshore current meter array was approximately 500 m in 
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Figure 9. SUPERDUCK field experiment site. 
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length and 25 m ofTshore. Its primary purpose was for the investigation of infragravity 
waves, however its data will also be used in this research to investigate temporal and 
spatial variability of longshore momentum flux and to provide an extra data station for 
sled transects of the surf zone. Unfortunately, the cross-shore array of current meters 
and pressure sensors was not functional for the experiment. The third cross-shore sen- 
sors were sonar bottom profilers. A high-resolution linear pressure sensor array in 8 m 
depth provided offshore directional wave data. High resolution bathymetry data were 
collected on a nearly daily basis between profile lines 165 and 275 by laser positioning 
the CRAB with the Zeiss surveying system as it methodically traversed the area along 
the profile lines. Example products of the surveying effort are included in Figure 1 1 and 
Figure 12. A major difference between this and previous large-scale nearshore exper- 
iments was the extensive bathymetric data acquired using the CRAB. 

C. EQUIPMENT DESIGN AND OPERATION 

Only a brief description of the sled and its instrumentation is provided here. 
Martens and Thornton (1987) describe in detail the design, instrument calibration, and 
specification of the sled system. The all-aluminum sled (Figure 13) was constructed by 
the CERC Field Research Facility stall in Duck, N.C., based on a Naval Postgraduate 
School design. The sled foundation is two runners 3 m in length on which was placed 
a 3 m x 3 m aluminum sled frame with a 9 m aluminum sailboat mast. The frame was 
approximately one quarter of a meter above the sled runners. All the instrumentation 
was then attached to this frame. Significant additional weight was placed on the runners 
to provide stability and prevent dragging in the robust surf zone. The sled was designed 
with a low profile to minimize flow disturbance. 

The differential pressure slope array concept was patterned after Bodge and Dean 
(19S4) with some major modifications. Two-inch O.D. stainless steel threaded pipe and 
pipe fittings were assembled in the shape of a capital "H" to form two adjacent triangular 
slope arrays with a common vertex at the center (see Figure 13). The five differential 
pressure ports were located at the four ends of the "I" and one at the center. The pipe 
cavity was filled with a premium grade vacuum pump oil, similar to mineral oil, which 
had a low viscosity coefficient and was assumed incompressible. Four Sensym model 
LX06002D monolithic differential pressure transducers were installed at the ends of the 
"1", just inboard of the diaphragms which were exposed to the sea. These differential 
pressure transducers provide an output voltage proportional to the applied differential 
pressure with a response time of 0.1 seconds. The pressure at the center diaphragm was 
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Figure 10. SUPERDUCK field experiment instrument locations. 
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Figure 11. Three-dimensional bathymetry example, 16 October 1986. Short ver- 
tical lines represent locations of cross-shore and longshore instruments. 
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Figure 12. Two-dimensional bathymetry example, 16 October 1986. 
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Figure 13. Side and top view of the SUPERDUCK sled 



hydraulically transmitted to the positive port of each differential pressure sensor, while 
the sea pressure at that location was sensed by the negative port of the sensor. A 
Paroscientific model 8020 digiquartz pressure sensor was installed adjacent to the center 
diaphragm to measure the absolute pressure. 

Current velocities were measured using three Marsh-McBirney bi-directional 
electromagnetic current meters with spherical (4-cm diameter) probes mounted on the 
sled. These current meters are ideally suited for surf zone studies by providing accurate 
measurements in free stream turbulence as well as combined steady, oscillatory currents 
(Guza et al., 1986a). Both variable vertical stacking and slope array comparisons were 
considered in positioning the current meters on the sled frame. Two current meters were 
mounted at the centroids of the two triangular-shaped slope arrays at elevations 0.7 and 

l. 0 m above the bed, with the third current meter midway between at an elevation of 1.5 

m. During sled deployment, the current meters were always positioned on the "up- 
current" side of the sled so as to avoid flow contamination by the sled mast. 

The current meters' directional electrodes were aligned with the longitudinal axis of 
the sled frame by eye and are considered to be within 1 — 2° of the true axis. During 
data analysis, a computational rotation is done for true alignment based on conservation 
of momentum flux outside the surf zone. 

A single-point R.M. Young model 05103 wind monitor was mounted atop the mast 
to provide "in-situ" measurement of vector wind in the surf zone with simultaneous wave 
momentum measurements. 

All equipment were battery-powered by 12 sealed lead-acid batteries housed in a 
watertight canister mounted on the sled. The analog data were PCM (pulse code mod- 
ulated) encoded and then telemetered ashore by a system described by Lowe ct al. 
(1972). The data were then decoded and recorded on 9-track magnetic tape. Data ac- 
quisition times for all the various SUPERDUCK experiments were referenced to a pre- 
cise Eastern Standard Time (EST). 

Visual observations of waves breaking on the sled mast were marked by an electrical 
trigger which registered a pulse on the data tapes. Although extremely tedious, this ef- 
fort allows determination of the occurrence of wave dissipation events. 

Precise orientation of the sled, with respect to true north and the SUPERDUCK 
experiment area, was obtained with the use of the FRF's Zeiss Elta-2 electronic survey- 
ing system. Two reflective prisms (5-cm diameter) were mounted at each end of a 2.4 
m longitudinal mast spreader located approximately 7 m above the sled. From triangu- 
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lation. the relative positioning of these two prisms and their absolute positioning within 
the FRF coordinate system provided the required precise orientation (within 0.5°) of the 
current meters and slope array. Relative three-dimensional positions between various 
points on the sled and the prisms were obtained with separate Zeiss shootings prior to 
and after the sled's deployment in the SUPERDUCK experiment. 

Prior to selecting an area of the surf zone to transect, the previous day's bathymetric 
profiles and contours were studied and changes in morphology noted. The geographic 
area with the highest degree of straight-and-parallel isobaths was selected for each day's 
operations. The runs were made during daylight and scheduled for 1-2 hours before and 
after a tidal change, so as to minimize the tidal effect and variation in mean shoreline 
location. This timing also coincided with the recording of all other SUPERDUCK data, 
which was accomplished within + 2 hours of a tidal change on a 24-hour basis. The sled 
was initially towed offshore by the FRF's CRAB to a position just beyond the breaker 
zone for the first run of a surf zone transect. Subsequent runs were made at the point 
of maximum breaking, on top the nearshore bar, and in the nearshore trough. Each run 
was at least 35-minute in duration. The sled was brought inshore for each subsequent 
run of a transect by towing the sled with a chain tether attached to a four-wheel drive 
forklift located on the beach. Zeiss Elta-2 laser surveying shots of the sled's prisms were 
taken at the beginning and end of each run, with many runs having additional surveying 
shots during the run itself. 

In addition to the sled's anemometer, simultaneous meteorological data of air tem- 
perature, sea surface temperature, atmospheric pressure, and wind velocity were re- 
corded by other SUPERDUCK investigators at the end of the pier and near the FRF 
laboratory. 
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V. PRELIMINARY DATA ANALYSIS 



Data analysis of the SUPERDUCK field experiment is addressed in the next four 
chapters. Preliminary data analysis is addressed in this chapter and includes data 
screening for data gaps and signal interference, instrument orientation, three- 
dimensional depiction of the mean current structure over the bar, and a discussion of 
homogeneity and stationarity for the experiment. Since surface wind stress is a major 
concern and few wind stress studies have been conducted in the surf zone, Chapter VI 
is entirely devoted to the parameterization of surf zone surface wind stress. The 
longshore momentum balance is addressed in Chapter VIII. The relative importance 
of each momentum term is calculated. Spatiallv-variable cj s are determined as residuals 
from the momentum balance. A longshore current model for planar beaches is modified 
to account for barred topography and its results are compared with observations in 
Chapter VIII. 

A. INITIAL DATA SCREENING 

The analog data were PCM decoded and then digitized at S Hz. Time series plots 
for all 16 channels of information were produced for the entire data collection period. 
These records were visually perused for data gaps and signal interference. Basic statistics 
(mean, variance, skewness and kurtosis) were calculated and formated in tabular form. 
Temporal plots of two-minute means and variances were produced. This information 
was used to eliminate bad data collection periods and data channels. Data collection 
achieved a 9S.5% successful collection rate with over 40 hours of data collected over a 
10-day period. Nearly all of the eliminated data were identified at the beginning or end 
of a data run and therefore temporal continuity was maintained. 

The wind direction portion of the R.M. Young wind monitor was not calibrated 
before its unexpected destruction by wave forcing on IS October. Therefore the mean 
wind direction measured at the sled is unavailable. The mean wind direction measured 
at the end of the FRF pier (SethuRaman et al., 19S7) was used for this research. Mean 
wind speed after IS October was also taken from the end-of-pier data. 

Occasionally, one or two current meters would be partially or totally above water 
during a sled data collection period, such as when the sled was atop the bar at low tide. 
These periods were identified from visual inspection of the time series, supplemented by 
field logbook entries. These data were not utilized in the radiation stress calculations. 
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B. SLED INTER-ORIENTATION. 



A portable Zeiss laser ranging system was used to precisely measure the sled 
instrumentation positions relative to the sled. A hand-held reflective prism was placed 
at a designated measurement point on the sled (e.g.. the right rear corner of the sled 
frame top) and that position designated as a three-dimensional reference point from 
which all other sled points of interest (prisms, current meter probes, slope array 
diaphragms, etc.) were measured. A complete laser measurement of the sled was con- 
ducted on 26 September 19S6 and 20 October 1986. Sled modifications, such as instru- 
ment adjustment, mast replacement, etc., were noted in the field logbook so as to enable 
selection of the proper survey day for subsequent calculations. The Zeiss measurements 
were supplemented with engineer tape measurements. 

The sled instruments were aligned to a three-dimensional local cartesian coordinate 
system referenced to the top of the sled frame by a rectangular coordinate rotation: 

*new = -'old cos 0 + y old sin 0 

J new -'old 8 "b.Vo]d COS 6 

where 0 is the angle of rotation. A table of absolute three-dimensional distances be- 
tween all points of interest was then generated. These distances were repeatedly used 
throughout the analysis. 

During field operations, three-dimensional instrument coordinates, with respect to 
the FRF coordinate system were obtained from the sled prism three-dimensional coor- 
dinates by: 

• relating the mast spreader, slope array, and current meter system longitudinal and 
vertical axes to the sled frame longitudinal and vertical axes, 

• determining instrument pitch and yaw angles ever} - 64 seconds from linear in- 
terpolation of the two-to-flve survey fixes for a sled run collection period. All data 
were then adjusted every 64 seconds for the current sled pitch and yaw angles, 

• adjusting current meter and anemometer heights above the bottom even 64 sec- 
onds for sled pitch and yaw angles, and 

• correcting the distance between current meters for the sled's pitch and yaw angles. 

Since all data were taken from a mobile sled rather than from fixed-bed instrumen- 
tation, it was imperative that the data were accurately translated into a coordinate sys- 
tem relative to the local bathymetric contours. Use of the Zeiss laser ranging system to 
accurately measure the sled points of interest, to fix the sled in the experiment area, and 
to provide comprehensive bathymetric data coverage, enabled this translation. 
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C. BEACH ANGLE ORIENTATION. 

Location coordinates obtained from Zeiss-sighting of the sled's prisms were con- 
verted into sled pitch and yaw angles so as to correct the sensor data with respect to the 
horizontal and vertical and the FRF coordinate system. Pitch and yaw angles were 
calculated every 64 seconds from linear interpolation between the two most current Zeiss 
sightings to account for any sled movement during the data run. The instrument array 
orientation was then numerically rotated relative to the measured bathymethv. The se- 
lection of 64 seconds coincided with the ensemble averaging time for the spectral calcu- 
lations. 

Sled data were acquired from 11-18 October and 20-21 October 1986. Locations of 
sled data runs within the SUPERDL’CK experiment site are illustrated in Figure 14. 
Incomplete survey sightings of the sled prisms due to salt incrustation on the prism lens 
and Zeiss operator error prevented determination of precise sled location on 13, 14 and 
20 October. Therefore these three days of data were not used. 

A complete bathymetric survey was not done on 11, 17 and 21 October due to high 
surf conditions or equipment availability. For sled collection days when no coincident 
bathymetry was surveyed, the previous or next day's bathymetry was appropriately se- 
lected based on minimizing any storm-induced morphological change (see Table 5). 

A mean bathmetric contour direction was determined by applying a least-squares 
regression fit to the survey data offshore and located within approximately 100 - 150 m 
of the nearshore bar. This is the area of maximum bathymetric gradient with minimal 
slope farther offshore. The bathymetry inside the bar was more complex than offshore 
the bar, however the shoreward distance over which it could affect wave action was small 
0(30-50 m). The most influential bathymetry was assumed to be that located on and 
offshore the bar. 

The longshore length scale over which this regression fit was applied was weighted 
(2 3 - 1/3), with emphasis on the portion of the beach which favored the approaching 
waves (e.g., for waves approaching from the northeast, the northern offshore contours 
were favored by a factor of 2). The total distance assessed was varied from 100 to 250 
m (or to maximum extent of precision bathymetry) with resultant mean orientation 
varying less than 0.5° between various distances considered. The mean orientation of 
the four "best-fit" straight lines (corresponding to 0.5 m contour intervals) was chosen 
as the mean bathymetric contour. The beach orientation angle was the angle between 
this mean contour and the coordinate system's y-axis. The beach orientation angle was 
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Figure 14 . Sled data run locations. Open circles are sled run locations and closed 
circles are fixed-bed current meter locations. 
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then assigned as a constant for that transect and used in conjunction with instrument 
orientation which is described later. 



Table 5. BATHYMETRY SELECTION AND BEACH ORIENTATION AN- 
GLES 



Sled data dav 


Bathvmetrv dav Profile 1 


Beach orientation ancle (°) 2 


11 


12 


250 


0.5 


12 


12 


250 


0.5 


15 


15 


255 


4.0 


16 


16 


197 


-0.5 


17 


16 


195 


-0.5 


IS 


IS 


195 


1.5 


21 


22 


170 


0.8 



1 Profile locations are shown in Figure 10. 

2 Beach orientation angle is measured between the FRF coordinate system's 

y-axis and the mean bathymetric contour from on top and 
olfshore of the bar. The angle was measured over a longshore length 
scale varying between 100-250 m. 



A combined angle ( 0 C ) was used to convert from the yaw angle measured by the 
Zeiss to the correct rotation angle upon which to coordinate transform the u,v velocities. 
This combined angle was the sum of the sled yaw angle, the beach orientation angle, the 
angle between the current meter longitudinal axis and the sled centerline, and the angle 
between the mast spreader longitudinal axis and the sled centerline. This combined an- 
gle (0 C ) was used in equation (47), restated here with velocities substituted for coordinate 
locations, to perform the coordinate transformation: 

“new “old 4" v old ^itl B c 

(hX) 

V'new = "“old °c + Void cos °c 

Fixed-bed current meter velocities were also converted to the FRF coordinate sys- 
tem using equation (48) where 6 C is field-measured orientation angles with respect to the 
FRF coordinate system. 

D. CURRENT METER ORIENTATION PRIOR TO A TRANSECT. 

Since the radiation stress gradient is conserved outside the surf zone, any difference 
in measured S yx (J) between the current meters on the sled is considered to be due to er- 
rors in current meter alignment with respect to each other and to the sled, current meter 
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calibration, or to imprecise orthogonal orientation of the dynamic measurements with 
respect to the assumed straight-and-parallel bathymetric contours. The contours outside 
the bar were nearly straight-and-parallel (see Figure 12 and Appendix B). Therefore, 
current meter alignment or gain was considered in error. Orientation of the current 
meters' axes with the sled's axes was considered to be within 1 — 2°. 

To correct for the radiation stress difference between two current meters outside the 
surf zone, one current meter was assumed to be true in both gain and alignment. The 
other current meter was then rotated numerically using equation (48), or adjusted in 
gain, until the radiation stress gradient was minimized. 

Two current meters and two adjustment criteria (gain and rotation) represent four 
initial parameters for adjustment: 

• assume the front current meter is correct and adjust the rear current meter in gain 
(RGA1N) or rotation (RTWIST), or 

• assume the rear current meter is correct and adjust the front current meter in gain 
(FGA1N) or rotation (FTW1ST). 

Because the sled instruments undergo significant stress in the surf zone, all sled transects 
commenced with an orientation run outside the surf zone to allow for daily reorien- 
tation. A numerical orientation prior to each transect's calculations was employed dur- 
ing data analysis. The resulting possible orientation adjustments for each transect and 
parameter (RGA IN, RTWIST, etc.) are shown in Table 6. 

The adjustments and resulting total radiation stress (S y J) were found to be linearly 
related, therefore an iterative method of linear interpolation, or "regula falsi" method 
(Gerald, 197S) was utilized to determine the orientation required to minimize S yx be- 
tween the current meters. Convergence to within a radiation stress gradient tolerance 
of 10 -3 Jnf 2 was generally obtained within six iterations. S yx (J) spectra with these ad- 
justments were perused to ensure that the strict numerical definition of a conserved ra- 
diation stress gradient achieved a realistic and comparative spectral shape for the current 
meters. 

The selection of orientation by gain or rotation does have constraints when dealing 
with bi-directional radiation stress spectra, which commonly occurred during SUPER- 
DUCK. Examples of typical bi-directional spectra are schematically illustrated in 
Figure 15. The spectra indicate swell with a 10-second mean period from the southeast, 
while wind waves with a 5-second period and significantly more energy approached from 
the northeast. Spectra from the front current meter are assumed correct and are de- 
picted by a solid line, whereas the spectra from an incorrectly oriented or calibrated rear 
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current meter (and requiring gain or rotation to make its spectra coincide with the front 
current meter spectra) are depicted by a dashed line. In the first example, gain is the 
obvious correction whereas in the second example, rotation is the obvious correction. 

All three current meters were calibrated in a flow tank before and after the exper- 
iment by different laboratories. The agreement in pre- and post-calibration was + 3% 
in gain and 0.01 — 0.04 ms -1 in offset. This agreement is considered excellent for meters 
exposed to the robust surf zone. With this agreement, coupled with the subjective ana- 
lyses of the different spectral shapes corresponding to different orientation combina- 
tions, and the field accuracy of current meter alignment (+ 2°), it was decided to utilize 
current meter rotation as the orientation criteria. Rear current meter rotation 
(RTWIST) was arbitrarily selected over front current meter rotation for the entire data 
analysis, except for Run 18-1 for which front gain was used. Resultant RTWIST angles 
are within the values of field accuracy alignment. Sensitivity to this rotation and the 
effect of other rotations (FTWIST only or both FTWIST and RTWIST) is addressed in 
the error analysis discussion in Chapter VII. 

After completing this rotation for the rear current meter, the side current meter was 
rotated (STWIST) in a similar manner and its S v [ values served as a first-order com- 
parison to the other two current meters values for S y [. These orientation rotations were 
performed after the 6 C coordinate transformation were fixed for the calculations inside 
the surf zone for that particular transect. This post-calibration procedure provided a 
means of improving the accuracy of the current meter calibration and alignment. The 
procedure is essential to obtaining realistic physical results from the data and is a dis- 
tinguishing feature of this research. 

E. MEAN CURRENT OBSERVATIONS 
1. Overview 

Three-dimensional portrayal of the mean current structure is illustrated in Fig- 
ure 16 and Appendix A (Figure 35 through Figure 40). All current vectors represent 
34.1 -minute means. The numbers indicate the data run location. The short dotted lines 
above each run location indicate the mean water level (MWL) measured by the sled's 
pressure sensor for that location and time. The different arrowheads delineate the height 
of measurement as illustrated in the three-dimensional box. MSL is indicated by the 
long dashed line. Only current meters which were completely submerged for the entire 
run duration are depicted. Thus some run locations have only one or two current vec- 
tors. 
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Table 6. INSTRUMENT ORIENTATION PARAMETERS. Rear and front cur- 
rent meter gains and rotations ("TWIST") required to conserve the radi- 
ation stress gradient outside the surf zone. 



Run 


Rear Gain(%) 


Front Gain (%) 


Rear Twist (dee) 


Front Twist (dee) 


12-1 


-5.25 


5.54 


-0.S1 


O.SO 


15-1 


1.04 


-1.03 


-0.29 


0.27 


16-2 


4.57 


4.3S 


-2.86 


2.34 


17-1 


-0.5S 


0.5S 


0.12 


-0.11 


1S-1 


5.20 


-4.94 


-2.14 


1.89 


21-1 


1.40 


-1.44 


0.02 


-0.02 


21-6 


- 


- 


-0.41 


0.40 



• indicates no iterative convergence or acceptable spectrum shape. 

• 'Twist's', or current meter rotations, are in units of degrees. 



The sled was first deployed on 1 1 October during the first of two SUPERDL CK 
storms. High surf (H mo greater than 3 m) prevented sled deployment beyond the surf 
zone, therefore there was no instrument orientation run for that day. The high surf also 
prevented completion of a ntinigrid survey by the CRAB. Therefore the bathymetric 
profile, as well as the instrument orientation parameter was obtained from the 12 Octo- 
ber data. The unrealistic c.j values later calculated for 1 1 October highlight the necessity 
for an instrument orientation prior to a transect. The 1 1 October data set was therefore 
deleted from further consideration. 

The 12 October data revealed a meandering rip current. Two transects were 
completed with Runs 12-1 and 12-4 being the orientation runs outside the surf zone. 
Run 12-4's data were not used because of a Zeiss system malfunction. Therefore Run 
12-1 is the only data run outside the surf zone in Figure 36. The strong offshore flow, 
0(0.5 ms -1 ), at all three levels above the bed was substantiated by comparison with 
current vectors from the longshore array. Although this rip current eliminates homo- 
geneity for S } J calculations, it provides a unique data set for a rip current study separate 
from this research. 

October 15 is the first day with an established bar-trough system. Run 15-1 is 
the orientation run outside the surf zone. Strong flow, 0(1.0 ms -1 ), is apparent inside 
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Figure 16. Mean current structure across the bar on 16 October 1986. Numbers 
represent data runs, short dotted lines are MWL's, and distance 
alongshore is arbitrary. 
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the surf zone with a 50% reduction in speed outside the surf zone. All speeds are nearly 
uniform in depth with slight speed reductions for sensor elevations closer to the bed. 
Multi-level current directions in the trough and outside the surf zone are co-linear. 
Angular deviation of currents with height (maximum 20°) are apparent at the location 
of maximum breaking, with the lower sensor elevations exhibiting a bias toward offshore 
flow near the bed. 

October 16's data have the strongest velocities of all the sled data with speeds 
exceeding 1.0 ms -1 . Run 16-2 was the orientation run outside the surf zone. Current 
directions and speeds were nearly uniform with depth, with current speeds exhibiting a 
slight decrease closer to the bed. Strongest flows were on the offshore side of the bar. 

Significantly weaker currents were measured on 17 October. Run 17-1 was the 
orientation run outside the surf zone. The three trough runs all produced negligible 
current. A combination of low tide and low wave heights [// rms 0(O.2ms -, )3 on this day 
may have caused most waves to break on the bar and thus cause the negligible current 
in the trough. The onshore Hows at Run 17-4 and 17-5 were due to the current meters' 
location in the upper water column of the shallow water. 

Strong flows, 0(1.0 ms -1 ), were again measured on IS October. Runs 18-1 and 
1S-6 were orientation runs outside the surf zone for two separate transects. Run 18-6 
was the start of a second transect which was abruptly terminated when the sled mast 
snapped in the high surf during Run 1S-7. A nearly uniform horizontal velocity gradient 
(dV/cx) is observed with a maximum current located offshore and just before the bar 
crest. Vector uniformity is observed. 

Two transects were conducted on 21 October. Runs 21-1 and 21-6 were orien- 
tation runs outside the surf zone. Unfortunately current velocities during both transects 
were too weak and variable in direction for meaningful analysis. Hence this data set was 
eliminated from further consideration, leaving only four (15- IS October) of the original 
nine transects for further analysis. 

2. Discussion 

Profiles of longshore currents across the barred beach generally showed signif- 
icant spatial variability. Maximum flow is observed at the steepest offshore slope of the 
bar, generally just below the bar crest, while weaker currents are found outside the 
breaker zone and in the trough. A typical distribution of longshore currents across a 
barred beach is illustrated in Figure 16. Longshore currents are generated after wave 
breaking (Run 3 and 4) and are extended outside the surf zone, possibly by turbulent 
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momentum exchange. The vertical distribution of longshore and cross-shore velocities 
are compared to normalized water depth (Figure 17). Current measurements for one 
run are denoted by a continuous line connecting only those current meters which were 
entirely submerged for the entire run. Thus some shallow locations have only one or two 
current means illustrated. The plot of October 16 longshore velocities serves as a legend 
and has numbers beside the mean current symbols. The numbers represent the chrono- 
logical order of runs in a transect, with the most offshore run being labeled as number 
1. 

The magnitude of the cross-shore flow is generally much less than the longshore 
flow and is of 0(0.2 ms -1 ). The longshore flow velocities on three of the four data days 
are of 0(1.0 ms -1 ). The cross-shore velocities exhibit a relative onshore flow in the upper 
third of the water column and relative offshore flow in the lower third of the water col- 
umn. The lower mean currents after wave breaking tend to flow more offshore, thus 
exhibiting a limited "spiral" distribution with height and indicating the presence of an 
undertow flow (an offshore flow near the bed). The observations are in agreement with 
cross-shore theory (e.g., Svendsen et al., 1987). The longshore velocities indicate nearly 
depth uniform flow with a slight velocity increase with elevation. This depth uniform 
flow confirms the appropriateness of the depth-integration assumption for longshore 
flow. Longshore vertical current structure did not exhibit a logarithmic profile, therefore 
the profile method for the determination of bed shear stress coefficients is not applicable. 

3. Stationarity and Homogeneity 

Stationarity over the transect was investigated by examination of S y [ calculated 
from measurements outside the surf zone by a fixed-bed pressure and velocity ('puv') 
sensor designated "south tripod". The formulation for this calculation is discussed later 
in equations (7S) and (79). Sequential 34.1 -minute spectra and resulting S y ][ were cal- 
culated for 15- IS October (Figure IS and Appendix A (Figure 41 through Tigure 43). 
The short lines between the arrowheads at the top of the S y [ graphs represent sled run 
durations, with the numbers above the lines identifying the sled run number of that day. 
These comparisons are considered a strict indication of stationarity because of the sen- 
sitivity of to changes both in wave height and wave incidence angle. 

The 'puv' data were collected by the FRF for four-hour periods, therefore these 
data were not completely synoptic with sled deployment times. The 'puv' data were 
digitized at 2 Hz. Spectral analysis was accomplished by Bartlett ensemble averaging 
which resulted in 31 degrees of freedom. Chi-squared confidence limits for 95% confi- 
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Figure 17. Mean cross-shore U and longshore velocity V profiles versus normalized 
water depth for 15-18 October. The 16 Oct V profile serves as a legend 
for the data symbols, where the numbers denote chronological order 
of runs in a transect with the most offshore run labeled as number 1. 
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dence yields a range for spectral estimates (S yx (/)) between 0.66 and 1.S1 of the spectral 
estimate calculated. Ninety-five percent confidence limits for the integrated or total ra- 
diation stress (S v [, see equation SO) were based on the Fisher-Z transformation (Miller 
and Freund, 19S5): 




where the correlation coefficient (r) is: 



y = 



< 7 ,( 7 , 



(51) 



and the covariance is calculated by the integrating the co-spectrum of u,v. 

The October 15 runs inside the surf zone (15-2 through 15-4) appear to have a 
reasonably constant (within the bounds of the confidence limits) S y £ value with no more 
than a 10% variation about the mean. Both October 16 and 17 data sets have a constant 
S y T with minimal variability. October IS data also appears to have a reasonably con- 
stant S } J, with the exception of Run 1S-5 which occurred at the time of an S y £ increase. 
Fortunately, this run is within the trough, where How conditions are not as rapidly af- 
fected by offshore change as are locations on the offshore side of the bar. Even with this 
S..T increase, the variation is no more than 20% about the mean. 

Tidal variation during an individual data run was approximately 5% of the 
mean water depth. Tidal variation over an entire transect ranged from 14% - 50% of 
the mean transect depth, with the longest transect duration having the largest tidal var- 
iation. The tidal range for each transect is illustrated by the short dotted lines in 
Figure 16 and Appendix A (Figure 35 through Figure 40). 

Another test for stationarity is illustrated in Figure 19 where ll mo and wave in- 
cident angle (a) calculated from the 'puv' sensor are plotted from 0500 - 2000 E.S.T. for 
each day and compared with sled run times. // mo values are denoted by an "x" and a 
values are denoted by an "o". Each short line at the top of the graphs represents the 
duration of a sled run, with the number above the line denoting the sled run number of 
that day. Comparison of Figure 19 with Figure 18 and Appendix A (Figure 41 through 
Figure 43) offers some insight into the nature of radiation stress. For example, the de- 
crease in S V J observed on 15 October from 1000 - 1200 was due to a slight reduction in 
wave height and a significant decrease in wave incident angle. In conclusion, no signif- 
icant change in S V J, H mo , or a. was observed during the times of sled deployment. 
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Ocl 16: 1020 - 1054 Oct16: 1054 - 1120 Oct16: 1128 - 1202 OctIG: 1202 - 1236 




i n 



Figure 18. Offshore S yx spectra and SJ from 'puv' gage 621 for 16 October 
19S6. Numbers indicate data runs. Intervals are 95% confidence. 
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Stationarity requirements are therefore assumed satisfied. It should be noted that the 
stationarity being investigated here is measured over relatively long temporal scales 
(measured in hours) whereas the temporal term examined later in Chapter VII is meas- 
ured over much shorter temporal scales (measured in minutes). In other words, the as- 
sumption of stationarity at this stage of the data analysis does not preclude 
contributions from the temporal term (e.g. due to edge waves) to the momentum bal- 
ance. 

Similarity of mean current vectors in the experiment area was considered a test 
for homogeneity inside the surf zone. Similarity was considered satisfied when all cur- 
rent directions were in the same quadrant with speeds ± 50% of an overall mean for the 
same relative on-offshore position. This definition was chosen considering the bounds 
of instrument accuracy, current vertical variability, the number of instruments in the 
experiment area, and the variable elevations of the longshore array current meters. No 
visual evidence of rip currents was observed during 15- IS October. Two-dimensional 
plots of sled-determined current vectors were created for each transect, along with the 
fixed-bed array of mean current vectors, temporally centered for the sled transect. These 
comparisons were then overtaxed on the coincident bathymetry (Appendix B). The 
bathymetry allows for isolation of cases where bathymetric forcing is dominant. Vari- 
able current direction is observed in the trough on 17 October. However this is not 
considered violation of homogeneity because of the extremely weak current speed in the 
trough on that day. Perusal of this data indicates the assumption of homogeneity for 
15- IS October is reasonable. Data sets from 15- IS October arc considered to meet the 
dual assumptions of homogeneity and stationarity. 
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Figure 19. H m0 ('x') and wave incidence angle ('o') versus time for 15 - 18 October 
1986. Numbers indicate data runs. 
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Tima E.S.T. (hr) Time E.S.T. (hr) 



VI. SURFACE WIND STRESS (tp 
A. INTRODUCTION 

The objective of this chapter is to determine representative surface wind stress for 
the surf zone based on atmospheric drag coefficients calculated from surface wind stress 
measured outside the surf zone and wind speeds measured by the sled anemometer inside 
the surf zone. Dr. S. SethuRaman and an investigative team from North Carolina State 
University calculated surface wind stress for the end of the FRF pier during SUPER- 
DUCK. Wind speeds were measured with a three-component Gill anemometer at a 
height of IS. 7 m above mean sea level (MSL) and wind stress was calculated using the 
eddy correlation method. Although this wind stress measurement was outside the surf 
zone, it is a simultaneous and uncontaminated measurement over water and within 5u0 
m of the sled transect areas. 

A synopsis of the synoptic scale meteorological conditions during October 1 9S6 is 
listed in Table 7. Sled deployments were on 11-12, 15- IS. and 21 October. Predominent 
meteorological features were onshore winds and an unstable atmospheric boundary layer 
(as indicated by the air-sea temperature difference in 'fable 7). 

The Gill anemometer has three polystyrene propellers which are rigidly mounted 
with their axes at 90° from each other. These propellers drive three d.e. generators 
whose output is directly proportional to wind speed. The anemometer has a threshold 
of 0.1 ms -1 . a distance constant of 0.S m. and can be sampled at 100 IIz. Although 
additional wind measurements were acquired at 22 m and 13.7 m, the eddy correlation 
method is selected over the profile method for wind stress determination due to its in- 
creased accuracy. Air temperature ( T a ) and relative humidity (RH) measurements were 
acquired by a YSI 440 IS Thermilinear thermistor and a slow response relative humidity 
sensor (Model 207) manufactured by Campbell Scientific, respectively. Both instruments 
were mounted at IS. 7 m above MSL. A YSI temperature probe measured sea temper- 
ature (T s ) at a depth of 0.6 m below MSL. 

Wind stress determined by the eddy correlation method is computed by: 

r = p a [iu'n-'f + (i-V/jT (51) 
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Table 7. SUPERDUCK SYNOPTIC-SCALE METEROLOGICAL DATA (OC- 





TOBER 1986) 


(from SethuRaman et al.. 


1987). 






Date 


Wind Field 


Brief Description 


T a ~T s 


Clouds 


Rema: 


1 


SW 3-7.5 


High dominates 


6.4 


Clr 




2 


Var. 3-7 


High dominates 


5.5 


Cu.As S-B 


Hazy 


3 


SE-SW 1-6 


High dominates 


1.7 


Ac.As S-B 




4 


SW 7-10 


High dominates 


7.9 


Clr 




5 


SW-NE 6-7 


Frontal passage 


3.3 


Cu,Sc S-B 




6 


XW-XE sust 13 


Post-frontal region 


-2.1 


Cu.Ac S 


Swell 


7 


XE-SE 3-5 


High dominates 


-3.5 


Ac.Cu, B-S 




8 


SE-SW 


High dominates 


2.1 


Sc.As B-0 




9 


XE-SE 1-5 


Cold front approach 


0.7 


Cu,Sc S-B 


Rain 


10 


XE 13- 16* 


Cold front passage 


-2.3 


St. Sc O 


Rain 


11 


XE 11 - 16* 


Post-frontal region 


-0.S 


St, Sc O-B 


Swell 


12 


XE 6-10 


High dominates 


-0.5 


Ac.As S 




13 


SE 3-8 


Cyclogenesis along shore 


1.0 


Cu.As B 


Rain 


14 


SW 6-S 


Cold front approach 


1.4 


TCu.St O 


Rain 


15 


XE 7-11 


Trontal passage 


-6.8 


St.Cu 


Swell 


16 


XE-XW 7-13* 


Post-frontal region 


-3.3 


St. Sc O-B 


Swell 


17 


XE-XW 4-S 


High dominates 


-2.3 


Cu.Ac S 




IS 


XW-XE 9-14* 


High dominates 


-3.7 


Cu S-B 


Swell 


19 


XE-XW 6-13* 


High dominates 


-2.9 


Cu S-B 




20 


XE 5-S 


High dominates 


-2.0 


Clr 




21 


XE-SE 1-3 


High dominates 


-1.3 


Cli- 




22 


SW 3-5 


High dominates 


1.7 


Ch- 




23 


SW 4-6 


High dominates 


2.5 


Cs S-B 




24 


XW-XE 3-9 


1 ligh dominates 


0.0 


Ac.Cs S-B 




25 


XE-E 7-10 


Warm front section 


0.0 


Sc. Ac O 


Rain 


26 


SE-SW 5-7 


Cold front approach 


4.1 


St.Cu O 




27 


SW-XE 3-5 


Cold front passage 


2.2 


Sc.As B-0 


Rain 


2S 


XE-SE 2-8 


High dominates 


-0.6 


Cu 




29 


SW 1-5 


Transition of 2 highs 


0.7 


Sc.Cs B-0 




30 


SW-XE 4-7 


Front sweeps through 


-1.2 


Sc.Cs B-O 




31 


XE gust 12 


High dominates 


-0.8 


Cu.Sc S-B 





Comments: 

1. S: scattered, B: broken, Clr: clear, O: overcast 

2. Cyclogenesis occurred in nearby region on 11th and 16th 

3. Strong upper- level trough accompanied cold front during 14- IS Oct 

4. * indicates high-energy cases of SUPERDUCK experiment 

5. T a — Tj was measured at 1400 EST which was the approximate mid-time 
of the sled transects. 
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where u.v.w denote wind speed in the horizontal and vertical directions and the primes 
(') denote deviations from the mean. Hourly means of wind stress data (SethuRaman 
et ah, 1 9 S 7 ) allow determination of an atmospheric shear velocity («,): 

".=( 77 )^ (52) 

where r is the surface wind stress. 

Inherent in the u, determination is the inclusion of atmospheric stability. The de- 
scription of the turbulent regime of a temperature-stratified friction layer is based on the 
Monin-Obukhov similarity theory, which defines the stability-dependent vertical wind 
speed gradient as: 



cu_ 

cz 



u. 




where z 
L 

K 

<i>m 



= height above the surface (m), 

= Monin-Obukhov length. 

= von Karman's constant (0.4), and 
= a stability function. 



(53) 



The non-dimensional ratio z L is a measure of atmospheric stability with negative 
values denoting unstable conditions and positive values denoting stable conditions, z L 



is delined as: 



L 



T v ut 



(54) 



where T v = virtual temperature (degrees Kelvin), and 
7^ = virtual temperature scaling parameter. 



T y is often considered to be the mean virtual potential temperature within the surface 
layer (Geernaert et al.. 19S7a) rather than its more strict definition of: 

T v = r.. ( 1 + .61</) (35) 

where q = specific humidity. In this research, T a in degrees Kelvin was substituted for 

Ty. 
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Integrating equation (53) from the surface (z = : 0 , where u = 0) to an arbitrary 
height (z = z. u = //,), an expression for a stability-dependent logarithmic wind velocity 
profile is obtained: 



= 



u* 



K 



( In 



(56) 



Determination of z L requires an iterative process, since 7 V * is also a function of 
stability. With all z L's initially set to zero, and the humidity roughness scale (z oq ) and 
temperature roughness scale (z ol ) set equal to 10 -6 , an iterative approach (Appendix C) 
based on the stability flux profile relation of Businger et al. (1971) is utilized to determine 
z/'L. and x ¥ q are assumed equal and formulated as (Large and Pond, 19S1): 



,v = (l - 16( -£-))« 

> P = 2 i n[ i,l± - l - l . ] JL <0 

V = ~7(f) f>0 



(57) 



The air at the air-sea boundary is assumed saturated. 

The stability function. </>„„ is described by the Businger-Dyer formula for unstable 
stratifications (Businger et al.. 1971) as: 

0 m = (l-16 -j-ff (58) 



The expression for T,,, is (Paulson, 1970): 

1 + 1 + 1 1 ~ 
Vm = 2 ln[ f2- ] + ln[ ] _ 2 + — 



(59) 



which is an empirical representation of comprehensive observations over smooth terrain. 

B. DETERMINATION OF AN ATMOSPHERIC DRAG COEFFICIENT FOR THE 
SURF ZONE 

Atmospheric shear velocity and L are independent of height in the surface boundary 
layer (Panofsky and Dutton. 1984) and are assumed to apply for the sled transect lo- 
cations in the surf zone located within approximately 500 m of the surface wind stress 
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measurements. Atmospheric shear velocity and L for the sled runs are linearly interpo- 
lated from the hourly pier data for the temporal midpoints of each sled run. 

Boundary layer investigators universally convert drag coefficients to be applicable 
at 10 m height to allow intercomparison between different investigations, l or equation 
(56) to be valid, a stability-dependent wind speed at 10 m ( L' 10 ) is required. L" 10 is ob- 
tained from equation (56) as: 

t\o = t*: m - [ ln( -jg- ) + ¥„( -J - ) - Vmi -77 }] (60) 



L' 10 is first calculated from the hourly pier data. Then two drag coefficients for a 
height of 10 m at the pier are determined — the stability-dependent Q 10 (from equation 
1 1 with l' l0 — U 0 substituted for W) and the stability-independent or neutral drag coef- 
ficient. Q,, 10 (from equations 1 1. 52, and 56): 



Qio — 



Pa (6 10 ~ 



1 0 





( 01 ) 



(62) 



The surface boundary condition ( i' 0 ). represented by a mean current, was taken into 
account to determine the wind speed relative to the moving water surface. i' 0 is assumed 
zero for the end of the pier. U 0 is represented by the mean surf zone current for the sled 
locations. 

Wind, wave height, atmospheric stability, and initial pier C d data versus time for the 
seven days of sled deployment are shown in Figure 20 and Appendix A (Figure 4-1 
through Figure 46). Drag coefficients and pier and sled wind speeds are adjusted to 10 
m elevation. Sled anemometer winds are represented by the first, mid. and last two- 
minute mean wind of the run. The clear arrowhead depicts the hourly mean wind vector 
for the pier wind and the solid arrowhead depicts a 34-minute mean wave direction vec- 
tor for a location near the pier end. The point of the arrowhead indicates the time for 
which the vector is valid. Wind and wave directions are relative to the shoreline which 
is represented by the y-axis. 
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arrowheads denote wind direction and closed arrowheads denote wave 



incident angle. Winds are adjusted to z = 10 nr. 
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The absence of hourly 10 L (z L stability indicator for z = 10 in) data and wind data 
indicates no boundary layer information was available. Absence of hourly 10 L data 
when a coincident wind speed is depicted indicates a 10 L value more negative than the 
limit plotted. 

Boundary layer stability is a primary concern in interpreting the wind data. There 
were no stable atmospheric conditions during SL’PERDL’CK. Q„ ]0 was plotted so as to 
remove the effect of stability in an effort to relate the drag coefficient variation to an- 
other physical parameter, such as wave height, wind direction, or wind speed. Q„ 10 is 
the standard drag parameter utilized for wind stress comparison between different in- 
vestigators and serves as a basis for examining wave influence on wind drag. Neither 
wave height, wave direction, nor wind direction exert obvious influence on the drag co- 
efficient as shown in Figure 20 and Appendix A (Figure 44 through Figure 46). 

Increased instability places the atmosphere in a convective instability or local free 
convection regime where buoyancy forces become equal or greater in magnitude than 
mechanical mixing forces. Wyngaard and Cote (1971) conducted measurements of ver- 
tical velocity and temperature under unstable conditions and observed local free con- 
vection to occur at z L. < -0.5. Using the Businger et al. (1971) stability flux profile 
relation. Garrett 1 1977) showed that as instability increases, not only does Q increase, 
but roughness becomes significantly more important. The Businger ct al. (1971) iterative 
approach is inapplicable in a local free convection regime. Therefore an empirical ap- 
proach is deemed more appropriate for surf zone Q determination under local free 
convection. 

Sled runs of z 1. < -0.5 generally occurred during periods of low onshore wind speed. 
Therefore the drag coefficients in Figure 20 and Appendix A (Figure 44 through 
Figure 46) depicted for low wind speeds and very unstable conditions are inaccurate. 
Onshore winds of low speed and diabatic (unstable) conditions are representative of a 
sea breeze. A sea breeze is a local coastal wind caused by the temperature difference 
when the sea is colder than the adjacent land. It usually occurs on relatively calm and 
sunny days (Iluschke, 1959). Review of surface weather analysis maps (Appendix D) 
and wind speed, wind direction, and stability (Table 7 and Table 9) indicate sea breeze 
conditions for most of the sled data run periods. 

Measurements of surface wind stress above the swash zone under the effect of a sea 
breeze were made by simultaneous wind and temperature profiles over a shore transect 
near Fort Walton Beach. FF (Hsu, 1970). The sea breeze in the surface boundary layer 



68 



was determined to be in the atmospheric free-convection regime. A least-squares fit to 
the data (Figure 21) yields: 

Q i0 X 10 3 = 10.06 - 1.30 (Z7 10 - U 0 ) (63) 

where L' 10 is the mean wind speed at 10 m above the sea surface. It is noted that the 
drag coefficient under the sea breeze condition is higher than that under more stable 
conditions which is in agreement with e.g.. Roll (1965). 




Figure 21. An empirically-derived stability-dependent drag coefficient obtained un- 
der a sea-breeze condition (from Hsu. 1970). 

To determine appropriate drag coefficients for the sled transects based on the 
SUPERDL'CK wind stress data, the following procedure is employed: 

• When the sled's anemometer failed (Runs 1S-1. 18-2. and all runs on 21 October), 
a linearly-interpolated hourly mean i) 0 from the pier is used as the sled L\ 0 . 

• A linearly-interpolated C d]0 from the pier is assigned for the temporal mid-point 
of the sled run. 

• When the actual z L was more negative than -1 or the wind speed was less than 
6.5 ms -1 , Q ]0 is determined from the empirical relation of equation (63) where the 
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sled's L] 0 is determined from equation (60) and the z represents the run-varying 
height of the sled anemometer above mean water level for that particular run. 

• When z L was more negative than -1 and the wind speed greater than 6.5 ms -1 , 
Q l0 is linearly interpolated between the value preceding and following the one in 
question. This is felt to be a better approximation than not using stability. 

C. ADJUSTMENT OF Q FOR THE SURF ZONE ENVIRONMENT 

Since the surface wind stress was measured outside the surf zone, a Q adjustment 
is required to more realistically reflect the surf zone environment. Early investigators 
of oceanic Q's sought a dependence on wind speed. Subsequent investigators (e.g., I lsu. 
1974. 19S6a; Byrne, 1982: Gccrnaert et ah. 19S6) began to search for a sea state de- 
pendence for Q, since surface gravity waves are assumed to act as roughness elements 
(z 0 ) for the momentum transfer process. Most of these investigations were done for open 
ocean conditions and therefore are not directly applicable to the surf zone. 

1. Uniqueness of surf zone air-sea interaction 

T he uniqueness of surf zone air-sea interaction as opposed to open ocean air-sea 
interaction can be summarized by three factors -- shoaling of the bottom, wave refrac- 
tion. and wave breaking. 

a. Shoaling effect 

Momentum transfer is based on the relative wind speed difference between 
the absolute wind and wave speeds. During shoaling, there is a wave steepening and 
celerity decrease (relative wind increase) of surface gravity waves as they shoal. 
Gccrnaert et al. (1987b) state that Q increases with decreasing values of the ratio c 0 m, 
where c 0 is the phase speed of the dominant locally-generated grav ity wave. Thus, larger 
than expected drag coefficients are observed over the sea during shoaling conditions 
(Smith. 1980: Gccrnaert et al., 19S6). 

b. Refraction 

In most cases, the wave propagation direction is at some angle (up to 
180°) to the predominant wind direction due to refraction of the waves in shoaling water. 
Over the open ocean. C d is found to be larger if the wind is blowing at an angle different 
than the wave direction (Large and Pond, 1981: Smith. 1980). Geernaert et al. (19S6) 
proposed wind stress enhancement in advance of a cold front due to surface gravity wave 
convergence. Swell generated behind the cold front would travel through the front and 
interact with waves generated ahead of the front and propagating parallel to the front. 
This interaction would result in a chaotic sea of larger wave amplitudes, steeper gravity 
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waves, and more wave breaking -- very similar to a robust surf zone. This enhancement 
was observed by Davidson et al. (1988a, 19SSb) during the FASINEX 1986 experiment. 
c. IF are Breaking 

Moderate to high surf conditions undeniably result in increased surface 
roughness, which would result in increased drag coefficients. Although some work 
(Amorocho and DeVries, 1980; Melville, 1977) has quantified relations between breaking 
waves and increased drag coefficients, these investigations have dealt with the open 
ocean where the onset of breaking waves is produced by a dramatic increase in wind 
speed. 

2. Surf Zone C d Adjustment Formulation and Results 

The wind stress data from the end of the FRF pier is assumed to include both 
the shoaling and refraction effects on C d as one transitions from the open ocean to the 
surf zone. The increased C d effect due to shoaling from the end of pier to the surf zone, 
as well as the additional refraction 0(15 — 25° maximum), is considered to be small 
compared to the breaking wave effect. 

The approximate mean distance from the swash line to the end of the pier is 500 
m. Each day at 0700 EST, the distance from the mean swash line to the outer extent 
of the surf zone was "stepped off' along the pier. These rough values for surf zone width 
on sled deployment days are listed in Table 8. The surf zone never extended to the end 
of the pier, although on 1 1 October, the farthest surf zone extent nearly did. With 
onshore winds, the increased roughness would not take effect until one was a 
measureable distance within the surf zone. Therefore, the increased roughness due to 
the breaking waves of the suif zone is not reflected in the end of the pier wind stress 
data. 

The purpose of this research is not to explicitly characterize the C d for the surf 
zone, but to intelligently select an appropriate C d with which to address longshore sur- 
face wind stress as measured by a single anemometer in the surf zone. The increased 
surface roughness due to surf zone breaking waves may be conservatively likened to the 
chaotic seas observed by Davidson et al. (19SSa, 1988b) in advance of storm fronts. 
From a composite of wind stress measurements taken from multiple chaotic sea cases 
found in advance of storm fronts, the drag coefficient was found to be increased by ap- 
proximately 33% over that expected for non-chaotic seas and equal wind speed 
(Davidson et al., 19SSa. 1988b). Therefore an increase of 33% in the drag coefficient 
(determined from equation 61) was deemed appropriate for all runs except those runs 
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Table 8. SUPERDUCK SURF ZONE WIDTH FOR SLED DEPLOYMENT 
DAYS (from Crowson et al., 1988). 



October 


Surf zone width Cm) 


11 


500 


12 


320 


15 


120 


16 


1 00 


17 


80 


18 


80 


21 


100 



occurring under strong diabatic conditions where equation (03) was utilized. Hsu's 
(1970) relationship is empirical and therefore the Q values determined by equation (63) 
would already have the increased roughness considered. Drag coefficient values for each 
run and procedural step are listed in Table 9. All the values are of O(10~ 3 ). 

Surface wind stress is determined from equation (11) using the C/s listed in 
Table 9. Atmospheric density is assumed constant at 1.25 kgm -3 and 0 is obtained from 
the end-of-pier wind direction. U I0 — U 0 is substituted for W where U 10 is the stability- 
dependent mean windspeed measured by the sled anemometer for the entire data run 
(approximately 34 minutes), corrected to an elevation of 10 m. U 0 is the mean surf zone 
current measured by the sled current meters. The term L' ]0 — V 0 is obtained from 
vectorial subtraction. 
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Table 9. ATMOSPHERIC DRAG COEFFICIENTS (C„ I0 x 10 3 ) FOR EACH 
SLEI) DATA RUN 



Sled 


Sled 


Pier 






^ler 


'Sled 


3 Sled 


Run 


^ 10 — Q 


^10 ~~ ^ o 


C. 


10 

l 


Qio 


Qio 


Qio 




(ms -1 ) 


(ms -1 ) 


(ms -1 ) 










111 


14.44 


11.13 


0.56 


-0.025 


2.49 




3.31 


112 


13.69 


11.33 


0.57 


-0.023 


2.53 




3.36 


113 


13.2S 


11.25 


0.55 


-0.023 


2.41 




3.21 


121 


10.70 


7.28 


0.30 


-0.06S 


1.72 






122 


9.62 


8.15 


0.32 


-04*48 


1.54 




2.05 


123 


9.14 


7.86 


0.29 


-0.055 


1.41 




1.S8 


125 


6.3-1 


4.79 


0.15 


-0.256 


1 .07 


1.82 




126 


5, N N 


4.34 


0.16 


-0.230 


1.48 


2.85 




127 


5. S3 


4.36 


0.09 


-0.640 


0.49 


2.4S 




1-1 


9.49 


8.78 


0.30 


-0.337 


1.18 




1.57 


152 


8.00 


7.76 


0.24 


-0.600 


0.95 




1.26 


153 


7. S3 


7.14 


0.27 


-0.44 S 


1.44 




1.92 


154 


7.57 


6.33 


0.14 


-1.410 


0.55 




1.92 


162 


12.51 


10.99 


0.37 


-0.1 2S 


1.15 




1.53 


163 


11.14 


11.30 


0.39 


-0.1 16 


1.22 




1.62 


164 


10.42 


11.25 


0.39 


-0.123 


1.22 




1.62 


165 


9.76 


10.72 


0.39 


-0.130 


1.33 




1.77 


166 


9.89 


9.67 


0.36 


-0.158 


1.40 




1.S6 


171 


6.25 


6.37 


0.20 


-0.473 


1.08 


1.94 




172 


7.21 


4.77 


0.20 


-0.471 


1.31 




1.74 


173 


7.56 


3.11 


0.21 


-0.4"6 


1.35 




1.S0 


174 


6.50 


5.27 


0.21 


-0.476 


1.39 




1.85 


175 


5.S7 


5.55 


0.18 


-< *.464 


1.10 


2.43 




176 


5.63 


5.18 


0.1 1 


-1.414 


0.49 


2.74 




177 


5.S6 


4.85 


0.08 


-2.505 


0.3O 


2.44 




17S 


5.86 


4.54 


o.os 


-2.664 


0.33 


2.44 




1S1 


in. 34 


10.61 


0.35 


-0.1 63 


1.11 




1.4S 


182 


10.39 


10.87 


0.40 


-0.13O 


1.34 




1.78 


183 


11.74 


11.20 


0.39 


-0.140 


1.17 




1.56 


184 


12.17 


11.42 


0.40 


-0.127 


1.20 




1.60 


1S5 


13.60 


11.50 


0.42 


- 0. 1 1 1 


1.32 




1.76 


186 


13.94 


11.80 


0.43 


-0.091 


1.32 




1.76 


211 


3.94 


3.96 


O.OS 


-2.147 


0.46 


4.94 




212 


4.49 


4.24 


0.15 


-0.72S 


1.20 


4.22 




213 


3.66 


3.5 1 


0.14 


-0.690 


2.10 


5.30 




214 


2.26 


2.24 


0.13 


-0.775 


3.38 


7.12 




215 


2.36 


2.25 


0.08 


-2.068 


1.25 


6.99 




216 


2.90 


2.89 


0.03 


-14.286 


0.12 


6.29 




217 


1.94 


1.93 


0.04 


-6.645 


0.49 


7.54 




2 1 S 


1.61 


1.57 


O.05 


-2.939 


1.38 


7.97 




219 


1.56 


1.39 


0.07 


-1.67S 


3.69 


8.05 





'eqn. (61) : eqn. (63) 3 33% increase due to surf zone roughness 
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VII. LONGSHORE MOMENTUM BALANCE 



The longshore momentum balance is calculated based on field measurements of 
wind, waves, and currents acquired across a barred beach. Spatially-dependent bed 
shear stress coefficients and the relative importance of the various terms in the momen- 
tum balance are calculated. Wind force to wave force ratios arc calculated from the field 
data and compared to the results of the wind forcing model in Chapter III. This chapter 
concludes with an error analysis of these calculations. 

The slope array did not perform as required for determination of the wave-induced 
longshore momentum flux (S yx ) in equation (1). Pre- and post-calibration values for the 
slope array differed as much as 10%. Assuming that breaking wave turbulence outside 
the surf zone is negligible, the gain required to match the slope array's radiation stress 
(S vx (fj) to the current meter's radiation stress (S yx (f)) was 0(30 - 60%). Therefore, slope 
array data were judged in error and were not used in this analysis. Studies addressing 
turbulent momentum flux were suspended. 

The slope array loss necessitates a rewriting of the longshore momentum balance 
equation as: 



c.\L 



ci 



cSj 



cx 



_J> I 
L y + L >* 



(64) 



The turbulent momentum flux gradient is combined with the wave-induced momentum 
.flux gradient to form a combined momentum flux gradient as measured by the current 
meters. The S V J term implies integration over frequency (equation 80). Determination 
of surface wind stress is addressed in the Chapter VI. Calculation of the remaining terms 
in equation (64) will be addressed separately in this chapter. 

A. TEMPORAL, OR ACCELERATION TERM ( dM y /dt ) 

The time-averaged depth-integrated mean momentum per unit area in the longshore 
direction (M y ) is composed of a steady (M v ) and unsteady (M y ) flow component: 



My = My + My' = pV{t, + 




(65) 
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The steady component is calculated every 64 seconds with the term, ?/ + h , being taken 
as the mean water level (MWL) measured by the pressure sensor for that period of time. 
This resulted in thirty-two mean values of 64-second duration for M y . 

The unsteady flow contribution to the mean momentum, MJ, is attributed to the 
mass flux of the wave motion: 



M; = JJ p v' dz (66) 

where v' is only known at z = 0 from wave theory. By expanding in a Taylor series 
about z = 0 and retaining only the first order term, Md may be expressed independent 
of height in the crest-trough region: 



Aly = J 0 P k(0) dz — p k(0) tj 

Using linear theory and substituting for v'(0) and i/: 



m; = 



Pga 



_i_ 

c 



A 

k 



E_ A 
C k 



(67) 



(6S) 



where a is wave amplitude. Rewriting equation (68) in terms of a finite frequency band, 
a significant wave incident angle (£(/)) , total wave energy (E(f)), and wave celerity (C(f)) 
yields: 



My ' = t "iw (69) 

where 

C(/) = tanh kh (70) 

Throughout this analysis, wave number (k) is calculated non-iteratively by the formu- 
lation of Wu and Thornton (19S6). 

Spectra of "significant wave angle" (Appendix F) are calculated as (Higgins et ah, 
1981): 
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(71) 



sin «(/) cos o.[f) = 



Sy x {f) 

[£(/)*(/)] 



The energy spectrum is calculated from velocity measurements. Time-averaged 
depth-integrated kinetic energy (KE) may be written as: 



KE = 0.5 pj h (u 2 + v 2 + \ 2 )dz (72) 

where w is the component of the flow velocity in the vertical direction. To the first order 
of approximation, there is equipartition of energy between potential energy and kinetic 
energy for linear progressive waves (Kinsman, 1965). The total energy, E(f), is calculated 
as twice the kinetic energy: 

E(C) = 2\H l (f)\ 2 [G u (f)+G v (m (73) 

where 

Pi sinh(2A/i)] 

(/ 4 ) 

4A[ cosh‘A(/i + zj] 



l^iOOl 



and G u and G v are the one-sided spectral density functions for the horizontal velocity 
components in the x and y directions and z m is the measurement depth. Calculations 
of equations (69) through (73) indicated that the time-varying momentum of the waves 
for this data set was one-to-two orders of magnitude less than the time-varying mo- 
mentum of the steady flow. Therefore the temporal term may be simplified as: 



cM y A M v 
3t ' Ai 



(75) 



where 

MyStJIy = pV(7j + h) (76) 

and At = time step over which the finite difference is calculated, and 
tj + h - MWL as sensed by the sled's pressure sensor. 



A value of XL is first determined for each ensemble of 64 seconds in the data run. 

y 

Since the temporal scale of the longshore current is unknown, a range of ensembles 
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(comprised of 1 to 1 1 segments of 64-second duration) are defined using a running mean 
interval, At'. Then a finite difference is calculated from these larger ensembles for a 
specified At. Finally, a mean temporal term is calculated: 




(77) 



for the k increments of AM y IAt. This mean temporal term is assigned as the temporal 
term for the sled run. Calculations of the mean temporal term (equation 77) with and 
without a running mean, as well as variation of the time step (with subsequent change 
in the k increment), made negligible difference for this data. Therefore At = 64 seconds 
is used. The selection of At' is discussed in paragraph D below. 

B. RADIATION STRESS GRADIENT (dS yx T / ex) 

Longshore momentum flux, as defined in equation (2), is calculated from current 
meter measurements and applying linear wave theory: 

SyJJ) = \H 2 (J)\ 2 C uv (J) (7S) 



where | H 2 (J) 1 2 is a linear transfer function: 
I |2 pt sinh(2/v7i) + 2kh~\ 

I HyJ) I — 

4A'L cosn k{h + z m ) J 



(79) 



and C uv is the real part of the cross-spectrum of u' and v\ It is assumed in the applica- 
tion of the transfer function (equation 79) to the measured covariance that the vertical 
distribution of turbulence in the water column is the same as in the waves, and is de- 
scribed by the transfer function which is a function of depth and frequency 
(Figure 22). 

The total radiation stress (S y J), where total indicates frequency summation, is de- 
termined from: 

h 

- y>>-(/)A/ (SO) 

A 
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Figure 22. Linear wave theory transfer function used to calculate S yx {J) from cur- 
rent velocity, showing decreasing frequency dependence with depth 
across a barred beach. 
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where /, and f 2 represent the lower and upper frequency limits of the sea-swell band un- 
der investigation and Af represents the resolution of the spectral estimates. is the 
quantity required for the momentum balance of equation (64). 

Analysis of the longshore component of momentum flux (S yx (f)) was accomplished 
using Fast Fourier Transforms and ensemble averaging. Ensemble averaging requires 
segmentation of a time series into several shorter time series of equal length, and calcu- 
lation of a sample spectrum for each of these sub-series. The mean of these spectra is 
a smoothed spectral estimate. The advantages of ensemble averaging are minimal com- 
puter storage and computational ease. 

The data were initially demeaned and detrended. The 34.1 minute data runs were 
segmented into 32 sub-series of length 64 seconds. The S yx (J) used in equation (SO) is 
the Bartlett smoothed spectral estimate (Jenkins and Watts, 1968) from 32 sample 
spectra. All spectral estimates ( S yx , > 7 , etc.) have 63 (v = 2m - 1) degrees of freedom, a 
spectral resolution (A/) of 0.015 Hz, and Chi-squared 95% confidence limits equal to 
0.73 and 1.45 times the spectral estimate itself (Jenkins and Watts, 1968). Since sea and 
swell are of interest, lower and upper frequency limits are established at 0.06 and 0.44 
Hz, respectively. These limits are used whenever frequency summation was required for 
any parameter. 

Radiation stress spectra may be negative at one frequency while positive at another 
frequency, as in the case of sea and swell approaching from different quadrants. S yx (J) 
spectra transitioning from positive, through zero, to negative, and vice versa, may not 
be accomodated with a log-linear plot, therefore S yx (f) spectra are plotted on linear plots. 
Related spectra are also plotted in this manner for consistency. 

The calculation of the gradient of radiation stress requires numerical differencing 
over an appropriate length scale: 



cS yx ^ A S yx 
cx ~ A.v 



(81) 



The spatial step (Ax) was initially designated as the distance between the front and rear 
current meters after adjustment for sled yaw angle. A radiation stress gradient (equation 
81) determined with this definition of A.v is labeled to have been determined by the "lo- 
cal" method. The A.v distance on the sled was approximately 2 m and the "local" finite 
differencing scheme is: 
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( 82 ) 



<?S 



yx 



(Syl) r 



cx 



- (S v ,) 



vx /front c m 



Ax 



A second approach to the determination of the radiation stress gradient, called the "in- 
tegrated" method, utilized the distance between sled run positions for Ax. This distance 
was 0(20-40 m) and this finite differencing scheme is: 



6S V 



ox 



4 



(Sj\ - (S V J) 4 



Ax 



+ 



1-4 



(V 7 ). - (s,I)j 

A- V l-3 



+ 



(s/) 3 ~(sju (Sj,-(S„ T ) , ] 

7 + - Z - r (83) 

Ax 2 _4 ^- x 2-3 

where the numbers 1 through 4 denote the two pair (front and rear) of current meter 
data available from two sled run positions (see Figure 23). The estimate of the S y J 
gradient is based on four independent current meter locations covering slightly different 
transect segments and slightly different Ax's. These segment and Ax differences are 
generally less than 10% of the distance between sled mast locations. It is therefore as- 
sumed that the four estimates represent four sample means of the true S y J gradient be- 
tween the two sled locations. From the Central Limit Theorem, the standard deviation 
of sample means is written as: 




or o asJ 

Ax 



a as t , T 



dx 




(84) 



where n is the number of sample means. Thus the standard deviation of the sample 
mean is reduced by 50% with the use of four sample means rather than one. This fact 
will be used in determining the measurement error for 5S yx T / cx in paragraph E. Se- 
lection of the "local" method versus the "integrated" method proved to be of significant 
importance. 

C. SEA SURFACE ELEVATION (»/) 

Although not specifically required for the momentum balance, sea surface elevation 
spectra provided a useful reference against which to compare radiation stress spectra. 

Sea surface elevation spectra (G^(/) ) are determined from linear wave theory using 
data from both the sled's pressure sensor and current meters: 
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g,(J) = um \ 2 G Ph (j ) 



(85) 



where 



\Hi(f )\ 2 



cosh 2 (A7t) 
cosh" C A (A + zj: 



( 86 ) 



and G ft (/) is the one-sided energy density spectrum for the pressure head (/y,) measured 
by the pressure sensor. The pressure head is determined from the hydrostatic equation: 



Ph = 



P ~ Pam 

Pg 



(87) 



where p is the absolute pressure sensed by the pressure gage and p alm is the atmospheric 
pressure. 

The formulation for determining G, ? (/) from current meter velocities is: 

G,jf) = 1 // 4 OO 1 2 IGJJ) + GM (88) 

where 



I !Uf) 1 2 



sinh(2A/fi 

2gk cosh 2 [A(// + zj] 



(89) 



G r {J) produced by these two formulations are in good agreement for the three 
vertically-stacked current meters and the pressure sensor, thereby providing confidence 
for the >7 spectra. Sea surface spectra produced in Appendix F are formulated from 
equation (S5). 

Wave height parameters are calculated from pressure-derived ?/ spectra: 

Hi ms = (8 ofi (90) 

H mo = 4 ( 91 ) 

where the total variance is expressed as: 

h 

A 
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D. RESULTS OF THE MOMENTUM BALANCE 

Field measurements of wind, waves, and currents are used to calculate the local 
longshore momentum balance at various locations across a barred beach. From the 
momentum balance, bed shear stress coefficients and the relative importance of the 
various terms comprising the balance are calculated. Mean values are determined rela- 
tive to the bar. Wind force to wave force ratios are calculated from the field measure- 
ments and compared to the results of the wind forcing model in Chapter III. 

I. Bed Shear Stress Coefficient 

A spatially-dependent bed shear stress coefficient (cf) may be determined from 
equations (18) and (64) as: 



c f - 



„ c T 

oS vx 



ox 



+ t; - 



8.\fy 

Cf 



(93) 



Resultant spatially-dependent bed shear stress coefficients are listed in Table 10 and il- 
lustrated in Figure 24. The column titled "cfjf are the cf s determined utilizing the "lo- 
cal" method for the radiation stress gradient. These coefficients are representative for 
the sled run position itself. The surface wind stress and temporal variability terms were 
34.1 -minute means taken only from the sled run being investigated in anticipation of 
gaining insight into the temporal variability of all terms in the longshore momentum 
balance. A running mean of 12.S minutes was selected for determination of both the 
temporal term and the radiation stress based on radiation stress statistics of Freilich and 
Pawka (19S7). They showed that, given n independent realizations of a time series of 
u(t) and v(t), that the Gaussian approximation for true S vx statistics is valid for n > 11 
independent realizations. An ensemble length of 64 seconds and a minimum number of 
ensembles (12) was selected for the analysis. Unfortunately, this procedure resulted in 
significant cfl) variability. Utilization of a 34.1-minute mean dS yx T / 5x made negligible 
difference in resultant cfl). Two possible reasons for the large cf!) variability may be 
that the radiation stress is much more variable across the surf zone on a small spatial 
scale than previously imagined and the Ax is too small to resolve the difference between 
two large numbers (see equation (S2) and paragraph 2 of the Error Analysis section in 
this chapter). 

The column titled "cfi)" are those cf s determined using the "integrated" method 
for the radiation stress gradient. These coefficients are representative for the portion of 
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the bathymetric profile between sled run positions. The wind stress and temporal vari- 
ability terms are the 34.1 -minute means from the two appropriate sled run means. The 
"beach" position along the transect is defined as the intersection of the bathymetric 
profile with the MWL as determined by the sled pressure sensor on the run closest to 
the beach. Extremely weak longshore currents inshore of run 17-6's location explain the 
variable Cj(i)'s. They are disregarded. The bed shear stress coefficients determined by the 
"integrated" method exhibit less variability than those from the "local" method and are 
all of O(10 -3 ). The "integrated" method is considered to be a reasonable method for the 
determination of the radiation stress gradient. 

2. Relative Contributions of the Longshore Momentum Balance Terms 

Relative contributions of wave forcing (8S y J / 8.x), wind forcing (r^), and 
temporal variability of mean momentum (8M V / cl ) to the total momentum balance are 
determined by first defining their total contribution (TOT) as the cumulative sum of 
their absolute values: 



TOT 




c.x 



+ 




BAty 

cl 



(94) 



Their individual relative contributions (rc) are then determined by dividing by the total 
contribution: 



wave rc 



cS yx T 

c.x 

TOT 



wind rc 



TOT 




6 My 

Cl 

TOT 



(95) 

(96) 



(97) 



Relative contributions from the wave, wind, and temporal terms, as formulated 
in equations (95) to (97) and utilizing the "integrated" method, are calculated. Results 
are presented in order of consecutive runs (Table 1 1) and by relative location along the 
transect with respect to the bar (Table 12). Relative locations are delineated as "off- 
shore" of the bar, "on the bar" and immediately before the bar, and in the "trough" re- 



84 



Table 10. SUPERDUCK BED SHEAR STRESS COEFFICIENTS (OCTOBER 
15-18, 19S6) 



October 15. 1986 






October 17. 


1986 




Run 


cj{0 


Cj( 0 


Run 


9(0 


9(0 


15-1 


* 


0.006 


17-1 


0.002 


0.003 


15-2 


0.008 


0.002 


17-2 


0.003 


0.005 


15-3 






17-3 


0.004 


0.003 


15-4 


0.012 


0.001 


17-4 


0.026 


0.005 


Beach 






17-5 


0.099 


0.003 








17-6 


0.119 




October 16. 19S6 






October 18, 


19S6 




Run 


Cj{ 0 


Cj(l) 


Run 


9(0 


9(0 


16-2 


0.001 


0.003 


1S-1 


0.001 


0.003 


16-3 


0.004 


0.004 


18-2 




0.004 


16-4 


0.003 


0.001 


1S-3 


0.007 


0.001 


16-5 


0.001 


0.001 


18-4 


0.043 


0.001 


16-6 


0.002 


0.002 


18-5 


-.001 


0.001 


Beach 






Beach 







• indicates cy value < 5. x 10 -J 

• Cj([) are determined by the "local" method 

• Cj{i) are determined by the "integrated" method 

• relative contributions determined by the "integrated" method 



gion. The mean wind stress remained nearly constant throughout each transect of the 
surf zone. However, the temporal term, as averaged over 34.1 minutes, was variable in 
sign and magnitude. The location of maximum wave breaking along the transect can 
be identified (Table 11) where wave forcing's relative contribution is maximum. Two 
peaks in wave forcing correspond to two breaking areas - waves breaking offshore of the 



85 



bar and waves reforming and then breaking on the nearshore beach face. The relative 
contribution of wind stress inside the surf zone generally increases in value as one pro- 
ceeds onshore from the location of maximum breaking, due to the decrease in wave 
forcing. The relative contribution of wind stress can then decrease again inshore at the 
location where waves re-break on the nearshore. The relative contribution of the 
temporal term appears to be largest in the trough where wave forcing is reduced and 
possibly where the effect of longer period waves (e.g., surf beat or edge waves) is more 
pronounced. 

The mean relative contributions for the wave, wind, and temporal terms (0.86, 
0.0S, and 0.06, respectively) are essentially identical for both "offshore" and "on the bar" 
(Table 12). However in the "trough", the mean relative contributions (0.64, 0.17, and 
0.19, respectively) are significantly different. The relative importance of wave forcing is 
reduced by 25%, whereas the relative importance of surface wind stress is doubled, and 
the temporal term relative importance is tripled. 

A comparison is made between the wind force to wave force ratios produced by 
the wind forcing model in Chapter 111 to those determined for the SUPERDUCK data. 
The wind to wave force ratios in Table 3 and Table 4 are based on a planar beach, 
therefore only the olfshore side to the SUPERDUCK bar (tan/? ~ 0.03) is considered. 
The SUPERDUCK force ratios are obtained from ratios of the relative contributions 
of surface wind stress to radiation stress gradient in Table 12 for the "offshore" and "on 
the bar" categories. The wind force to wave force ratio for 15-18 October field data is 
0.09. The mean wind speed was 9.6 ms -1 and the mean wave incident angle for the 
"offshore" region was 12°. Using linear interpolation in Table 3 and Table 4 results in 
a wind to wave force ratio of 0.1 2. Thus the field measurements are in good agreement 
with the results of the wind forcing model in Chapter 111, with differences attributable 
to non-fully arisen seas and the presence of swell from a different direction than the wind 
waves. 

Bed shear stress coefficients are calculated for cases without the surface wind 
stress and without both surface wind stress and the temporal term. The percentage 
change in cy ranges from near zero to 100%, with a mean change of 10-30%. Thus under 
certain conditions for this data set, the surface wind stress and temporal terms can be 
of first order in the momentum balance. 

A mean cj based on relative postion along the transect with respect to the bar 
is determined. Bed shear stress coefficients determined from the "integrated" method are 
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Table 11. RELATIVE CONTRIBUTIONS OF THE LONGSHORE MOMEN- 
TUM BALANCE TERMS FOR SUPERDUCK 



Run 


6S yx T 1 dx 


T*7 

y 


cMyfdt 


15-1 


0.97 


0.02 


0.01 


15-2 


0.93 


0.04 


0.04 


15-3 


0.77 


0.17 


0.06 


15-4 


0.79 


0.14 


0.07 










16-2 


0.87 


0.08 


0.04 


16-3 


0.95 


0.03 


0.02 


16-4 


0.77 


0.19 


0.03 


16-5 


0.70 


0.18 


0.12 


16-6 


0.81 


0.02 


0.17 










17-1 


0.63 


0.17 


0.20 


17-2 


0.94 


0.03 


0.04 


17-3 


0.88 


0.03 


0.09 


17-4 


0.89 


0.05 


0.07 


17-5 


0.72 


0.16 


0.12 










1S-1 


0.83 


0.13 


0.04 


18-2 


0.93 


0.05 


0.02 


1S-3 


0.S1 


0.12 


0.07 


1S-4 


0.46 


0.33 


0.21 


1S-5 


0.21 


0.19 


0.60 


• Relative contributions determined by "integrated" method. Run number indi- 
cates offshore run of the required run pair. 

• Runs 17-6 through 17-8 had negligible currents of variable direction. Their cy 
results are disregarded. 
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Table 12. RELATIVE IMPORTANCE OF MO- 
MENTUM BALANCE TERMS FOR LO- 
CATIONS RELATIVE TO THE BAR 



Offshore 


Run 


cS vr T j d.v 


t? 


cMJdt 


15-1 


0.97 


0.02 


0.01 


16-3 


0.95 


0.03 


0.02 


16-2 


0.87 


0.08 


0.04 


17-1 


0.63 


0.17 


0.20 


17-2 


0.94 


0.03 


0.04 


18-1 


0.83 


0.13 


0.04 


18-2 


0.93 


0.05 


0.02 


Mean 


0.87 


0.07 


0.05 


On the Bar 


Run 


c'S rr r 1 c.x 


n 


cftljct 


15-2 


0.93 


0.04 


0.04 


16-4 


0.77 


0.19 


0.03 


17-4 


0.89 


0.05 


0.07 


17-3 


0.88 


0.03 


0.09 


IS- 3 


0.81 


0.12 


0.07 


Mean 


0.86 


0.09 


0.06 


Trough 


Run 


cS vr 7 / c x 




cMJct 


15-4 


0.79 


0.14 


0.07 


15-3 


0.77 


0.17 


0.06 


16-5 


0.70 


0.18 


0.12 


16-6 


0.81 


0.02 


0.17 


17-5 


0.72 


0.16 


0.12 


18-4 


0.46 


0.33 


0.21 


18-5 


0.21 


0.19 


0.60 


Mean 


0.64 


0.17 


0.19 



assigned and a mean determined for each of these relative locations (Table 13). The 
results indicate a decrease in mean Cy as one proceeds from offshore to onshore with 
mean cy values equal to 0.004 ("offshore"). 0.002 ("on the bar"), and 0.001 (trough"). 



8S 




The shoreward decreasing Cj was initially thought to be the result of decreasing bed 
roughness. As previously stated, bed roughness is a function of both sand characteristics 
and bed form. Beaches generally exhibit a shoreward gradation of sand grain size from 
large to small. Waves tend to sort the sand with the largest sand residing in the most 
energetic area of the surf zone. This gradation would be commensurate with the maxi- 
mum Cj found just inside the surf zone and a decreasing Cj further onshore. Unfortu- 
nately, field measurements of bed roughness were overlooked during the SUPERDUCK 
experiment and cannot be used to substantiate this hypothesis. Interviews with 
SUPERDUCK divers were conducted several months after the experiment, however 
their information about bed roughness was not specific enough to draw conclusions. 
An October 19S2 sampling of sand grain size across the FRF beach showed shoreward 
increasing grain size (Birkemeier et ah, 19S5). It is reasonable to assume that grain size 
gradation would not change over a few years, therefore bed roughness is discarded as the 
reason for the Cj- trend. Field and laboratory studies have shown that bed shear stress is 
enhanced when the near-bottom How field is composed of mean current and oscillator.' 
(i.e. wave) components (Grant and Madsen, 1979, 19S2; Grant et ah, 19S4). Thus 
shoreward decreasing wave action would result in shoreward decreasing Cj-. 



Table 13. MEAN BED SHEAR COEFFICIENTS DETERMINED FOR LO- 
CATIONS RELATIVE TO THE BAR 



Trough 




On the Bar 




Offshore 




15-4 


.001 


15-2 


.002 


15-1 


.006 


15-3 


<.001 


16-4 


.001 


16-3 


.004 


16-5 


.001 


17-4 


.005 


16-2 


.003 


16-6 


.002 


17-3 


.003 


17-1 


.003 


17-5 


.003 


18-3 


.001 


17-2 


.005 


1 S — 1 


.001 






18-1 


.003 


18-5 


.001 






18-2 


.004 


Mean Cj= 


.001 




.002 




.004 



The variation in Cy due to varying At' in the temporal term was examined by 
analyzing sample groups where At' was specified as a running mean of (approximately) 
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1, 5, 10, 12, 14, 18, and 22 minutes. The mean variation in Cj due to varying A t' for the 
SUPERDUCK data was negligible. 

E. ERROR ANALYSIS 

An error analysis is conducted to assign appropriate error bounds on the resultant 
spatially-dependent values. Approximated errors are the use of a linear wave theory 
transfer function to integrate the velocity covariance measurement over the water col- 
umn and errors in stationaritv due to tidal variation. Three sources of error which are 
specifically quantified are incident wave angle error, propagation of measurement error, 
and finite differencing error for the radiation stress gradient. 

1. Incident Wave Angle Error 

Radiation stress (equation 3) is sensitive to the incident wave angle (a). The 
determination of an accurate a is dependent on several other angles. The Zeiss elec- 
tronic surveying system was required to measure the sled yaw angle. The Zeiss error 
(Table 14) is assumed negligible in the a determination. The adjustments (Chapter V) 
to translate the sled instrument's three-dimensional coordinates into the FRF coordinate 
system, and the measurements by the portable Zeiss system which made them possible, 
are considered to have negligible error. Three angular sources of error which cannot be 
neglected are beach angle orientation error, current meter orientation error, and current 
meter directional gain error. 

a. Beach Angle Orientation Error 

An analysis is conducted to assess sensitivity to beach angle orientation er- 
ror. The beach angle orientation is calculated from a mean bathymetric contour deter- 
mined in a least squares sense as described in Chapter V. For this analysis, 
monochromatic, linear, unidirectional wind waves are assumed for a straight-and- 
parallel planar beach. Using Snell's Law, the wave incident angle (a 2 ) at the inshore 
current meter may be determined given the wave incident angle (oq) at the offshore cur- 
rent meter: 



*2 



sin 



-l 



sin a, C 2 



(98) 



where wave celerity (C) is determined from equation (70). The radiation stress gradient 
(from equation 3) in the cross-shore direction may be stated as a finite difference be- 
tween radiation stress calculated at the two on-offshore current meters where the sub- 
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script "rear" refers to the rear or offshore current meter and "front" to the front or 
inshore current meter of the sled: 




cx 



\_E n sin(a + Aa p ) cos(a + Aa p )J rcar — \_E n sin(a + Acy) cos(a + Aa e )3 

Ax 



f/Jfrom 



(99) 



where A is the beach angle orientation error. Total energy (E) is calculated from: 



Root-mean-square wave height (// rms ) for monochromatic shallow water waves can be 
related to depth (h) in the inner surf zone by equation (39). For SLPERDLCK, ratios 
of /7 rms to h ranged from 0(0.3) outside the surf zone to 0(0.4) inside the surf zone. 



resentative depth, slope, and incident wave angles is illustrated in Figure 25. The three 
angles plotted for each depth are centered about a typical incident wave angle encount- 
ered during SLPERDL'CK for that particular depth. For example, an incident wave 
angle of 25°. a depth of 3.5 m, and a beach orientation angle error of —5°, would result 
in an overestimation of the radiation stress gradient by 10%. On the other hand, a ra- 
diation stress gradient measured in a shallower depth with a reduced incident wave angle 
(h = 1.5 m, a = 15°) and the same error angle of —5° would result in an overestimation 
by 30%. Note that the overall governing parameter is incident wave angle. Therefore, 
for a fixed beach orientation angle error, the relative radiation stress error increases with 
decreasing wave angle, and hence with decreasing depth. For SLPERDL'CK, errors in 
mean beach orientation angle are estimated to be less than + 2° and 1.0 m was the 
minimum depth in which data were acquired. Therefore radiation stress gradient errors 
for 3.5. 2.5, and 1.5 m depths would be + 5%, ± 7-8%, and ± 10%, respectively. A 10% 
radiation stress error can be equated to approximately a 10% Cy error, as long as wind 
stress and temporal variability of mean momentum are small. As the magnitude of wind 
stress and temporal variability increase, the effect of a fixed radiation stress error on Cy 
becomes less. Varying beach slope from 0.01 to 0.10 or varying Ax from 2 m ("local" 
method) to 25 m ("integrated" method) had negligible effect on the relative radiation 
stress gradient error due to an incorrect beach orientation angle. Therefore Cy error due 
to inaccurate beach orientation angle is estimated to be less than 10%. 




(100) 



The radiation stress gradient (equation 99) versus orientation error for rep- 
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Figure 25. Sensitivity of cS iX / cx to beach angle orientation error for varying inci- 
dent wave angles and depths and tan /? = 0.03. 
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b. Current Meter Orientation Error 

Orientation of the current meters axes to the sled frame axes was considered 
to be within 1 — 2° (Chapter V). Assuming conservation of the radiation stress gradient 
outside the surf zone, all current meters were individually oriented by a previously de- 
scribed selection process (Table 6). This numerical rotation is considered to have elim- 
inated any current meter orientation error. 

To assess the effect of other rotations which minimized the radiation stress 
gradient and which had realistic spectra shapes, rotation subsets within Table 6 combi- 
nations were also examined. For example, 

• in run 16-2. rotating both current meters by a lesser amount instead of just the rear 
current meter by —2.86°, or 

• in run 15-1. rotating the front current meter by 0.27° instead of the rear current 
meter by — 0.29 0 , 

yielded no substantial difference in the magnitude of <y along that transect. In other 
words, by using the concept of conservation of S vx to initially orient the current meters, 
the resultant cj s were insensitive to the manner in which the orientation was done. 

c. Current Meter Directional Gain Error 

Current meter directional gain error is the error introduced by using a single 
gain and bias factor to calibrate current flows for all directions. The current meters were 
initially calibrated in a laboratory flow tank where gains and biases were determined for 
meter probe angles of 0°, 90°, 180°. and 270°. These angles correspond to the four 
sensor probes protruding from the current meter ball. Using the gain and bias factors 
for the u and v sensors determined for these four probe angles, the current meter was 
then checked against a controlled How of 0(0.5 ms -1 ) for meter probe angle increments 
of 15° . The standard deviations from the controlled flow for sled current meters 1, 2, 
and 4 were 0.012, 0.012, and 0.009 ms -1 , respectively. This error is included within the 
measurement error for current velocities addressed in the next section. 

2. Measurement Error 

To ascertain the effect of direct measurement error on the final determination 
of cy, an error propagation analysis is conducted for both the "local" and "integrated" 
methods for cj determination. Standard deviations (a) of directly and indirectly meas- 
ured variables used in equation (93) were estimated or determined from laboratory cali- 
bration results (Table 14). 

The coefficient of variation (Barn - . 197S), or fractional standard deviation 
(Young. 1962), is defined as: 
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Table 14. MEASUREMENT ACCURACY 



Direct Measurement 



parameter 


a 


Z X (%1 


tvpical value 


water density 


± 2 kam -3 


+ 0.2 


1022 kam- 3 


Zeiss vert control 1 


+ 0.03 m 


+ 1.5 


2 m 


Zeiss horiz control 1 


± 0.03 m 


+ 0.075 


40 m 


c/m velocity 


+ 0.05 ms -1 


± 5.0 


u = 0.25 ms' 


time 

wind stress 


insignificant error 
± 0.04 Jm~ 3 


+ 40.0 


v = 1. ms -1 
0.1 Jm- 3 


MWL (Paro) 


x 0.0003 m 


+ 0.015 


2 m 


wind speed 


+ 0.35 ms -1 


+ 5.0 


7 ms -1 


atmos. press. 


+ 2 mb 


+ 0.2 


1010.0 mb 


atmos. temp. 


±2°K 


+ 0.7 


2S8°K 


Indirect Measurement 


atmos. drag coeff. 


+0.5 x 10" 3 


± 30. 


1.5 x 10 -3 


atmos. density 


± 0.00S kgm -3 


+ 0.6 


1.25 kgm -3 



• where ± refers to an approximate standard deviation. 
1 from Birkemeier and Mason (19S4) 



v x =Y (ioi) 

where X is the arithmetic mean of replicate measurements of a variable. For the coef- 
ficient of variation determination (Table 14), the strict X has been replaced by typical 
(order of magnitude) values, which are listed. 

Calibration of the sled current meters was performed before and after the ex- 
periment by two different calibration laboratories. Pre-calibration gains and biases are 
listed in Table 15. Post-calibration gains varied by ± 3.0% and 0.01 - 0.04 ms -1 in bias 
or offset. A conservative estimate of measurement standard deviation for the current 
velocities is ± 5% to include the unmeasured environmental source of error between field 
and laboratory' conditions as well as the uncertainty (+ 3%) associated with using a 
single gain factor for all frequencies (Guza and Thornton, 1980) and directions. 
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The general rules of measurement error propagation can be stated for four 
arithematic operators as (Barr}'. 197S): 



2 , 2 , 2 , .2 
^sum or diff = ± \ °A + C B + °C + ••• + a n 



^ product - ± ABC ... n \/(~7~y + ( Y + + ••• + ~0~) 



g B \2 



7 C x 2 



V 



B 



C 



7 n s2 



(102) 

(103) 



^quotient — \J 



(■x) 2 + 'l» 2 



^povver i G £ 



(104) 

(105) 



where the standard deviations o A o B etc., refer to the + error of the parameters (e.g., 
A. B) being analyzed (Young, 1962). 



Table 15. SLED CURRENT METTR CALIBRATION 



CM 


Sled 


Sled 


PCM 








Correlation 


Ser. = 


Lctn list 


CM A 


Channel 


Vel 


Gain 


Bias 


CoefC (r 2 ) 


SI 14 


rear mid 


1 


6 


u 


0.951 


0.000 


0.99944 








5 


V 


0.916 


-0.00S 


0.99970 


S762 


front lower 


4 


13 


u 


1.156 


0.015 


0.99S99 








12 


V 


1.147 


-0.012 


0.99S7S 


S354 


side upper 


D 

L 


S 


u 


1.457 


-0.030 


0.99949 






7 


V 


1.457 


0.017 


0.99955 



• current meters were Marsh- McBirncy. Model 551 

• calibrations by U.S.G.S. Gulf Coast Hydroscience Center, NSTL, MS 

• formulation: velocity in ms -1 = (gain x voltage) + bias 



The wind stress measurement error has the largest V x of all the parameters in 
Table 14. V x can reach as high as 40%, depending on the manner and conditions under 
which the measurement was taken. Conditions during SUPERDUCK were conducive 
to accurate measurements of wind stress. These conditions were: 

• over-water measurements taken during onshore winds, 

• at the end of a 500 m pier, 
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at a height of 18.7 m MSL which is 11.1 m above the pier deck, 
using the eddy correlation method, and 

using a Campbell CA27T Sonic Anemometer as backup for the Gill anemometer. 












Standard deviation of atmospheric density (p a ) is estimated from equations 
(103), (104), and the equation of state for a perfect gas (Haltiner and Martin, 1957). 
Therefore, using a subjectively-determined and conservative value of ± 30% for the wind 
stress measurement error, a Q coefficient of variation of 30% is determined from 
equations (11) and (103). 

Rewriting equation (93) as : 



Cr = 



8Syx 

Sx 






P I w r I 



p\ u X 



dMy 



cr 



p l«/l’ 



(106) 



or 

c f = -A + B - C , (107) 

the measurement error in Cj. or the standard deviation of Cy, is the result of the additive 
o ' s of the three terms on the right-hand-side of equation (106), abbreviated as A, B. C 
in equation (107). Applying equations (102) to (104), results in standard deviations for 
the four terms in equation (107) as listed in Table 16. 

Standard deviation determination of the radiation stress gradient {o A ), requires 
o of (dSyJ / d.v), which in turn requires g of S y £. S } J or S vx are calculated from equation 
(78) using the directly measured parameters u', v',p, and h . An upper bound for o of 
S y [ can be formulated by approximating for shallow water waves: 



T cf~ 

S yx = ph C uv df = xphcov uv < pho u o v (108) 

where cov uv is the covariance of the u and v velocities and is bounded by (e.g. Davis, 
1986): 

cov uy < a u o v (1 09) 

Therefore o s t can conservatively be approximated as: 
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( 110 ) 



O T = 



± x /(^) 2 + <^) ! + (>) ! + < x ) ! 



o p and o h are negligible compared to o u and o v (Table 14) and may be neglected. 
One may then write: 

° = — \ S } )rear T { G S yj )front 



( 111 ) 



Sx 



+ 



A.y 




r + <^r» 



( 112 ) 



The standard deviation of Ax is much less than the standard deviation of A S y [ and 
therefore may be neglected. The standard deviation determined by equation (112) for 
the "integrated" method is reduced by a factor of one-half to account for that method's 
increased sample size (see equation S4). Standard deviation of the radiation stress gra- 
dient for each of the sled's data runs, using both the "local" and "integrated" methods 
for codetermination, were calculated. The mean a for radiation stress gradient of all four 
transects is listed in Table 16. 

The standard deviation for the wind stress (tJJ) was determined from equations 
(11) and (103), and Table 14. Standard deviation of 6M y / dt was determined by as- 
suming: 



6 M y _ A M v _ A(pDV) 
ci * A/ * At 



(113) 



since M y ' .\I y , and then applying equations (102) through (104) and Table 14 values. 

A typical mean Cj- value of 0.003 was used in equation (101) to determine the cj 
coefficient of variation in Table 16. Because of the relatively large incident wave angles 
and the measurement precision utilized in the SUPERDUCK experiment, propagation 
of measurement error resulted in reasonable cy values for both methods (Table 16). The 
"integrated" method's o C/ was superior to the "local" method because of the larger rela- 
tive error for the radiation stress gradient (equation 112) experienced when Ax is small. 
The integrated method's error bounds for the Cj range in Table 13 are + 4-15%. 



98 



Table 16. ERROR PROPAGATION ANALYSIS 



Term (eqn 107) 


a 


r*(%i 


c f Method 


A 


0.0015 


24 


local 




0.00013 


10 


integrated 


B 


0.00003 


31 




C 


0.00007 


70 




c f 


0.0015 


50 


local 


0.00015 


5 


integrated 


where term 








pc 

4 ‘'Ryjc 


T ldx 


t* 


cMJdi 


n 


/ 


r 


P\ 


u , 1 v 


p\u t \v 


p\u,\ V 



3. Finite Differencing Error 

Thornton and Guza (19S1) examined tire error involved in representing an ar- 
bitrary continuous function by discrete measurement points. Using a Taylor series ex- 
pansion for S' 



S yx {x + Ax) = S vx (x) + Ax 



dSyxi*) (Ax) 2 S 2 S yx (x) 



cx 



+ 






Ax 2 



+ 



(114) 



and then solving for the true radiation stress gradient: 
ASy X 



3Syx 

CX 



XI c 

A.v 0 



Ax 



21 



cx 



A.v 2 " 3 V 



- 3 
CX 



3! 



(115) 



or 



CX 



A Syx 
Ax 



Error 



(116) 



where the sum of the second and infinite terms on the right-hand side of equation (115) 
is the error of the finite difference space approximation. 

Relative error (R) may be defined as: 
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R = 



Error 



(117) 



6S 



VX 



cx 



Solving for R where only the first error term in equation (115) is considered, for 
monochromatic shallow water linear waves inside the surf zone may be expressed as: 

S yx = E sin a cos a a; AD SI ‘ (US) 

j- pgy 2 D 2 , 

) / + /;, and 

complicated constant 
Substituting equation ( 1 1 S) into (117): 

Rv 4 Ax IQjJL (i 19) 

4 h 

where D a h. (Thornton and Guza. 19S1) 

This approximation (equation 1 19) is adequate providing Ax and tan p are rea- 
sonably small. However, A.r and tan P are not always small in this application. There- 
fore R is redefined to include a sufficient number of terms on the right-hand side of 
equation (115) to allow convergence to within 5% of R. Resultant computations re- 
vealed that equation (119) was sufficiently accurate (within 5%) for most of the 
SUPERDL'CK data. However, using the "integrated" method between the most inshore 
data run and the beach required additional iterations (up to five) to stabilize R. 

The relative error in equation (117) is biased in one direction because the term 
of largest magnitude is written as: 



~r = slope (120) 

ch 

5x 



A-V ^ Syx 

21 



cx 



°°v.x 



cx 



where 



D 

A 
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and assuming negligible bottom curvature. Therefore the sign of R is fixed as the same 
sign as beach slope. The finite approximation for the radiation stress gradient may be 
written as: 



ASy X 

Ax 



(1 + /?) 



*S yx 

cx 



(121) 



where R represents the overestimation of dS yx / dx for positive slopes and the underes- 
timation for negative slopes. Finite difference relative error for the four transects of this 
research are listed in Table 17 and Table 18. They are listed for both methods of co- 
determination. The larger error for the "integrated" method value between the most in- 
shore data run and the beach was due to the combined effect of increased Ax, decreased 
depth, and steeper slope, all which combine to increase R. The relative error is signif- 
icantly reduced when using the "local" method due to its much smaller Ax. "Integrated" 
method's R ranges from -27% on the nearshore side of the bar to nearly 200% on the 
nearshore slope. However, except for the gross error on the nearshore slope, "inte- 
grated" method's relative error is less than 50%. This would cause Cy error due to finite 
differencing to be less than 50% (assuming negligible wind stress and steady state con- 
ditions in equation 93). 

4. Combined Error 

The various sources of error are now combined into a total, or combined, error 
relative to the bar. To establish a conservative combined error, finite differencing error 
is assigned a range of error rather than a bias and unquantified errors are estimated at 
± 20%. Mean finite difference error for the trough does not include the gross error. 
0(200%), between the most inshore sled run position and the beach. Therefore the final 
cj- value for the trough does not reflect the nearshore slope. 

The maximum combined error is determined by adding the four maximum pos- 
sible errors listed in Table 19. However it is not reasonable to define maximum error 
in this manner since each of the four errors have a range of uncertainty associated with 
them and they are uncorrelated. It is very unlikely that the errors would all be a maxi- 
mum and in the same direction at the same time (Barry, 1978). Hence, assuming the 
errors are independent of one another, the combined error is obtained from equation 102 
and results in final momentum balance cj values of 0.004 ± 0.0010, for offshore the bar, 
0.002 ± 0.0006 for on top and immediately before the bar, and 0.001 + 0.0003 for the 
trough region. 
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Table 17. FINITE DIFFERENCE RELATIVE ERRORS 



October 15 1986 Relative Error (R) 






Sled Run 


A.v (m) 


tan /? 


MWL (m) 


R (%) 


15-1 


1.99 


0.0199 


2.80 


l. 




27.99 


0.0238 


2.29 


23. 


15-2 


1.99 


0.0287 


1.78 


2. 




24.24 


0.0241 


1.72 


27. 


15-3 


1.99 


-0.0031 


1.67 


0 . 




44.65 


-0.0127 


1.99 


-20. 


15-4 


1.99 


0.0554 


2.30 


4. 




42.03 


0.0549 


1.15 


197. 


Beach 










October 16 1986 Relative Error (R) 






16-2 


1.99 


0.0158 


3.33 


1. 




50.52 


0.0290 


2.61 


46. 


16-3 


1.99 


0.0349 


1.89 


3. 




17.69 


0.0341 


1.70 


28. 


16-4 


1.99 


0.0309 


1.51 


3. 




12.99 


0.0198 


1.52 


13. 


16-5 


1.99 


0.0145 


1.54 


1. 




33.49 


-0.0081 


1.72 


-12. 


16-6 


1.99 


0.0029 


1.91 


0 . 




33.88 


0.0564 


0.96 


196. 


Beach 
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Table 18. FINITE DIFFERENCE RELATIVE ERRORS 



October 17 1986 Relative Error (R) 



Sled Run 


A.v(m) 


tan /? 


MWL (m) 


R (% 


17-1 


1.99 


0.0182 


3.S5 


1. 




34.60 


0.0224 


3.22 


19. 


17-2 


1.99 


0.0397 


2.59 


2. 




26.80 


0.0364 


2.01 


* 39. 


17-3 


1.99 


0.0309 


1.43 


3. 




9.54 


0.0265 


1.31 


15. 


17-4 


1.99 


0.0221 


1.19 


3. 




4.09 


0.0232 


1.14 


6. 


17-5 


1.99 


0.0166 


1.09 


2. 




31.56 


-0.0153 


1.25 


-27. 


17-6 


1.99 


-0.0036 


1.40 


0. 




4.88 


0.0109 


1.48 


3. 


17-7 


1.99 


0.0137 


1.55 


1. 




6.27 


0.0276 


1.49 


9. 


17-8 


1.99 


0.0495 


1.43 


5. 




20.06 


0.0716 


0.72 


197. 


Beach 










October 18 1986 Relative Error (R) 






18-1 


1.99 


0.0113 


3.80 


0. 




48.33 


0.0259 


3.09 


32. 


18-2 


1.99 


0.0405 


2.38 






18.85 


0.0369 


1.94 


29. 


18-3 


1.99 


0.0334 


1.49 


j. 




21.60 


0.0258 


1.31 


34. 


18-4 


1.99 


0.0083 


1.13 


1. 




27.67 


-0.0096 


1.30 


-15. 


18-5 


1.99 


-0.0223 


1.48 


-2. 




33.71 


0.0438 


0.74 


196. 



Beach 



103 



Table 19. COMBINED ERROR RELATIVE TO THE BAR ON THE BED 
SHEAR STRESS COEFFICIENT AS DETERMINED FROM A MO- 
MENTUM BALANCE 





Trough 


On the Bar 


Offshore 


beach angle orientation error 


± 8% 


± 10% 


± 10% 


measurement error 


± 5° o 


i 5% 


i 5% 


finite differencing error 








mean R 


-19% 


19% 


31% 


^ error = 1/(1 + R) 


23% 


-16% 


-24% 


unquantified errors (est.) 


± 20% 


± 20°. o 


± 20° o 


combined <y error 


± 32% 


± 28° o 


+ 33% 


final <y 


0.001 ± 0.0003 


0.002 ± 0.0006 


0.004 ± 0.0010 
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VIII. LONGSHORE CURRENT MODEL 



A. INTRODUCTION 

A first-order longshore current model for a barred beach is proposed. The model is 
based on balancing the wave-induced radiation stress gradient with longshore bed shear 
stress, turbulent momentum exchange, and longshore surface wind stress. Stationary 
wave conditions and straight and parallel bottom contours are assumed. The waves are 
described as random, but narrow-banded in frequency and direction. Model results are 
compared to field observations and bed shear stress coefficients are determined by 
model-fitting. 

B. MODEL THEORY 

There are no field- verified longshore current models for a barred beach. Therefore 
a model which has been field-verified for a planar beach was selected for this analysis. 
The model (Thornton and Guza, 19S6; hereafter referred to as TG86) is based on a 
simple balance between the gradient of the wave-induced radiation stress, the bed shear 
stress, and the gradient of turbulence-induced radiation stress. The model's assumptions 
are stated in the previous paragraph. The longshore current model is based on a random 
wave height transformation model (Thornton and Guza, 1983; hereafter referred to as 
TG83). Modifications to the TGS6 model are made based on the availability of addi- 
tional data from SUPERDUCK and with the purpose of improving and broadening the 
model's applicability to barred beaches. These modifications are the inclusion of surface 
wind stress in the momentum balance, retention of incident wave angle dependency, and 
the use of a breaking wave height distribution based on a slope-dependent y and an im- 
proved empirical expression for the breaking wave percentage as a function of II, h, and 
y. Whereas TG86 found turbulent momentum exchange to be unimportant for their 
planar beach topography, turbulent momentum exchange is shown to be required in 
order to model current velocities in the trough region of a barred beach. The TG86 
longshore current model is briefly summarized below with modifications described later. 

1. Wave Height Transformation 

This portion of the model is an extension of work by Battjes and Jannsen (1978). 
The wave height transformation is conceptualized as an energy flux balance given by: 
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<h> 



c(EC s COS 7.) 
c.x 



( 122 ) 



where C, cos si is the shoreward component of wave group velocity, a is the mean inci- 
dent wave angle of the narrow-banded waves, and < e b > is the ensemble-averaged dis- 
sipation due to wave breaking. Analysis of SUPERDL'CK data, as well as results from 
TGS3. indicated that the ensemble-averaged dissipation due to bottom friction < e,-> is 
negligible compared with < i b >, and therefore < is omitted in equation ( 122 ). 

Based on analysis of field data from Torrey Pines Beach, CA, the Rayleigh wave 
height probability density function (pdf). /?(//), was shown to qualitatively describe the 
measured wave heights across that planar beach. The Rayleigh pdf is described by: 
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where II is wave height. The Rayleigh distribution describes all waves. 

Only the breaking waves are dissipative, resulting in changes in the momentum 
flux, i.e. wave forcing. Therefore a distribution describing the breaking waves alone is 
required. Wave height probability density functions were calculated from field meas- 
urements taken at Ft. Ord's Soldiers Beach. Monterey. CA. These pdfs were used to 
more realistically model random wave height transformation across the surf zone. 

Within the surf zone, breaking wave height distribution is approxi- 

mated by weighting the Rayleigh pdf: 

P b m = W(lh P [H) (124) 



where the weighting function Will) < 1. It is noted the breaking wave distribution is 
not a pdf. but that the area of p b {H) gives the percentage of the waves that are breaking. 
A simple form of W(H) which approximates the fraction of waves that break from off- 
shore to saturation conditions is: 



W(H) 



Ijrms 

yh 



4 



(125) 



where 




( 126 ) 
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for saturated wave breaking conditions. Based on the Ft. Ord data, a modification to 
equation (125) which more heavily weights the larger waves is given by: 



\V(H) = 



H r ms 
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(127) 



Wave breaking dissipation is modeled as a simple periodic bore: 

T (bh? 

£ » = 4 ** ~1T 



(12S) 



where B is a breaker coefficient of 0(1) and / is the mean wave frequency which is the 
peak frequency for narrow-band wave spectra. The B coefficient accounts for differences 
in various breaker types and is a function of the proportion of foam region on the 
breaker face (Figure 26). The peak frequency in the SUPERDL’CK S yx spectra is used 
for/. 




The ensemble-averaged dissipation due to wave breaking is calculated by 
multiplying the dissipation for a single broken wave of height H by the probability of 
wave breaking at each height and integrating over all H: 
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By substituting from equations (124) and (127) and evaluating the definite integral, 
equation (129) may be written as: 




(l+(H r m S h/h) 2 f 2 



(130) 



Energy flux [EC g from equation 122), is approximated using linear theory for C g at fre- 
quency/. and energy is calculated based on the Rayleigh pdf: 



Substituting equations ( 1 30 ) and (131) into (122) yields a first-order differential equation 
for the linear wave height transformation model: 



A spatially-dependent wave height [//, m |.r)] can then be determined for an arbitrary 
bottom profile across the surf zone by numerical integration, given an initial value of 
II rms . The H rms profile is used as input for the longshore current calculation. 

2. Longshore Current 

The longshore current distribution across the surf zone is based on conservation 
of momentum. Rewriting equation (3). wave-induced radiation stress may be written: 
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S, Y = EC 0 cos y. 
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by Snell's Law of linear wave refraction. The subscript "o" refers to any arbitrary lo- 
cation. Therefore the wave-induced radiation stress gradient may be written as 



4— [EC a cos a) = 
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where < > is substituted from equation (122). Assuming steady state conditions and 

negligible wind stress, the longshore momentum balance is written as: 

. A 



sin a. 0 



<c b > 




(136) 



TGS6 used equation (136) and model-fitting to solve for the horizontal profile 
of longshore current across a planar beach. They obtained a hierarchy of solutions by- 
choosing combinations of linear or non-linear bed shear stress, with and without turbu- 
lent momentum exchange. Turbulent momentum exchange was found not to be im- 
portant for their planar beach topography. They concluded that a non-linear bed shear 
stress without turbulent momentum exchange was their preferred choice for general ap- 
plications, based on its fit to their data. A mean c^-for their four days of data was 0.006 
± 0.0007. 

C. MODIFICATION OF MODEL BASED ON FIELD DATA 

Modifications to the TGS6 longshore current model are the inclusion of surface 
wind stress, incident wave angle dependency, and use of a breaking wave weighting 
function based on a slope-dependent y and an improved empirical expression for break- 
ing wave percentage as a function of H, h, and y. The use of turbulent momentum ex- 
change is found to be important for modeling longshore current velocities in the trough. 
The inclusion of surface wind stress and incident wave angle dependency are covered in 
the Model Results section. The other modifications are addressed next. 

1. A Slope Dependent y 

Since the TG83 model was developed for planar beaches, y was assigned a con- 
stant value of 0.42 (equation 39) and is the ratio H rm Jh for saturated wave breaking 
conditions. However. Sallenger and Holman (1985) found y to be a function of beach 
slope based on data from several field experiments. They determined an empirical re- 
lation for a slope-dependent y based on their data: 

y = 3.2 tan/? + 0.30 (137) 



The ratio H rms jh is compared to beach slope (Figure 27) for both the SL'PER- 
DUCK data and the data used by Sallenger and Holman (1985). The solid line in 
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Figure 27 is Sallenger and Holman's (1985) empirical relation for y. Most of the 
SUPERDUCK data has the same magnitude of scatter as previous data sets about the 
solid line, with even negative slopes (found on the inshore side of the bar) in general 
agreement. The exception is for SUPERDUCK data sets located in the trough and on 
the foreshore slope (Runs 11-3, 15-4, 17-6, 17-7, 17-8, 21-4). They have a significantly 
lower value of H rm Jh and are indicative of non-saturated breaking wave conditions. 
Obviously a slope-dependent y relationship is required for a barred beach. Based on 
Figure 27 and the general fit of the SUPERDUCK data, equation (137) is used for the 
SUPERDUCK analysis. 

2. Breaking Wave Weighting Function 

The weighting function, \V(H), used in TG83 was based on only four data runs 
for a single transect of a beach with a tidal plateau. To improve on the formulation of 
\V(H). breaking wave pdfs are calculated for the SUPERDUCK data. Time series of 
pressure head. are linearly transformed into >/(/) time series by first calculating 

the complex Fourier spectra of />/,(/), then multiplying by a linear transform function 
[11(0], band pass filtering between 0.03 - 1.00 Hz, and finally applying an inverse 
Fourier transform to the ?/(f) function. 

The transfer function 11(0 (from equation 38) is held constant at the 0.5 Hz 
value for the interval 0.5 Hz to 1.0 Hz so as to include the additional higher frequency 
variance for the H rms calculations, but not to overly bias them with the exponentially 
increasing 11(0 values. Filtering is accomplished by zeroing out the Fourier amplitude 
coefficients of the undesired frequencies before applying the inverse Fourier transform. 
The entire 34. 1 min record is transformed at one time. As an example, the initial three 
minutes of Run 15-2's p h {t) time series is compared with its transformed ?/(/) series in 
Figure 28. The transfer function H(0 and the G^(/) spectra determined from the entire 
;/(/) series are also included. Wave periods and heights are determined from the ?/(r) time 
series by the down-crossing method. Root-mean-square wave height is determined sta- 
tistically from the calculated ?/(/) time series as: 




(138) 



The statistical H rms from the >/(r) series are compared with the energy-based H rms 
(equation 90) determined from both the G^if) spectra discussed in Chapter VII and the 
G n (f) spectra calculated from the t](i) record (Figure 28). Root-mean-square wave 
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Figure 27. 



Ratio H rm Jh versus beach slope for various data sets, 
represent trough and foreshore slope locations. 



Circled data 
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height from all three sources compared favorably [within (10%)], with the difference 
indicating a measure of non-linearity in the surface elevation record. 

As discussed in Chapter IV, visual observations of wave dissipation (i.c., wave 
breaking) at the sled were electronically marked on the data set records. These marks 
are then superimposed on the calculated ?/(/) time series to identify which waves were 
breaking during each 34.1 min data run. From this procedure, wave heights and periods 
for all waves, and a subset for just breaking waves, are obtained. Pdfs of all waves and 
breaking waves for 15-17 October are illustrated in Figure 29 and Appendix A 
(Figure 47 through Figure 49). Since the data were measured sequentially and not si- 
multaneously, wave heights are normalized by the mean offshore wave height to account 
for any variability in the incident wave conditions. The pdfs are normalized to the mean 
offshore wave height measured by a WAVERIDER buoy located 6 km offshore. 
Equipment malfunction precluded breaking wave identification on other days of the ex- 
periment. 

Breaking wave data from SLPERDLCK, as well as the four data runs from Ft. 
Ord (TGS3), are compared in Figure 30. The ordinate is the fraction of waves which 
are breaking (Q) and the abcissa is the ratio, H rtr Jyh , where y is slope-dependent and 
calculated from equation (137). Thus for saturated conditions, both the ordinate and 
abcissa are equal to unity. The breaking wave fraction is obtained by dividing the 
number of waves identified as breaking by the total number of waves. The total number 
of waves in a data run was approximately 400. The Ft. Ord data are indicated by an , 
SL'PERDL'CK data acquired in the trough are solid circles, and all other SL'PERDL'CK 
data are indicated by open cirlcles. Three breaking wave weighting functions are illus- 
trated. The dotted line is the TGS3 weighting function (equation 127). The chain-dot 
line is from a theoretical expression (Battjes and Janssen, 197S; Battjes and Stive. 19S5) 
where the cumulative probability distribution of all wave heights (breaking or non- 
breaking) is the Rayleigh distribution, but abruptly terminated at the breaking wave 
height: 



Q = 



^_rms_ 

yh 



In 2 + 1 



(139) 



Neither of these two functions provide a reasonable fit to the data. However, a function 
which substantially improves the fit to the data is: 
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the solid line and empirical breaking wave distributions given by 



equations (127) and (140) are the dotted and dashed lines, respectively. 
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O = W{H) = 



(140) 
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and is delineated by a dashed line. This function (equation 140) has no physical justi- 
fication but is merely an empirical expression derived to better fit the non-trough data. 
It is noted that none of the three functions fit both the trough and non-trough data. 
The right half of the right-hand side of equation (140) is taken from equation (127) and 
is used to more heavily weight the larger waves towards breaking. Observations show 
that the largest waves at a particular location are more likely to break (TG83). This 
skewing of the pdf is substantiated by the breaking wave observations (hatched areas) 
seen in the pdf figures. Very little wave breaking occurs until H rms jyh reaches 0.8, then 
the percentage of breaking waves increases rapidly. The addition of the SUPERDUCK 
data suggests a more rapid increase in the breaking wave percentage within a smaller 
H rm Jyh band — almost to the point of describing a "step" function indicative of all waves 
breaking at a fixed breaker point. 

The ratio H rm Jyh in Figure 30 has a certain amount of error associated with it. 
The error in y is the result of an imperfect fit to the data in Figure 27. Tidal fluctuations 
also affect the h term. Thus there are runs which have H rm Jyh ratios greater than one, 
yet were unsaturated. 

It is also noted that during field identification of the breaking waves, breaking 
wave percentages were visually estimated to be often in the range between 60-90%. Yet 
subsequent data analysis reveals the percentages to be much lower. This may perhaps 
be similar to the phenomena associated with visually estimating wave heights from ships, 
where it has long been recognized that the human eye perceives an average of the highest 
13 wave heights and not the true mean wave height. This lead to the designation of ship 
wave height observations as "significant wave heights" and not mean wave heights. In 
a similar manner, SUPERDUCK visual observations of breaking wave percentages were 
overestimates. Thus visual observations of saturated wave conditions from other ex- 
periments may need to be reviewed. 

D. MODEL RESULTS 

1. Wave Height Transformation 

As input to the wave transformation model, H rms and / measured during the 
orientation runs outside the surf zone (Runs 15-1, 16-2, 17-1, 1S-1) are used as initial 
conditions. Mean wave frequency is assigned as the frequency of the spectral peak in 
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Figure 30. Breaking wave fraction (Q) versus the ratio H rm Jyh. 
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sea surface elevation. Mean incident wave angle is then assigned as the angle for/ from 
the significant wave spectra for that day. TG83 and TGS6 assumed small wave incident 
angles and set the cos a term in equation (122) to unity. Because of the large incident 
wave angles during SUPERDUCK, incident wave angle dependency is retained 
throughout the analysis by inclusion of the cos a term. A forward stepping scheme is 
used to evaluate equation (122): 

(EC g cos a) 2 = (ECgCOsa)] + (<£ 6 >)]A.v (1-11) 



All variables in equation (141) are known except B, which is determined by model-fitting 
the calculated H rms to the field-measured H rms . Model-fitting is accomplished by mini- 
mizing the least square error between observed and modeled data. Results of the 
model-fitting are shown in Figure 31 and Appendix A (Figure 50 through Figure 52). 
Field H rms is represented by an "x" for each of the sled runs and is determined from 
pressure sensor measurements. The spatial step size, A.v, was varied from 0.5 m to 10.0 
m with negligible effect on the H rms profile. Therefore A.v = 5.0 m was arbitrarily se- 
lected to allow suitable cross-shore definition along with reduced computational time. 

A constant B coefficient of 0(1) is determined by a least square fit for all four 
days (Table 20). Field and modeled H rms show good agreement except for the trough 
area on 15 and 16 October. Root-mean-square wave height is undermodeled by about 
35% for these exceptions. Spatially-variable B coefficients were also determined which 
gave exact fits to the field data. However their drastic spatial changes in B appeared 
unrealistic since the breaker type was observed to be similar for most of the transect, 
except on the very nearshore slope inshore of the last sled runs. A constant B coefficient 
for a transect was deemed appropriate. 

2. Longshore Current Model Results 

Longshore current is modeled by including surface wind stress in the momentum 
balance equation of equation (136): 
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(142) 



and using various combinations of the following approaches: 

• linear or non-linear formulation of xj!, 

• with and without turbulent mixing, 

• spatially-variable or constant Cy and X. 
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Figure 31. Model-generated H rms and V compared m ith field observations versus 
distance for invariant Cj- and B and including turbulent momentum ex- 
change, 16 October. Model results are solid lines and field observa- 
tions are "x'"s. Elevation is referenced to MSL. 
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In all cases, the spatiallv-invariant B coefficients listed in Table 20 are used. As with the 
B coefficient determination, Cj and X are model-fitted by a least square error criteria. 

A spatiallv-invariant cj- with a non-linear (and presumably more realistic) t* and 
without turbulent mixing results in poor agreement between modeled and field-measured 
longshore current velocities (Figure 32 and Appendix A, Figure 53). Although the 
general shape of the observed longshore current profile is modeled, there are large dis- 
crepancies between observed and modeled current velocities, especially in the trough re- 
gions. Current velocities in the trough region are in error by as much as a factor of 7. 
Order of magnitude for Cj- is similar to the momentum balance determination of cj- 
(Chapter VI 1) for all 4 transects. 

Allowing Cj to vary spatially provides some improvement in the model-fit. 
However major discrepancies still occur on the nearshore slope where a steep increase 
in velocity is modeled along with significant underestimation of velocity in the trough. 
The spatially-variant cr solution is dismissed because of these discrepancies. 

Although TGS6 determined turbulent momentum exchange to be unimportant 
for their planar beach data, turbulent momentum exchange is investigated for this anal- 
ysis as a possible mechanism to transfer momentum across the bar and into the more 
placid trough region where wave breaking is reduced. Turbulent momentum exchange 
is expressed as the last term in equation (142) and is approximated using an eddy 
viscosity formulation: 



S' " = 
°>\ x 



-h 



pu'v'dz = - pcD 
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where D is total depth and the eddy viscosity coefficient (c) is parameterized by a 
dimensionless and adjustable coefficient (X) from equation (21). A first approach for 
the inclusion of turbulent momentum exchange was to determine the minimum least 
square error for spatiallv-invariant cj and X across a reasonable range of cj and X values. 
Bed shear stress coefficients were varied from 0.0005 to 0.0200 and X was varied from 
0.000 to 0.016. A three-dimensional depiction of their solution surface is illustrated in 
Figure 33 for a typical data run. Minimum least square error is represented by a peak 
in the solution surface closest to the top of the three-dimensional box. Xote that there 
are many combinations of cy and X that minimize the least square error to nearly the 
same degree. Thus too many degrees of freedom result in multiple mathematical sol- 
utions with a corresponding loss of insight for the physics involved. 
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Figure 32. Model-generated V compared to field observations versus offshore dis- 
tance for constant no turbulent momentum exchange, 15-16 

October. Model results are solid line and field observations are "x"'s. 
Error is mean least square error. 



120 




almost constant ridge for optimal Cj- and N indicating a non-unique 
solution. 
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The spatiall} -invariant Cy approach is therefore considered a baseline upon 
which to use turbulent momentum exchange to improve the model results. The proce- 
dures are to first determine an invariant Cy, and then to vary N until the least square error 
is minimized. This results in a substantial improvement in model results over previous 
approaches. The IS October data run required an additional Cy iteration after a value 
of N was determined. The model result prior to this additional iteration is represented 
by a dashed line in Figure 52. With longshore momentum being transferred over the 
bar and into the trough, the previous velocity underestimation in the trough has been 
significantly reduced. Velocities are modeled within 20% of measured values, with the 
exception of the trough region on 17 October. This case is unique and is addressed be- 
low. Interval bars in Figure 31 and Appendix A (Figure 50 through Figure 52) repre- 
sent 10% and 20% of the observed value and are provided solely as a visual aid for 
interpreting the data. Mean least square error ranges between 7-15% for all four 
transects. 

Current velocities in the trough (Runs 17-6, 17-7, and 17-S) on 17 October are 
much lower (and actually of opposite direction to the longshore current offshore the bar) 
than on the other three days. Guza et al. (19S6) showed that it is S y J, and not the 
structure of S vx (f), that is important for longshore current response. It is therefore pos- 
sible that the low current velocities in the trough on 17 October are the result of bi- 
directional wave forcing where breaking swell and sea are driving longshore currents in 
opposing directions and cancelling each other out. Radiation stress spectra (Appendix 
F) show weak wave forcing for the trough region on 17 October. Radiation stress 
spectra for runs 17-1 and 17-2 show sea and swell approaching from different directions 
with relatively large incident angles (1S° and —9°, respectively). The cross-over wave 
period separating swell from sea is 6.7 seconds. A scatter plot of breaking wave heights 
versus breaking wave periods for a trough run (17-6) shows that approximately half of 
the breaking waves consisted of swell and half consisted of sea (Figure 34). A scatter 
plot for Run 17-4 (located offshore the bar) does not show as much breaking wave ac- 
tivity for the swell. Thus the minimum combined velocity in the trough on 17 October 
and the poor model results for that particular day and location are possibly a result of 
opposing longshore currents. 

3. Bed Shear Stress Coefficients 

Spatially-invariant Cy's are determined from model-fitting using a non-linear 
formulation of t* (TGS6) and turbulent momentum exchange (Table 20). The ratio of 
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Figure 34. Scatter diagram of breaking wave periods and heights for 17 

October. Run 17-4 is offshore the bar and 17-6 is in the trough. 
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Table 20. MODEL-FITTING RESULTS FOR B, cy AND N 





15 Oct 


16 Oct 


17 Oct 


18 Oct 


B 


0.74 


0.78 


1.10 


1.06 


N 


0.006 


0.015 


0.004 


0.016 


Error 


0.146 


0.071 


0.078 


0.114 


c /' 


0.003 


0.003 


0.009 


0.003 




0.003 


0.002 


0.004 


0.002 


Error is mean least square error for the transect 







'transect mean Cf deterntined by model-fitting 



2 

transect mean cy determined by momentum balance 



wind force to wave force for the SUPERDL'CK data was 0.09 (see Chapter VII), there- 
fore the inclusion of surface wind stress in this model did not have significant impact on 
the magnitude of cy from model-fitting for this particular data set. However, as discussed 
earlier, the inclusion of surface wind stress under certain wind, wave, and beach slope 
conditions can have significant impact. Bed shear stress coefficient values for 15-18 
October are 0.003, 0.003, 0.009, and 0.003, respectively. N coefficient values are 0.006, 
0.015, 0.004, and 0.016, respectively. These Cy values are in good agreement (see 
Table 20) with the mean <y s determined for each day by the momentum balance 
method. This agreement is encouraging and lends credibility to the Cj results considering 
the complex physics being addressed and the different and independent methods em- 
ployed. 
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IX. SUMMARY AND CONCLUSIONS 



Previous investigators of longshore currents have used many simplifying assump- 
tions which often preclude meaningful and realistic results. The assumptions of a planar 
beach, a steady and depth uniform flow, invariant bed shear stress, and negligible wind 
stress are quantitatively examined by measuring the relative importance of each term in 
a local longshore momentum balance based on measurements acquired during a field 
experiment at a barred beach. Wind and wave forcing of longshore currents is investi- 
gated with a numerical model and from field measurements. Measurements of pressure, 
current velocities at three levels, wind speed, and sea surface slope were acquired from 
an instrumented sled positioned at various locations across a transect of a barred beach. 
Unfortunately the slope measurements were found during data analysis to be unreliable 
and were disregarded. 

Three-dimensional depictions of the current structure over a nearshore bar are de- 
veloped. The cross-shore velocities C0(O.2ms _, )3 indicate a relative onshore flow in the 
upper third of the water column and a relative offshore flow in the lower third of the 
water column, irrespective of bar location and in agreement with cross-shore theory. 
The longshore velocities C<?(1.0ms -1 )] indicate nearly depth uniform flow with a slight 
velocity increase with elevation. This depth uniform flow confirms the appropriateness 
of the depth-integration aspect of radiation stress for a barred beach. The longshore 
vertical current structure did not exhibit a logarithmic profile. 

A numerical model is developed to compare wind and wave forcing for varying wind 
conditions and beach slope. The finite-depth TMA spectral wave model is used as input 
for the energy spectrum. Waves are assumed to be locally generated and fully arisen. 
Wind force is proportional to surf zone width and therefore inversely proportional to 
beach slope. Wave force is independent of surf zone width and beach slope. Wind force 
is found to increase more rapidly than wave force for an increase in wind speed or an 
increase in the wind and wave incident angle relative to beach normal. Ratios of wind 
to wave force are calculated. For a given wind direction, a wind speed increase of 10 ms 
results in a 50% increase in this ratio. In other words, wind force becomes increasing 
significant relative to wave force, as wind speed and incident wind direction increase. 
For obliquely onshore winds and a steep beach slope (tan /? = 0.10), wind force is gen- 
erally less than 10% of the wave force and thus may be considered negligible. For ob- 
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liquely onshore winds and a gentle beach slope (tan /? = 0.01). wind force can be 
significant, with wind to wave force ratios ranging from 15% to 100% for winds between 
10-30 ms -1 and wind and wave directions from 10 — 70°. Thus under certain wind and 
wave conditions, wind force can be a first order term in the longshore momentum bal- 
ance. For offshore winds, longshore currents are totally wind-driven. In summary, the 
model shows that wind force should be considered as a forcing term for longshore cur- 
rents and that it is inappropriate to assume that wind force can be omitted. 

There are no field studies addressing surface wind stress for robust surf zones ( H mo 
greater than 0.5 m). Surf zone atmospheric drag coefficients are determined from surface 
wind stress measurements acquired immediately outside the surf zone and wind speed 
measurements acquired by the sled anemometer inside the surf zone. Stability- 
dependent drag coefficients are calculated using the Businger et al. (1971) iterative ap- 
proach. A conservative increase of 33% for Q was made to account for the significantly 
increased surface roughness of the surf zone due to breaking waves based on the results 
of Davidson et al. (19SS) which indicate an increased Q for chaotic seas. During periods 
of atmospheric convective instability and reduced surface roughness (low wave heights), 
Hsu's (1970) empirical relationship for C d is used. Resultant C d 's range from 1.5 x 10 -3 
to 3.4 x 10 -3 for non-convective instability and from 4.6 x 10~ 3 to 8.25 x 10 -3 for 
convective instability. These C d magnitudes are consistent with previous field measure- 
ments of C d over the open ocean as well as the theoretical C d values of Geernaert et al. 
(1987) for shallow water depths. 

All terms in the longshore momentum balance are calculated from field measure- 
ments of pressure, current velocity, and wind speed. The sled was initially positioned 
beyond the surf zone for a complete data run. This additional data run allowed a daily 
numerical calibration of the current meters during data analysis based on conservation 
of SSyJdx outside the surf zone. This procedure is essential for obtaining realistic 
physical results and is a distinguishing feature of the research. Longshore momentum 
due to unsteady flow was determined to be less than 10% of the longshore momentum 
due to steady flow for this data set. Relative contributions of the radiation stress gra- 
dient (cS vx / c.v), temporal term ( dMJdt ), and surface wind stress (tjj) to the longshore 
momentum balance were calculated and analyzed by location relative to the bar. Mean 
relative contributions for the wave, wind, and temporal terms were 0.86, 0.08, and 0.06, 
respectively, for both offshore and on top the bar. However in the trough, the mean 
relative contributions were significantly different (0.64, 0.17, and 0.19. respectively). For 
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the trough region, the relative importance of wave forcing is reduced by 25%, whereas 
the relative importance of surface wind stress is doubled and the temporal term relative 
importance is tripled. All three terms are of first order in the trough. Although surface 
wind stress was nearly constant across the surf zone, the relative contribution of surface 
wind stress increases shoreward from the location of maximum breaking due to the de- 
crease in wave forcing. The relative contribution of the temporal term appears to be 
largest in the trough, possibly due to the presence of surf beat or edge waves. Wind to 
wave force ratios calculated from the field data (0.09) are in agreement with those de- 
termined from the model (0.12) for the same wind conditions. 

Spatially-depcndent bed shear stress coefficients are determined between sled run 
locations. The spatial step size (Ax of the radiation stress gradient) is determined by 
error analysis to be too small when assigned the 2 m distance between the current meters 
on the sled. It had originally been planned that the slope array, which did not work 
properly, would have provided the necessary precision to measure the radiation stress 
gradient over the 2 m distance. Reasonable Cj results are calculated when the spatial step 
size is increased to the distance between run positions 0(20 - 40 m). Resultant cj s range 
from less than 0.001 to 0.006. Spatially-variable mean cj s are found with values of 0.004 
+ 0.0010 for offshore the bar, 0.002 + 0.0006 for on top the bar, and 0.001 + 0.0003 for 
the trough (excluding the nearshore slope). A trend for a shoreward decreasing cy is 
noted, possibly due to the shoreward decrease in wave action. A mean <y for the 
nearshore slope is not specified since error analysis shows that a significant overcsti- 
mation of ty can occur for the nearshore slope due to the combined effect of decreased 
depth and steeper slope upon the finite difference approximation. The percentage 
change in <y due to neglecting surface wind stress of the temporal term in the momentum 
balance ranges from near zero to 100%, with a mean change between 10-30%. Thus for 
this data set and under certain conditions, the surface wind stress term and the temporal 
term can be of first order in the momentum balance. 

Longshore current velocities at a barred beach are modeled to within ± 20% of the 
observed values at SUPERDUCK by modifying the Thornton and Guza (19S6) model 
to include wind stress, wave incident angle dependency, slope-dependent y, and an im- 
proved breaking wave weighting function. Mean least square error between modeled 
and observed current velocity ranges from 7-15%. The inclusion of turbulent momen- 
tum exchange to transport momentum into the more placid trough region is required 
and is physically reasonable. Spatially-invariant ty values determined by model-fitting 
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for 15-18 October are 0.003, 0.003, 0.009, and 0.003, respectively. These cy values agree 
well with mean Cy values determined for each day by the momentum balance method. 
This agreement is encouraging and lends credibility to the cy results considering the 
complex physics being addressed and the different and independent methods employed. 
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APPENDIX A. 



DATA ANALYSIS FIGURES FOR DAYS OTHER THAN 



OCTOBER 16, 1986 
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Figure 35. Mean current structure on 1 1 October 1986. Numbers represent data 
runs, short dotted lines are MWL's, and distance alongshore is arbi- 
trary. 



130 



SUPERDUCK 
12 Oct 86 
Profile 250 




Figure 36, Mean current structure on 12 October 1986. Numbers represent data 
runs, short dotted lines are MWL's, and distance alongshore is arbi- 
trary. 
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SUPERDUCK 
15 Oct 86 
Profile 255 
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Figure 37. Mean current structure on 15 October 1986. Numbers represent data 
runs, short dotted lines are MWL's, and distance alongshore is arbi- 
trary. 
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SUPERDUCK 
17 Oct 86 
Profile 195 




Figure 38. Mean current structure on 17 October 1986. Numbers represent data 
runs, short dotted lines are MWL's, and distance alongshore is arbi- 
trary. 



133 



SUPERDUCK 
18 Oct 85 
Profile 195 




Figure 39. Mean current structure on 18 October 1986. Numbers represent data 
runs, short dotted lines are MWL's, and distance alongshore is arbi- 
trary. 
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SUPERDUCK 
21 Oct 86 
Profile T70 




Figure 40. Mean current structure on 21 October 19.86. Numbers represent data 
runs, short dotted lines are'MWL's, and distance alongshore is arbi- 
trary. 
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Figure 41. Offshore 5 >x spectra and 5,.J from 'puv' gage 621 for 15 October 
1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Figure 42. Offshore S >x spectra and S } J from 'puv' gage 621 for 17 October 
1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Figure 43. 



Offshore S yx spectra and S y J from 'puv' gage 621 for 18 October 
1986. Numbers indicate data runs. Intervals are 95% confidence. 
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Figure 44. Wind speed, wave height, atmospheric stability, and initial C d data for 
October 11-12, 1986. Beach is assumed parallel to the y-a\is. Open 
arrowheads denote wind direction and closed arrowheads denote wave 
incident angle. Winds are adjusted to z = 10 m. 
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Figure 45. Wind speed, wave height, atmospheric stability, and initial C' d data for 
October 17-18. 1986. Beach is assumed parallel to the y-axis. Open 
arrowheads denote wind direction and closed arrowheads denote wave 
incident angle. Winds are adjusted to z = 10 m. 
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Figure 46. Wind speed, wave height, atmospheric stability, and initial C, data for 
October 21, 1986. Beach is assumed parallel to the y-axis. Open 
arrowheads denote wind direction and closed arrowheads denote wave 
incident angle. Winds are adjusted to z = 10 m. 
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the solid line and empirical breaking wave distributions given by 



equations (127_) and (1-10) are the dashed and dotted lines, respectively. 
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Figure 48. Wave height and breaking wave height distributions (hatched area) nor- 
malized to offshore H rms for 17 October. The Rayleigh distribution is 
the solid line and empirical breaking wave distributions given by 
equations (127) and (140) are the dashed and dotted lines, respectively. 
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the solid line and empirical breaking wave distributions given by 



equations (127) and (140) are the dashed and dotted lines, respectively. 
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Figure 50. Model-generated H rms and V compared with field observations versus 
distance for invariant Cj and B and including turbulent momentum ex- 
change, 15 October. Model results are solid lines and field observa- 
tions are "x"'s. Elevation is referenced to MSL. 
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change, 17 October. Model results are solid lines and field observa- 
tions are "x"'s. Elevation is referenced to MSL. 
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Figure 52. Model-generated H rms and V compared with field observations versus 
distance for invariant Cj- and B and including turbulent momentum ex- 
change, 18 October. Model results are solid lines and field observa- 
tions are "x"'s. Elevation is referenced to MSL. Model result prior to 
the last Cf iteration is represented by the dashed line. 
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Figure 53. Model-generated V compared to field observations versus offshore dis- 
tance for constant Cy and no turbulent momentum exchange, 17-18 Oc- 
tober. Model results are solid line and field observations are "x"'s. 
Error is mean least square error. 
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APPENDIX B. 



CURRENT VECTORS OVERLAYED ON BATHYMETRY 
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Figure 54. Sled and alongshore array mean current vectors for 1 1 October 1986 
overlaved on bathymetry. Velocity scale is at upper right. 
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Figure 55. Sled and alongshore array mean current vectors for 12 October 1986 
overlayed on bathymetry. Velocity scale is at upper right. 
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Figure 56. Sled and alongshore array mean current vectors for 15 October 1986 
overlayed on bathymetry. Velocity scale is at upper right. 
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Figure 57. Sled and alongshore array mean current vectors for 16 October 1986 
overlayed on bathymetry. Velocity scale is at upper right. 



153 




Figure 58. Sled and alongshore array mean current vectors for 17 October 1986 
overlayed on bathymetry. Velocity scale is at upper right. 
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Figure 59. Sled and alongshore array mean current vectors for 18 October 1986 
overlaved on bathymetry. Velocity scale is at upper right. 
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Figure 60. Sled and alongshore array mean current vectors for 21 October 1986 
overlayed on bathymetry. Velocity scale is at upper right. 



156 



APPENDIX C. ATMOSPHERIC STABILITY DETERMINATION 



The following set of equations represent an iterative procedure to 
based on the stability flux profile relation of Businger et al. (1971). 

Initial equations: 

T dew =T air + 35( ln(.0 1 /?//)) 

AT= T air - T sur f ace — .009 SZ T 
T Kel =T air + 273.15 



Assign all — s = 0 
Determine V from iteration: 



T = 



AT k 






T* = TC 



x = (l- 16 (-^))T 



T = 2 ln[ 



( 1 + -v ) 



] 



V = -7(f) 



T<° 



L 



>0 



Q = Q s (T dew + 273.\5)(.0\RH) 



Qo = Q s (T sea + 273.15) 

Q s = .62510 " k = 23.84-^^-- 5.03 log(7\ v/ ) 

1 K el 




((?-{?> 

^oq L * 



Tl = T +6. lx \0~\T air Q) 



determine z L 

(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 
(C 7) 

(C8) 

(C9) 

(CIO) 
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II 

H f/ * 

N 


(Cll) 


( z ) * K2 T T * 
{ 1 )T T Kel ul 


(Cl 2) 


, 2 S SKZ q T v 

[ L q T u 2 

1 Ke! u * 


(C 13) 


Continue iterations from equation (C4) until ( ) converges. 





Definitions: 



z ? 

RH = 
T d - 

Z 0T = 

r 0t? = 

v* = 

^ = 
Q = 



temperature measurement height (m), 
humidity measurement height (m), 
relative humidity (%), 
dewpoint temperature (degrees C) 
temperature roughness length (2 x 10~ 6 ) 
humidity roughness length (2 x 10~ 6 ) 
integral diabatic term for temperature 
integral diabatic term for humidity 
specific humidity (g kg) 
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APPENDIX D. SURFACE WEATHER MAPS 
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APPENDIX E. ANALYSIS PARAMETERS FOR ALL SLED DATA RUNS 



Run 


start 

time 

(hhmniss) 


h 

(m) 


I_7 

1 1 rms 

(m) 


fm 

(Hz) 


V 


EMW 

(deg) 


^ anem 

(m) 


tan p 


11-1 


123504 


2.99 


0.92 


0.0938 


24.33 


1.4 


5.92 


.000 


11-2 


134359 


2.80 


0.S5 


0.0938 


IS. 86 


1.2 


6.11 


.000 


11-3 


143700 


1.95 


0.76 


0.093S 


18.51 


1.5 


6.95 


.106 


12-1 


074600 


2.97 


1.12 


0.0938 


-134.39 


-5.1 


5.95 


.059 


12-2 


092453 


2.03 


0.73 


0.093S 


13.23 


1.2 


6.88 


.000 


12-3 


100300 


2.29 


0.64 


0.0781 


-14.03 


-1.6 


6.64 


.000 


12-5 


145000 


2.93 


1.21 


0.07S1 


-56.10 


-1.8 


6.00 


.000 


12-6 


155700 


3.03 


0.93 


0.0781 


-28.95 


-1.6 


5.89 


.000 


12-7 


164200 


1.94 


0.82 


0.07S1 


13.61 


0.9 


6.97 


.106 


15-1 


103257 


2.79 


0.92 


0.1563 


150.03 


9.1 


6.13 


.029 


15-2 


120000 


1.78 


0.76 


0.1563 


60.60 


5.2 


7.15 


.029 


15-3 


124700 


1.67 


0.57 


0.1563 


16.76 


2.5 


7.25 


.029 


15-4 


133100 


2.30 


0.59 


0.203 1 


33.45 


5.0 


6.62 


.053 


16-2 


122200 


3.33 


0.98 


0.1875 


250.54 


14.7 


5.60 


.012 


16-3 


131S21 


1.S9 


O.SO 


0.1875 


151.02 


12.1 


7.03 


.034 


16-4 


140600 


1.51 


0.57 


0.1 S75 


57.20 


8.9 


7.41 


.034 


16-5 


144900 


1.54 


0.52 


0.1875 


47.99 


9.0 


7.39 


.015 


16-6 


153600 


1.91 


0.49 


0.1875 


25.94 


5.4 


7.01 


.005 


17-1 


092040 


3.85 


0.72 


0.1719 


61.53 


6.3 


5. OS 


.017 


17-2 


100036 


2.59 


0.69 


0.1719 


52.S4 


5.7 


6.32 


.034 


17-3 


104900 


1.43 


0.57 


0.1094 


24.14 


3.4 


7.50 


.034 


17-4 


113500 


1.18 


0.41 


0.1094 


27.96 


7.8 


7.74 


.034 


17-5 


122S57 


1.09 


0.36 


0.1094 


70.37 


31.4 


7. S3 


.034 


17-6 


134336 


1.41 


0.24 


0.0938 


1.45 


1.1 


7.52 


.017 


17-7 


141749 


1.55 


0.26 


0.1094 


3.99 


2.9 


7.38 


.017 


17-8 


150024 


1.43 


0.32 


0.1094 


2.60 


1.2 


7.45 


.017 
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18-1 


110730 


3.80 


0.93 


0.2031 


175.62 


11.8 


5.12 


.007 


18-2 


115000 


2.38 


0.84 


0.1250 


124.59 


8.6 


6.54 


.032 


18-3 


122555 


1.49 


0.67 


0.1094 


49.62 


5.3 


7.44 


.032 


18-4 


130900 


1.13 


0.38 


0.1094 


99.69 


45.0 


7.80 


.007 


18-5 


134937 


1.48 


0.28 


0.0469 


3.71 


2.2 


7.45 


-.025 


21-1 


110319 


3.55 


0.76 


0.0938 


-6.82 


-0.6 


5.31 


.021 


21-2 


121953 


1.65 


0.68 


0.0938 


2.3S 


2.4 


7.22 


.000 


21-3 


130448 


2.51 


0.49 


0.0938 


-17.72 


-3.4 


6.35 


-.029 


21-4 


134S00 


2.33 


0.45 


0.0938 


2.08 


0.5 


6.54 


.014 


21-5 


144000 


2.17 


0.38 


0.0938 


11.54 


3.7 


6.71 


.014 


21-6 


154200 


3.05 


0.76 


0.0938 


4.20 


0.3 


5.82 


.021 


21-7 


163059 


1.20 


0.59 


0.0938 


0.28 


0.0 


7.64 


.010 


21-S 


171000 


2.22 


0.45 


0.0938 


6.13 


1.5 


6.65 


-.029 


21-9 


174900 


1.51 


0.48 


0.0781 


- 7.72 


-1.6 


7.35 


.089 



• start time is referenced to Eastern Standard Time 

• H rms is measured from variance of the pressure sensor 

• 2 anem is the height of the anemometer above the mean water level as measured 
by the pressure sensor for that run 

• x,y coordinates are referenced to the FRF coordinate system 

• f m is the peak freq in the sea sfc elevation spectrum 

• 'EMW' is the "equivalent monochromatic wave" incident angle (Thornton and 
Guza, 19S6) 

• S yx r is the mean S yx r for the submerged current meters on that run 

• f m is the peak frequency of the sea surface elevation spectra 
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Run 


sled x 
coord 
(m) 


sled v 
coord 
(m) 


Pa tm 

(10 z Pa) 


0 

(° 77V) 


(ms- 1 ) 


boa 

(deg) 


# sub- 
merged 
c/ Ill'S 


11-1 


168.44 


899.90 


1023.6 


067. 


13.37 


0.5 




11-2 


139.66 


905.65 


1023.8 


066. 


12.84 


0.5 




11-3 


116.11 


907.30 


1023.9 


065. 


12.72 


0.5 


'S 


12-1 


247.15 


900.44 


1024.1 


049. 


10.00 


0.5 


3 


12-2 


198.94 


904.57 


1023. S 


047. 


S.S2 


0.5 




12-3 


154.69 


904.92 


1023.7 


047. 


S.71 


0.5 


A 


12-5 


198.55 


904.42 


1023.6 


052. 


5.66 


0.5 


3 


12-6 


159.41 


902.77 


1023.5 


047. 


5.04 


0.5 


-> 


12-7 


125.21 


904. 1 7 


1023.6 


041. 


5.43 


0.5 


3 


15-1 


244.01 


SS7.14 


1017.3 


037. 


9.37 


4.0 




15-2 


216.02 


S91.67 


1017.2 


037. 


8.11 


4.0 


-> 


15-3 


191.78 


895.64 


1017.2 


039. 


7.81 


4.0 




15-4 


147.13 


901. S9 


1017.3 


040. 


7.80 


4.0 


-s 


16-2 


269.28 


1157.39 


1016.8 


030. 


12.28 


-0.5 




16-3 


218.76 


1162.25 


1016.9 


027. 


11.43 


-0.5 


-S 


16-4 


201.06 


1165.24 


1017.0 


025. 


10.90 


-0.5 




16-5 


1SS.07 


1 166.64 


1017.1 


026. 


10.22 


-0.5 


2 


16-6 


154.58 


1 170.37 


1017.3 


026. 


10.04 


-0.5 


3 


17-1 


264.71 


1 1 79. S3 


1019.1 


001. 


6.23 


-0.5 




17-2 


230.10 


1 1S2.31 


1019.2 


276. 


7.24 


-0.5 




17-3 


203.50 


1184.07 


1019.3 


150. 


7.06 


-0.5 


2 


17-4 


193.76 


1184.49 


1019.4 


037. 


6.55 


-n.5 


i 


17-5 


1S9.67 


1184.39 


1019.5 


033. 


5.89 


-0.5 


i 


17-6 


15S.10 


11SS.54 


1019.9 


025. 


5.53 


-0.5 


2 


17-7 


153.23 


1 IS 9. 60 


1020.1 


026. 


5. SO 


-0.5 


2 


17-8 


146.96 


1190.44 


1020.2 


029. 


5.80 


-0.5 


2 


18-1 


275.56 


1179.98 


1026. S 


023. 


inop 


1.5 


3 


18-2 


227.23 


1178.08 


1027.1 


021. 


inop 


1.5 


3 


18-3 


20S.39 


1177.67 


1027.2 


01S. 


12.13 


1.5 


2 


18-4 


1S6.79 


1179.03 


1027.4 


017. 


12.59 


1.5 


1 


18-5 


159.11 


1 185.79 


1027.4 


018. 


13.75 


1.5 


1 
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21-1 


252.S9 


1282.94 


1023.1 


353. 


in op 


0.8 


-■> 


21-2 


204.39 


1 280. SO 


1022.5 


044. 


inop 


0.8 




21-3 


182.12 


1283.45 


1022.6 


035. 


inop 


0.8 


J) 


21-4 


173.04 


1283. SI 


1022.7 


041. 


inop 


0.8 


-> 


21-5 


157.84 


1285.67 


1022.7 


056. 


inop 


0.8 




21-6 


259.15 


1295.27 


1022.8 


083. 


inop 


0.8 


3 


21-7 


207. SS 


1292.45 


1022.8 


089. 


in op 


0.8 


i 


21-8 


180.63 


1291.45 


1022.8 


079. 


inop 


0.8 


2 


21-9 


147.96 


1290.75 


1022.9 


083. 


inop 


0.8 


2 



• sled x.y coordinates are referenced to the FRF coordinate system 

• 'boa' is beach orientation angle 

• MRS conversion for pressure: 1 mb = 100 Pa 



Run 


kfj 

n 


H 0 

(m) 


fn 

of 

S,,v 

(Hz) 


Q 


Hfms 

(m) 


ill 


-1.2 




.0156 






112 


-1.7 




.0781 






1 13 


1.0 




.0781 






121 


-4.3 




.0781 






122 


-0.4 




.1250 






123 


-4.3 




.078 1 






125 


-1.2 




.0781 






126 


0.0 




.0781 






127 


1.8 




.0625 






151 


10.2 


1.43 


.'1563 


0.00 


0.79 


152 


6.5 


1.40 


.1719 


0.44 


0.73 


153 


5.8 


1.30 


.1563 


0.11 


0.56 


154 


8.9 


1.14 


.1719 


0.04 


0.51 


162 


19.0 


1.59 


.1875 


0.00 


0.78 


163 


14.9 


1.53 


.1375 


0.13 


0.79 



167 



164 


12.6 


1.46 


. 1 S 75 


0.06 


0.55 


165 


IS. 6 


1.41 


. 1 S75 


6.64 


0.56 


1 66 


9.5 


1.35 


. 1 S75 


0.02 


< >.45 


171 


7 .6 


1.16 


.1 S~5 


(>.00 


0.56 


172 


s.s 


1.14 


.1875 


0.61 


o.59 


173 


-2.S 


1.1 1 


.1719 


0.19 


0.57 


174 


1.8 


1.15 


.1719 


0.19 


0.42 


175 


-2.1 


1 .09 


. 1 250 


0.21 


0.38 


176 


0.2 


1.09 


.0938 


0.17 


0.24 


177 


-1.9 


1.11 


.07S1 


0.06 


0.23 


178 


0.2 


1.11 


.0781 


6,07 


0.28 


181 




1.44 


.1875 


0.00 


0.73 


182 


-1.6 


1.46 


.18'5 


0.00 


0.73 


183 


-1.4 


1.46 


.1719 


1.00 


0.68 


184 


14.1 


1.4S 


.1094 


1 .00 


0.39 


1S5 


-0.7 


1.67 


.2506 


0.00 


0.27 


211 


-3.6 




.0938 






212 


4.5 




.09 3 S 






213 


-3.6 




.O'Sl 






214 -1.2 




.0156 






215 | 6.7 




.0625 






216 | -2.0 




.10^4 






217 


^ . 




.6938 






218 | 6.6 


i 


.0935 






219 | -1.6 




.6156 






• is the mean a of the rear and front current 
meter (if they were both submerged) and is ob- 
tained from the significant wave angle spectra for 
peak frequency in sea surface elevation spectra 
(see equation 71 ). 

• H 0 is determined from Sl'PERDUCK gage 630. 

• of S vx is the peak frequency from the rear cur- 
rent meter S vx spectra 

• H rms is determined from the ?/ time series using 
equation 13S 

• Q is the fraction of all waves which are breaking. 
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APPENDIX F. RADIATION STRESS, SEA SURFACE ELEVATION, 
SIGNIFICANT WAVE ANGLE, AND CURRENT VELOCITY SPECTRA 
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Figure 65. Radiation stress (S IX ), sea surface elevation ( 7 ), significant nave angle 
(a), and current velocity spectra for sled run 15-1. 
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Figure 66. Radiation stress (5^), sea surface elevation (rj), significant nave angle 
(a), and current velocity spectra for sled run 15-2. 
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Figure 67. Radiation stress (£,.*), sea surface elevation (tj), significant wave angle 
(a), and current velocity spectra for sled run 15-3. 
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Figure 68. Radiation stress sea surface elevation (tj), significant wave angle 

(z), and current velocity spectra for sled run 15-4. 
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RUN 16-2 — SUPlRDUCK 
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Figure 69. Radiation stress (S x ), sea surface elevation (tj), significant wave angle 
(a), and current velocity spectra for sled run 16-2. 
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Figure 70. Radiation stress (S >x ), sea surface elevation (rj), significant wave angle 
(a), and current velocity spectra for sled run 16-3. 
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Figure 71. Radiation stress (S.,*), sea surface elevation (r]), significant wave angle 
(a), and current velocity spectra for sled run 16-4. 
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RUN 16-5 — SUPFRDUCK 
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Figure 72. Radiation stress (5 >x ), sea surface elevation ( tj ), significant wave angle 
(a), and current velocity spectra for sled run 16-5. 
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Figure 73. Radiation stress (S yx ), sea surface elevation ( t ] ), significant wave angle 
(a), and current velocity spectra for sled run 16-6. 
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Figure 74. Radiation stress (5 >x ), sea surface elevation (»/), significant wave angle 
(a), and current velocity spectra for sled run 17-1. 
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RUN 17-2 — SUPERDUCK 
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Figure 75. Radiation stress (S yx ), sea surface elevation (>/), significant wave angle 
(a), and current velocity spectra for sled run 17-2. 
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Figure 76. Radiation stress (5 >x ), sea surface elevation (?/), significant wave angle 
(£), and current velocity spectra for sled run 17-3. 
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Figure 77. Radiation stress (5^), sea surface elevation (//), significant »ave angle 
(a), and current velocity spectra for sled run 17-4. 
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Figure 78. Radiation stress sea surface elevation (;/), significant wave angle 

(a), and current velocity spectra for sled run 17-5. 



183 



RUN 17-6 — SUPERDUCK 



95% x 2 (0.73,1.45) 

Syx OL 





Current Velocity 




Frequency (Hz) Frequency (Hz) 



Figure 79. Radiation stress (S }X ), sea surface elevation (;/), significant wave angle 
(a), and current velocity spectra for sled run 17-6. 
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Figure 80. Radiation stress (5 >JC ), sea surface elevation (//), significant \>ave angle 
(a), and current velocity spectra for sled run 17-7. 
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Figure 84. Radiation stress (5^), sea surface elevation (tj), significant wave angle 
(a), and current velocity spectra for sled run 18-3. 
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Figure 85. Radiation stress ( S x ), sea surface elevation (//), significant wave angle 
(a), and current velocity spectra for sled run 18-4. 
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Figure 86. Radiation stress (5 >x ), sea surface elevation (rj), significant wave angle 
(a), and current velocity spectra for sled run 18-5. 
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APPENDIX G. MEAN CURRENT VELOCITIES AND TOTAL 
RADIATION STRESS FOR EACH SLED CURRENT METER AND RUN 



Run 


C/M 1 (mid) 


C/M 2 (upper) 


C/M 4 (loner) 


Syl (•/'”--) 


V 

(ms- 1 ) 


U 

(nis- 1 ) 


V 

(ms- 1 ) 


U 

(ms- 1 ) 


V 

(ms -1 ) 


U 

(ms- 1 ) 


C/M 

1 


C/M 

4 


C/M 


111 


-0.41 S 


0.36" 


-0.512 


0.266 


-0.354 


0.407 








112 


-0.5 1 4 


0.195 


-0.5SS 


0.155 


-0.405 


0.154 








113 


-0.256 


0.0"’ 4 


-0.339 


0.130 


-0.164 


0.005 








121 


0.307 


0.252 


0.259 


0.325 


0.304 


0.173 


-134.40 


-134.40 


-134.40 


122 


0.165 


0.517 


0.166 


0.365 


0.069 


0.592 


S.8S 


13.01 


17.79 


123 


-0.246 


0.295 


-0.320 


0.31S 


-0.IS9 


0.193 


-21.60 


-8.81 


-11.67 


12? 


0.295 


0.476 


0.280 


0.471 


0.278 


0.449 


-66.46 


-39.30 


-62.55 


126 


0.219 | 0.323 


0.195 


0.2S7 


0.204 


0.270 


-25.27 


-34.55 


-27.02 


127 


0.544 ) o.l IS 


0.49S 


0.123 


0.476 


0.069 


5.78 


25.20 


9.84 


151 


-0.393 


0.06S 


-< >.470 


0.108 


-0.396 


0.153 


1504 »0 


150.00 


150.00 


152 


-0.7S6 


0.24" 


-0.S30 


0.141 


-0.798 


0.385 


59.47 


43.67 


7S.S7 


153 


-0.773 


0.207 


-0.SS2 


-0.045 


-0.723 


0.425 


16.84 


16.67 


45.94 


154 


-0.719 


0.053 


-0.S24 


0.073 


-0.669 


0.130 


36.62 


22.32 


41.41 


162 


-0.443 


0.01 1 


-0.532 


0.052 


-0.437 


0.074 


250.50 


250.50 


250.50 


163 


-0.9S5 


0.040 


-1.090 


0.050 


-0.993 


0.076 


153.50 


143.10 


156.40 


164 


-1.209 


-0.039 


-0.931 


-0.320 


-1.156 


0.111 


62.10 


52.29 


* 


165 


-1.131 


-0.0S9 


-0.S92 


-0.309 


-1.072 


0.036 


4S.40 


47.5S 


* 


166 


-0.63S 


-O.O0S 


-0.747 


-0.073 


-0.605 


0.070 


27.52 


25.64 


24.66 


171 


.0.219 


-0.024 


-0.301 


0.025 


-0.191 


0.054 


61.53 


61.53 


61.53 


172 


-0.270 


0.015 


-0.355 


0.078 


-0.254 


0.086 


53.75 


53.12 


51.66 
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173 


- 0.414 


- 0.037 


- 0.280 


- 0.291 


- 0.382 


0.074 


25.65 


22.62 


❖ 


174 


- 0.327 


- 0.335 


- 0.144 


- 0.144 


- 0.299 


- 0.022 




15.34 


* 


175 


- 0.299 


- 0.485 


- 0.114 


- 0.041 


- 0.376 


- 0.150 


* 


10.67 


* 


176 


0.003 


- 0.163 


- 0.106 


- 0.144 


0.023 


0.028 


0.95 


1.95 


# 


177 


0.031 


- 0.074 


- 0.078 


- 0.117 


0.047 


0.055 


4.03 


3.95 


❖ 


178 


0.038 


- 0.045 


- 0.091 


- 0.133 


0 . 02 S 


0.065 


3.56 


1.63 


* 


181 


- 0.402 


0.010 


- 0.479 


0.060 


- 0.376 


0.077 


175.60 


175.60 


175.60 


182 


- 0.703 


0.050 


- 0.795 


0.125 


- 0.693 


0.106 


122.00 


122.60 


129.10 


183 


- 0.981 


0.126 


- 0.708 


- 0.252 


- 0.978 


0.227 


58.04 


41.20 




184 


- 0.778 


- 0.290 


- 0.198 


- 0.092 


- 0.896 


0.073 


* 


14.04 


* 


185 


- 0.622 


- 0.022 


- 0.419 


- 0.181 


- 0.554 


0.123 


3.02 


4.39 


* 


186 


- 0.651 


0.054 


- 0.748 


0.075 


- 0.656 


0.118 








211 


- 0.023 


0.033 


- 0.086 


0.107 


- 0.012 


0.119 


- 6.82 


- 6.82 


- 6.82 


212 


- 0.002 


0.209 


- 0.071 


- 0.063 


0.010 


0.331 


20.18 


27.38 


* 


213 


- 0.120 


0.240 


- 0.177 


0.218 


- 0.098 


0.343 


- 22.19 


- 16.32 


- 14.64 


214 


- 0.256 


0.112 


- 0.313 


0.096 


- 0.276 


0.252 


0.83 


2.71 


2.67 


215 


- 0.183 


0.106 


- 0.285 


0.176 


- 0.121 


0.190 


11.20 


13.86 


9.56 


216 


0.003 


- 0.053 


- 0.068 


- 0.007 


0.005 


0.064 


4.20 


4.20 


4.20 


217 


0.060 


- 0.263 


- 0.082 


- 0.228 


0 . 0 S 3 


- 0.036 


❖ 


3.44 


❖ 


218 


- 0.150 


- 0.037 


- 0.204 


- 0.140 


- 0.156 


0.171 


7.63 


3.83 


6.93 


219 


- 0.125 


0.115 


- 0.128 


- 0.042 


- 0.087 


0.173 


- 11.35 


- 4.10 


1.84 



* indicates current meter was not underwater for entire data run 
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